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“GEOPHYSICAL  TESTS  FOR 
INTERMEDIATE-RANGE  FORCES” 


The  suggestion  in  1986  of  a  possible  “fifth  force”  in  Nature  led  to  consid¬ 
erable  experimental  and  theoretical  effort  to  detect  such  a  force.  Evidence 
for  the  presence  of  such  a  force  could  come  from  apparent  deviations  from 
the  predictions  of  Newtonian  gravity.  We  can  characterize  this  new  force  by 
writing  the  effective  gravitational  potential  between  point  masses  mx  and  m2 
in  the  form. 


V(r)  = 


-Gmim2 


(1  +  ae-r/x)  =  Visr(r)  4-  Vs(r) 


Here  r  —  |fi  —  r2|  is  the  separation  of  the  masses,  A  =  h/mc  is  the  Comp¬ 
ton  wavelength  of  the  particle  responsible  for  the  new  interaction,  and  a 
is  a  dimensionless  constant  which  characterizes  the  overall  strength  of  this 
interaction  relative  to  gravity.  The  results  of  any  experimental  search  for  non- 
Newtonian  gravity  can  thus  be  expressed  in  terms  of  the  strength  parameter 
a  and  the  range  A.  For  technical  reasons  a  given  experiment  can  set  sensitive 
limits  only  on  values  of  A  comparable  to  the  natural  size  of  the  experimental 
system.  For  distance  scales  of  order  10m  -  10  km  geophysical  experiments 
using  towers,  boreholes,  etc.  provide  the  natural  distance  scale,  and  this  ob¬ 
servation  has  motivated  much  of  the  recent  interest  in  non-Newtonian  gravity. 
At  the  same  time,  the  fundamental  unity  of  physics  demands  that  geophysical 
experiments  be  integrated  along  with  other  tests  for  non-Newtonian  physics. 
Thus  our  efforts  have  been  directed  both  at  geophysical  tests  as  well  as  at 
other  tests  for  evidence  of  non-Newtonian  gravity. 


1 


Since  the  work  done  here  has  been  published  in  various  journals  and 
conference  proceedings,  we  will  simply  summarize  the  topics  covered  in  this 
Introduction,  and  refer  to  the  full  papers  in  the  Appendices  for  complete 
details. 

2.  Geophysical  Tests  for  Non-Newtonian  Gravity 

A.  Finite-Size  Effects  in  Eotvos- Type  Experiments  [E.  Fischbach  and 
C.  Talmadge  in  “New  and  Exotic  Phenomena’  ,  Proceedings  of  the  XXVth 
Rencontre  de  Moriond,  Les  Arcs,  France,  20-27  January  1990,  edited  by 
O.  Fackler  and  J.  Tran  Thanh  Van  (Editions  Frontieres,  Gif-Sur- Yvette, 
1990)  pp.  187-196].  This  paper  investigates  the  possibility  that  all  non- 
Newtonian  effects  are  fundamentally  composition-independent,  but  that 
composition-dependent  effects  may  arise  due  to  the  fact  that  V2Vs(r)  /  0 
for  V5(r)  given  in  Eq.(l).  This  paper  is  reprinted  Appendix  A. 

B.  Is  the  Eotvos  Experiment  Sensitive  to  Spin?  [A.M.  Hall,  H.  Arm- 
bruster,  E.  Fischbach  and  C.  Talmadge,  in  Progress  in  High  Energy  Physics, 
Proceedings  of  the  Second  International  Conference  on  Medium  and  High 
Energy  Nuclear  Physics,  Taipei,  Taiwan,  ROC,  8-18  May,  1990,  edited  by 
W.-Y.P.  Hwang,  S.-C.  Lee,  C.-E.  Lee,  and  D.J.  Ernst,  (North-Holland  Else¬ 
vier,  New  York,  1991),  pp.  325-339.) 

This  paper  discusses  the  possibility  that  the  systematic  effect  uncovered 
in  1986  by  Fischbach  et  al.  in  the  Eotvos  data  could  represent  the  coupling 
of  a  new  field  to  a  novel  spin-dependent  change.  This  model  has  specific  im¬ 
plications  for  geophysical  experiments.  This  paper  is  reprinted  in  Appendix 

B. 

C.  Exponential  Models  of  non-Newtonian  Gravity  [E.  Fischbach,  C. 


2 


Talmadge,  and  D.  Krause,  in  Phys.  Rev.  D43,  460-467  (15  January,  1991).] 
The  non-Newtonian  potential  in  Eq.(l)  has  the  typical  Yukawa  form, 
which  gives  the  spatial  variation  (l//r)e-r/\  However,  in  some  models  the 
spatial  variation  of  the  potential  is  given  by  an  experimental  ~  ( l/\)e~rlx . 
This  changes  the  phenomenology  of  the  putative  fifth  force  in  a  significant 
way,  which  we  describe  in  detail  in  Appendix  C. 

D.  Non-Newtonian  Gravity  and  New  Weak  Forces:  An  Index  of 
Measurements  and  Theory  [E.  Fischbach,  G.T.  Gillies,  D.E.  Krause,  J.G. 
Schwan,  and  C.  Talmadge,  in  Metrologia  29,  215-260  (1992).] 

This  paper  presents  a  short  overview  and  review  of  the  fifth  force  as  an 
introduction  to  the  first  ever  bibliography  on  the  subject  of  new  weak  forces. 
The  bibliography  contains  813  entries  by  825  authors,  and  is  presented  in 
Appendix  D. 

E.  One  Hundred  Years  of  the  Eotvos  Experiment  [L.Bod,  E.  Fischbach, 
G.  Marx,  and  M.  Naray-Ziegler,  in  Acta  Physica  Hungarica  69,  335-355 
(1991).] 

The  original  Eotvos  experiment  provided  much  of  the  motivation  for  the 
current  interest  in  searches  for  non-Newtonian  gravity.  In  this  paper  we 
provide  additional  historical  details  on  this  experiment,  which  remains  one 
of  the  few  which  gives  any  evidence  for  the  possibility  of  a  new  force.  This 
paper  constitutes  Appendix  E. 

F.  Six  Years  of  the  Fifth  Force  [E.  Fischbach  and  C.  Talmadge,  in  Na¬ 
ture  356,  207-215  (1992).] 

This  review  was  solicited  by  the  Editor  of  Nature,  and  summarizes  the 
status  of  searches  for  the  putative  “fifth  force”  as  of  1992.  It  presents  the 
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exclusion  plots  which  summarize  what  had  been  learned  up  to  that  point. 
This  paper  is  presented  in  Appendix  F. 

G.  The  Second  Coming  of  Tower  Gravity:  An  Update  [D.H.  Eckhardt, 
A.J.  Romaides,  R.W.  Sands,  C.  Jekeli,  E.  Fischbach,  C.L.  Talmadge  and  H. 
Kloor,  in  “The  Second  Coming  of  Tower  Gravity:  An  Update”,  in  Per¬ 
spectives  in  Neutrinos,  Atomic  Physics  and  Gravitation.  Proceedings  of  ihe 
XHIth  Moriond  Workshop  (Villars,  Switzerland,  30  January  -  6  February 
1993)  edited  by  J.  Tran  Thanh  Van,  (Editions  Frontieres,  to  be  published.] 

The  pioneering  work  of  the  Air  Force  Geophysics  Laboratory  on  gravity 
measurements  carried  out  on  towers,  has  been  pursued  by  a  new  tower  ex¬ 
periment  using  the  WABG  tower  in  Inverness,  Mississippi.  The  first  results 
from  this  experiment  are  given  in  Appendix  G. 

H.  Global  Positioning  Systems  and  Television  Signals:  Are  the 
Two  Compatible  [A.J.  Romaides,  R.W.  Sands,  E.  Fischbach  and  C.L.  Tal¬ 
madge,  submitted  to  Surveying  and  Land  Information  Systems .] 

In  the  course  of  carrying  out  the  tower  experiment,  we  noticed  that  the 
measurement  of  ground  elevations  using  the  Global  Positioning  Satellite  Sys¬ 
tems  could  not  always  be  carried  out  in  close  proximity  to  active  television 
towers.  In  this  paper  we  discuss  this  problem  in  detail,  and  describe  our 
efforts  to  understand  its  origins.  This  paper  is  given  in  Appendix  H. 
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APPENDIX  A 


PURD-TH-90-1 


FINITE-SIZE  EFFECTS  IN  EOTVOS-TYPE  EXPERIMENTS* 

Ephraim  Fischbach  and  Carrick  Talmadge 

Physics  Department,  Purdue  University,  West  Lafayette,  IN  47907  USA 


in  New  and  Exotic  Phenomena 
Proceedings  of  the  Xth  Moriond  Workshop 
Los  Arcs,  France,  20-27  January,  1990 


ABSTRACT 

We  discuss  a  new  class  of  experiments  in  which  searches  for  composition-dependent  de¬ 
viations  from  Newtonian  gravity  can  be  used  to  set  limits  on  composition-independent  effects 
as  well.  These  experiments  utilize  the  observation  that  objects  with  different  charge  or  mass 
distributions,  will  couple  differently  to  the  (non- vanishing)  Laplacian  of  a  non-Newtonian 
potential. 


*  Work  supported  in  part  by  the  the  United  States  Air  Force  Geophysics  Laboratory  and 
by  the  United  States  Department  of  Energy 
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I.  INTRODUCTION 


Current  searches  for  deviations  from  the  predictions  of  Newtonian  gravity,1  _B  have  thus 
far  focussed  primarily  on  two  classes  of  experiments:  These  are  the  composition-independent 
experiments, which  look  specifically  for  a  breakdown  of  the  1/r2  force  law,  and  the  com¬ 
position-dependent  experiments  which  compare  the  accelerations  of  chemically-differcnt  test 
masses  towards  a  common  source.  Doth  classes  of  experiments  can  be  understood  in  terms 
of  the  conventional  parametrization  of  the  potential  energy  V (r)  for  two  point  masses  t  and 
j,  in  which  a  “fifth  force”  contribution  Vs(r)  is  added  to  the  Newtonian  term  VV(r)  to  give3 

V(r)  =  VN(r)  +  Vs(r)  =  ±  f*QQ±e~rl\  (1) 

r  r 


Here  Goo  is  the  Newtonian  gravitational  constant,  mi(mj)  is  the  mass  of  i(j),  X  is  the  range 
of  the  fifth  force,  and  Qi(Qj)  is  the  corresponding  charge  which  determines  the  fifth  force 
coupling.  The  overall  strength  is  set  by  the  constant  /,  which  is  the  fifth  force  analog  of  the 
electric  charge  e,  and  the  ±  sign  reflects  the  possibility  that  the  contribution  from  V5  can  be 
attractive  or  repulsive.  Using  (1)  we  can  write  V(r)  in  the  form3 


V(r)  =  -G001^-(l±aije- 
r 


Oti 


V 


-TO* 


r/A), 


(2) 


Goo  ml 


H 


In  Eq.(2)  p,  =  m./m#,  where  mg  =  and  £  gives  the  strength  of  V5  relative  to 

-* 

gravity.  Differentiating  Eq.(2)  we  can  write  the  force  F(r)  in  the  form 


F(r)  =  -  W(r)  = 


—G(r)mimjf 

r2 


G(r)  =  Goo[l  ±  «.»(!  +  r/X)e~r'x). 


(3) 


The  deviation  of  G(r)  from  G <»  is  an  indication  of  the  presence  of  a  non-Newtonian  coupling, 
and  we  sec  from  Eq.(3)  that  such  deviations  can  be  detected  in  two  ways:  A  dependence  of 
the  effective  strength  G(r)  on  the  separation  r  of  the  test  masses  would  reflect  a  breakdown 
of  the  expected  1/r2  law  of  Newtonian  gravity.  Secondly  the  constant  a,_,  in  (3)  is  in  general 
different  for  each  pair  i,j  of  materials,  and  hence  the  effective  strength  G(r)  depends  on 
the  compositions  of  the  interacting  test  masses.  Experiments  which  look  for  a  dependence  of 
G(r)  on  r,  and  on  a,j,  are  denoted  respectively  as  composition-independent,  and  composition- 
dependent,  searches  for  non-Newtonian  gravity. 

In  what  follows  we  point  to  the  existence  of  a  new  class  of  experiments  which  complements 
those  described  above,  and  which  are  distinguished  by  the  fact  that  they  depend  on  the  finite 
size  of  the  test  masses  relative  to  the  range  X.  The  principle  behind  these  experiments  was 
first  described  by  Stacey,  10  and  is  a  consequence  of  the  fact  that  V2V^  ^  0. 
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DERIVATION  OF  THE  FINITE-SIZE  EFFECT 
To  understand  the  origin  of  the  “finite-size  effect”  consider  the  interaction  energy  W  of 
a  test  object  whose  charge  distribution  is  p(f ),  with  the  potential  ♦(f )  of  an  external  field, 

W  =  J  d3rp(rMn ■  (4) 

♦(f)  can  be  expanded  about  the  center-of-mass  of  the  test  object  (f  =  0)  and,  retaining  the 
first  few  terms,  we  have 

W=  f  d3rp(f  )[«(0)  +  f  •  V*(0)  +  \  Y,  0)  +  •••]>  (5) 

^  ij 

where  d,  s  d/dx\  and  (f ),  =  x,-.  If  ♦(f)  arises  from  a  massless  field  (so  that  V2^(f )  = 
0),  then  the  usual  multipole  expansion  for  W  is  obtained  by  adding  to  Eq.(5)  the  term 
—  8  |f  I2^«jV2^(0),  which  allows  the  third  term  in  (5)  to  be  expressed  in  terms  of  the  quadru- 
pole  moment  tensor  Qij, 

Qij  =  d3rp(r)[2xiXj  -  |f|a$y].  (6) 

However,  in  searching  for  deviations  from  Newtonian  gravity  we  are  probing  specifically  for 
interactions  for  which  V2^(f )  ^  0.  It  follows  that  if  we  wish  to  expand  W  in  terms  of  the 
same  multipole  moments  that  arise  in  the  massless  case,  then  we  must  add  back  the  term 
proportional  to  V2^(0),  so  that  the  leading  rotationally-invariant  contribution  to  W  now  has 
the  form 

W  =  J d3rp(f)[  1  +  \\r  \2  V2  +  •  •  •]♦(()).  (7) 

Since  the  term  proportional  to  jf  |2  is  independent  of  the  orientation  of  the  test  object,  it  will 
behave  as  an  additional  contribution  to  its  mass.  To  see  what  effect  this  term  has  consider 
the  case  where  ♦(f)  is  given  by  a  Yukawa, 

$(r)  =  7T~~^~T exp(-|f -  f '  |/A),  (8) 

|r  —  r  | 

corresponding  to  a  point  source  of  mass  M  and  strength  aG  located  at  r' .  Since  ♦(f )  is  a 
solution  of  the  time-independent  Klein-Gordon  equation,  it  follows  that 

V2^(f )  =  (l/A*Wf),  (9) 

and  hence 

W  =  J  d3rp(r)[  1  +  +  •  •  j^(0)  =  m[l  +  +  •  •  •  -  (10) 
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Here  m  is  the  mass  of  the  test  object,  ( R 3)  =  (1/m)  f  d3rp(r)\f\2  is  its  mean-square-charge 
radius,  and  •  ■  *  denotes  the  remaining  terms  from  (5)  which  depend  on  higher  derivatives  of 
4(f).  We  see  from  (10)  that  for  a  Yukawa  potential  the  leading  correction  to  the  standard 
multipolc  formula,  arising  from  the  fact  that  X  is  finite,  has  the  effect  of  multiplying  the 
(inertial)  mass  m  by  the  expression  in  (  ]  in  (10).  As  we  now  show,  this  correction  leads  to 
an  apparent  violation  of  the  Equivalence  Principle. 

Let  Pi  denote  the  total  force  acting  on  a  test  mass  mi  in  the  combined  presence  of  the 
Earth’s  acceleration  field  g(r)  and  V$(f ), 

Pi  =  mj<7  -  mi(l  +  /<i)V$>,  (11) 


where  /cj  =  (R2)i/ 6A2.  The  acceleration  difference  of  objects  1  and  2  along  the  direction  n, 
A  a  —  (ai  -  a2)  ■  n,  is  then  given  by 


where  A/e  =  /cx  -  *2,  and  a  =  (at  +a2)/2.  It  follows  from  (12)  that  two  objects  with  different 
values  of  (R2)  will  experience  different  accelerations  in  the  presence  of  $,  irrespective  of 
whether  or  not  they  have  the  same  composition.  A/c  will  almost  always  be  different  from 
zero  in  an  Eotvos  experiment,  where  the  accelerations  of  two  objects  having  the  same  mass 
but  different  compositions  are  compared.  For  the  current  generation  of  Eotvos  experiments, 
where  test  samples  have  not  only  the  same  mass  but  the  same  external  dimensions  as  well, 
p(r)  must  be  different  for  the  two  test  masses,  and  hence  A/c  is  necessarily  different  from 
zero. 

The  finite-size  contribution  in  Eq.(7),  which  applies  to  any  matter  distribution  p(r)  in¬ 
teracting  with  an  arbitrary  potential  $(f ),  generalizes  a  result  originally  derived  by  Stacey3,10 
for  a  spherical  mass  in  the  field  of  a  point  Yukawa  source.  To  establish  the  connection  be¬ 
tween  (7)  and  Stacey’s  work,  we  consider  (as  he  did)  the  experiment  of  Thieberger11  in  which 
the  differential  acceleration  between  a  spherical  copper  shell  and  the  water  it  displaced  was 
measured.  For  a  uniform  sphere  of  radius  R  we  have 

(R2)  =  (13a) 


while  for  a  shell  with  outer  (inner)  radius  R2(R i)  the  corresponding  result  is 

/p2\  _  3  (R$  -  Rf) 

{“>-5(R*-R*y 


(136) 


From  Eqs.(13a,b)  and  (10)  we  find  immediately  that  the  (anomalous)  acceleration  difference 
between  the  shell  and  the  water  is  proportional  to 


A  -1$  (1  -R\IR%) 

A/C  -  Kgphere  Kshell  - 


(14) 
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which  is  what  Stacey  found.  In  addition  to  generalizing  Stacey’s  result,  Eq.(7)  also  demon¬ 
strates  that  the  leading  non-Newtonian  corrections  to  the  multipolc  formula  can  be  obtained 
directly ,  without  having  to  start  from  an  exact  analytic  result  as  Stacey  did.  This  means  that 
the  non-Newtonian  contribution  can  be  calculated  simply  for  test  masses  of  arbitrary  shape, 
where  an  analytic  expression  would  he  difficult  (if  not  impossible)  to  come  by. 

It  follows  from  the  preceding  discussion  that  any  experiment  which  measures  the  accel¬ 
eration  of  a  test  object  a,  will  be  sensitive  at  some  level  to  the  finite-size  anomaly  na.  This 
effect  will  be  the  dominant  signal  for  a  non-Newtonian  force  when  the  mean-square  charge 
radius  of  an  object  is  comparable  to  the  distance  scale  over  which  the  non-Newtonian  poten¬ 
tial  is  varying,  which  for  a  Yukawa  occurs  when  (R2)1^2  «  A.  This  observation  suggests  that 
the  finite-size  effect  may  offer  a  practical  means  of  adapting  Eotvos  experiments  to  carry  out 
high-precision  searches  for  composition-independent  non-Newtonian  gravity  over  distances  of 
several  centimeters,  which  is  the  characteristic  size  of  the  masses  that  are  typically  used.  More 
significantly,  by  appropriately  redesigning  the  test  masses,  it  may  be  possible  to  achieve  a  far 
greater  sensitivity  for  composition-independent  short-range  forces  than  is  currently  possible 
by  other  means. 

CONNECTION  TO  EXPERIMENT 

a)  Short-Range  Composition-Independent  Experiments 

At  present,  the  most  sensitive  composition-independent  experiments  at  short  ranges  are 
those  of  Hoskins,  et  al.12  (2  -  105  cm),  and  Chen  et  a/.13  (5-9  cm).  Hoskins  et  a/,  in  fact 
carried  out  two  separate  experiments  which  respectively  covered  the  intervals  2  — 5  cm  and  5- 
105  cm.  Although  various  clever  methods  were  used  to  enhance  any  possible  non-Newtonian 
signal,  none  of  these  measurements  is  a  null  experiment  in  the  same  sense  that  the  Eotvos 
experiment  is.  For  this  reason  an  Eotvos  experiment  utilizing  the  finite-size  effect  would 
possess  a  number  of  advantages  over  existing  experiments,  such  as  minimizing  the  sensitivity 
to  source  inhomogeneities,  and  to  the  source-detector  separation.  More  importantly,  since 
the  finite-size  effect  depends  on  the  Laplacian  of  the  non-Newtonian  force,  rather  than  on  the 
force  itself,  it  may  be  more  sensitive  than  existing  techniques  to  certain  types  of  couplings. 
For  the  special  case  of  a  point  Yukawa  source,  we  see  from  Eq.(9)  that  V2$(f )  has  the  same 
functional  dependence  on  f  as  does  $(f )  itself,  but  this  is  not  generally  the  case.  For  a  test 
mass  in  the  presence  of  an  arbitrary  force  f/(r)  we  have, 

V2(f/(r)]  =  -2 f  ^  +  2[V/(r)  •  V]f  +  fV2/(r),  (15) 

where  we  have  used  the  identity  V2f  =  —2 f/r2  in  the  first  term.  We  see  from  (15)  that  the 
functional  form  of  V2(f/]  can  be  quite  different  from  that  off  /  itself,  and  that  this  difference 
could  be  significant  for  small  r.  Consider  the  case  of  two  nearly  cancelling  Yukawas,  which 


9 


in  the  appropriate  limit  can  be  approximated  by  an  exponential  coupling,1 


$(r)  =  VE{r)  -  aE 


e-r/\ 


(16) 


FB(r)  =  -VVE(r)  =  raE- 


,-r/X 


(17) 


The  experiments  of  Hoskins,  et  al.,12  and  Chen,  et  al.,13  search  for  FE  directly,  and  hence 
are  sensitive  to  the  factor  exp(— r/A).  By  contrast,  the  finite-size  effect  is  determined  by 


( R 2) 
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V2FE(r)  =  f 


,(R2)  aBe-r'x 


6A2 


(18) 


and  hence  is  sensitive  to  a  different  radial  function.  It  is  easy  to  demonstrate  that  the  ex¬ 
periments  in  Refs.  (12  and  13)  would  have  very  limited  sensitivity  to  certain  combinations  of 
elementary  Yukawas,  whereas  the  finite-size  effect  would  be  quite  sensitive,  and  vice  versa. 
Thus  an  Eotvos  experiment  utilizing  the  finite-size  effect,  would  complement  existing  mea¬ 
surements  at  short  ranges  by  virtue  of  their  different  systematics.  The  same  observation 
applies  to  the  Laplacian  detector  of  Paik  et  al.,14  whose  greatest  sensitivity  will  be  to  forces 
of  somewhat  longer  range. 


b)  The  Eotvos  Experiment 

Since  (R2)  can  be  calculated  in  a  straightforward  way  for  any  test  mass,  it  is  natural 
to  apply  the  preceding  analysis  to  the  original  experiment  of  Eotvos,  Pekar,  and  Fekete 
(EPF)15  to  see  whether  there  is  any  suggestion  of  a  correlation  with  (R2)  in  their  data.  The 
torsion  balances  used  by  Eotvos,  Pek&r,  and  Fekete  (EPF)  in  performing  their  experiments 
were  originally  designed  to  measure  gravity  gradients,  which  was  the  purpose  of  the  set 
of  experiments  performed  by  Eotvos  in  the  mountains  of  Hungary15.  For  this  reason,  the 
ceutcrs-of-mass  of  the  test  bodies  used  were  separated  by  a  vertical  distance  of  about  21  cm, 
which  had  the  unfortunate  side  effect  of  making  Eotvos’  apparatus  significantly  more  sensitive 
to  vertical  gravity  gradients  than  to  anomalous  accelerations.  To  cancel  out  the  effects  of 
gravity  gradients,  EPF  combined  the  results  of  measurements  made  with  the  torsion  balance 
oriented  in  different  directions.  The  H2O-CU  comparison,  for  example,  was  carried  out  by 
measuring  the  difference  in  the  deflection  of  the  torsion  bar  when  it  was  oriented  North-South 
versus  South-North,  and  East- West  versus  West-East.  This  was  done  for  both  H2O  compared 
to  Pt  (where  the  Pt  was  loaded  at  the  upper  position  of  the  torsion  balance),  and  for  Cu 
compared  to  Pt.  The  effects  due  to  the  gravity  gradients  could  then  be  eliminated  by  taking 
the  appropriate  difference  of  the  deflections  measured  using  H20-Pt  versus  those  measured 
using  Cu-Pt. 

In  the  case  of  certain  comparisons  (e.g.,  magnalium-Pt  or  snakewood-Pt),  the  sample 
could  be  machined  to  the  desired  specifications,  and  hence  did  not  need  to  be  stored  in  a 
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Figure  1:  Model  of  container  +  sample  for  Eotvos,  Pckar,  and  Fcketc  experiment.15 


container.  Where  this  was  not  possible  ( e.g .,  H2O-C11  or  Asbcstos-Cu),  the  sample  was  placed 
in  a  brass  vial,  and  the  outside  dimensions  and  masses  for  the  various  containers  arc  given 
in  Ref.  15.  Since  the  densities  of  the  samples  could  estimated  from  the  composition  of  the 
respective  samples,  there  is  sufficient  information  to  allow  the  separate  determination  of  the 
holding  volume  of  the  brass  vials,  and  the  volume  of  each  sample.  In  many  cases,  the  volume 
of  the  sample  turned  out  to  be  smaller  that  the  volume  of  the  brass  vial,  and  in  this  case  the 
model  for  brass  vial+sample  shown  in  Fig.  1  was  used.  If  we  assume  that  the  thickness  of  the 
wall  of  the  vials,  t,  was  the  same  for  all  samples  contained  in  a  vial,  and  was  equal  to  that  of 
the  end  caps,  then  we  obtain  t  —  0.025  cm.  For  a  solid  cylinder  of  radius  R  and  length  L  we 
have 

<fl!>  =  (19) 

where  ( R '?)  is  measured  relative  to  the  center  of  mass.  For  the  composite  sample  given  by 
Fig.  1, 


(R2)  =  2^| PsL.Rl  +  l-L\  +  \h2  +Pc[(L2R*  -  Z-iRf) 

+  -  L\R\)  +  \h2(l2R*  -  LrR2)]  j, 


(20) 


where  mt  is  the  total  mass  of  the  brass  vial  +  sample,  La  =  m,/(irR2),  and  m„  is  the  mass  of 
the  sample.  Using  Eqs.  (19)  and  (20),  we  then  obtained  estimates  for  (R2)  for  each  sample, 
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Table  I:  Estimated  values  of  (iZ2)  for  the  various  samples  used  in  the  EPF  experi¬ 
ment.  Here  (R2)  is  given  in  units  of  cm2. 


Sample  I 

(RJ) 

Sample  II 

<K2> 

Magnalium 

12.0 

Pt 

3.0 

8.9 

Snakewood 

48.1 

Pt 

3.0 

45.1 

Cu 

3.5 

Pt 

3.0 

0.5 

Ag-Fe-S04 

17.3 

Ag-Fe-S04 

17.2 

0.0 

h2o 

15.5 

Cu 

3.6 

11.9 

•1 

15.8 

Cu 

3.5 

12.3 

CuS04  •  5H20 

46.4 

Cu 

5.5 

40.9 

II 

49.8 

Cu 

5.5 

44.3 

CuS04  Solution 

27.0 

Cu 

5.5 

21.4 

II 

26.7 

Cu 

5.5 

21.2 

Asbestos 

41.1 

Cu 

5.5 

35.6 

II 

40.5 

Cu 

5.5 

35.0 

Tallow 

19.4 

Cu 

5.5 

13.9 

II 

19.1 

Cu 

5.5 

13.6 

and  A(R2)  for  each  pair,  and  our  results  are  shown  in  Table  I. 

These  values  of  A(R2)  were  fitted  against  the  acceleration  differences  A/c  =  A a/h  ■  q 
mc;i8urcd  by  EPF  using  a  model  of  the  form 

A/c  =  aA(R2)  +  6,  (21) 

and  the  results  are  shown  in  Fig.  2  (solid  curve).  It  evident  from  Fig.  2  that  Eq.  (21)  is  a  very 
poor  model  of  the  EPF  data.  More  quantitatively,  we  find  x2  =  43  for  7  degrees  of  freedom 
for  this  fit.  Fitting  to  only  the  “Method  III”  data  points  (c./.,  Refs.  3  or  15),  a  somewhat 
better  fit  is  obtained  (dashed  line  in  Fig.  2).  However,  the  resulting  line  fails  to  pass  through 
the  origin,  and  hence  must  also  be  rejected  on  physical  grounds. 

SUMMARY 

The  principal  results  of  our  paper  are  presented  in  Eqs.  (7)-(12).  They  indicate  that, 
irrespective  of  the  functional  form  of  the  non-Newtonian  potential,  an  apparent  composition- 
dependent  effect  will  arise  in  an  Eotvos-type  experiment  by  virtue  of  the  fact  that  V2$(f )  /  0. 
Failure  to  detect  such  an  effect,  at  some  level,  is  thus  an  (almost)  model-independent  test  for 
the  presence  of  a  non-Newtonian  interaction.  In  practice  this  test  will  be  most  sensitive,  and 
hence  most  useful,  as  a  null  test  for  forces  which  vary  significantly  in  space  over  a  distance 
scale  comparable  to  the  size  of  the  test  masses.  We  note  in  passing  that  in  the  presence  of  a 
non-Newtonian  coupling,  the  finite-size  effect  produces  an  analog  of  the  Nordvedt  effe*  This 
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Figure  2:  Results  of  a  fit  of  A (R2)  from  Table  I  versus  Ak  obtained  by  Eotvos, 
Pekar,  and  Fekete.15  The  solid  curve  represents  a  fit  to  all  of  the  data  points,  and 
the  dashed  line  represents  a  fit  to  only  the  “Method  IIP  data  points.  Both  fits  fail 
to  give  a  satisfactory  result,  as  is  discussed  in  the  text. 

will  be  discussed  elsewhere,  as  will  be  the  details  of  possible  laboratory  experiments  utilizing 
the  finite-size  effect. 
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IS  THE  EOTVdS  EXPERIMENT  SENSITIVE  TO  SPIN? 

A.  M.  Hall  and  H.  Armbruater 
Department  of  Physics 
Virginia  Commonwealth  University 
Richmond,  VA  23284  USA 

E.  Fiachbach  and  C.  Talmadge 
Department  of  Physics 
Purdue  University 
West  Lafayette,  IN  47907  USA 


Recently,  Hall  and  Armbruster  have  introduced  a  phenomeno¬ 
logical  spin-dependent  charge,  which  they  have  suggested  may 
account  for  the  nonzero  acceleration  differences  reported  in  the 
Eotvos  experiment.  They  discuss  three  possible  explanations  for 
the  origin  of  this  charge.  This  paper  examines  specifically  one 
of  those  explanations,  namely,  that  this  charge  is  evidence  for  a 
“fifth  force”  in  Nature.  In  doing  so,  we  explore  the  implications 
of  fifth  force  models  based  on  this  charge,  and  study  the  question 
of  whether  this  charge  can  be  derived  from  a  more  general  theory. 


INTRODUCTION 


There  is  considerable  experimental  and  theoretical  interest  at  present  in 
the  possibility  of  deviations  from  the  predictions  of  Newtonian  gravity.  Part 
of  the  motivation  for  this  renewed  interest  was  a  reanalysia1,2  of  the  classic 
experiment  of  Eotvos,  Pek&r,  and  Fekete  (EPF)3  which  uncovered  in  the  EPF 
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data  »  correlation  suggesting  the  existence  of  ft  new  intermediate- range  coupling, 
Vs,  to  baryon  number  B.  The  potential  energy  V(r)  for  two  point  masses 
in  the  presence  of  the  “fifth  force”  potential  energy  Vs(r)  and  the  Newtonian 
gravitational  potential  energy  Vjy(r)  is  then  given  by3 

V(r)  =  VN(r)  +  V5(r)  =  -Goc^  +  /3^e-'/\  (1) 


where  m*  and  Bi  are  the  mass  and  baryon  number  of  *,  respectively,  G0 0  is 
the  Newtonian  gravitational  constant,  f  is  the  coupling  strength  of  the  new 
interaction  (the  analog  of  the  electric  charge  e  for  electromagnetism),  and  A  is 
the  range  of  the  proposed  new  force.  If  the  masses  are  expressed  in  units  of 
m^H1)  —  1.00782519(8)u,  so  that  m,  =  then  Eq.  (1)  can  be  written 

in  the  form3  _ . 

V(r)  =  -Goo^  (l  +  aye-''*)  , 


BiBic 
an  = - 


(  = 


f2 


(2) 


GaoTTl 


T 


Since  otij  depends  on  the  compositions  of  i  and  j  through  the  ratios  BJfii 
and  Bj/fij,  the  coupling  in  (2)  leads  to  an  apparent  violation  of  the  Weak 
Equivalence  Principle  (WEP),  which  postulates  the  equivalence  of  gravitational 
and  acceleration  effects.  The  acceleration  difference  A a,-3«  of  two  objects  j  and 
j'  towards  t  would  then  be  directly  proportional  to 


(3) 


and  it  was  this  specific  correlation  between  A  (B/p)jy  and  the  EPF  data  for 
Aajji  which  stimulated  discussion  of  a  “fifth  force.” 


The  model  defined  by  Eqa.  (l)-(3)  has  been  discussed  extensively  in  the 
literature1,3,4-9.  It  has  been  noted  in  Refs.  1  and  2  that  B/n  varies  across  the  Pe¬ 
riodic  Thble  in  a  way  that  is  quite  different  from  that  of  other  proposed  charges, 
such  as  lepton  number  L.  That  B/n  is  not  (even  approximately)  a  monotonic 
function  of  atomic  number  Z,  may  help  to  explain  why  attempts  to  understand 
the  correlation  in  the  Eotvos  data  in  terms  of  conventional  physics2’10'11  thus 
far  have  been  unsuccessful.  Experiments  aimed  at  reproducing  the  EPF  results 
have  also  been  unsuccessful,  however,  and  by  now  there  can  be  little  doubt  that 
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the  model  of  V5  given  in  Eqs.  (1 )— (3)  is  incompatible  with  existing  data.  Var¬ 
ious  attempts  to  generalize  this  model  have  thus  far  failed  to  explain  bow  the 
correlation  in  the  EPF  data  could  be  compatible  with  the  results  of  modern 
experiments.  The  implications  of  these  experiments  is  that  the  Eotvos  experi¬ 
ment  is  flawed,  yet  no  convincing  model  has  been  put  forward  to  date  which  can 
account  for  the  original  EPF  data  in  conventional  terms.  We  refer  the  interested 
reader  to  the  discussions  in  Reis.  2, 12,  and  13  for  further  consideration  of  these 
points. 

Before  turning  to  the  question  of  the  possible  relevance  of  spin  in  the  EPF 
experiment,  it  is  worth  emphasizing  that  both  the  modern  repetitions  of  the 
EPF  experiment  and  the  related  Galileo  (free-fall)  experiment14’15  have  been 
designed  to  optimize  their  sensitivity  to  the  coupling  in  Eqs.  (l)-(3),  or  to 
certain  alternatives  to  it.  These  experiments  may  not  be  equally  sensitive  to  a 
spin-dependent  charge,  and  thus  to  a  spin-dependent  “fifth  force”  V55. 

PHENOMENOLOGY  OF  A  SPIN-DEPENDENT  CHARGE 


Recently  Hall  and  Armbruster  (HA)  have  noted16  that  a  correlation  similar 
to  that  found  in  Ref.  1  arises  in  the  Eotvos  data,  if  the  baryon  number  B  of  a 
nucleus  is  replaced  by  a  charge  Q  defined  by 


Q  =  M6\ 


for  J  >  0, 
for  7  =  0, 


(4) 


where  M  is  the  mass  of  the  nucleus,  and  J  is  its  nuclear  spin.  (Note  that  we 
use  the  notation  J  rather  than  the  more  conventional  I  to  avoid  any  possible 
confusion  with  the  nuclear  isospin  which  was  another  suggestion  for  Q.)  For 
elements  with  more  than  one  isotope,  the  total  Q  is  obtained  in  the  usual  way3 
by  weighting  contributions  from  each  isotope  according  to  its  isotopic  abundance 
rk: 

Q  =  '£/Mk6krk.  (5) 

k 

Before  returning  to  the  question  of  the  physical  significance  of  Q,  we  present 
in  Tkble  I  the  values  of  Q/fi  for  the  natural  elements,  and  in  the  Table  II  the 
values  A (Q/p)  for  the  EPF  samples,  where  for  samples  i  and  j,  A(Q/fi)  = 
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Table  I:  Average  value  of  Q/n  using  Eq.  (4)  for  the  first  92  elements  of  the 
Periodic  Table. 


Element 

Q/p 

Element 

Q/p 

Element 

Q/p 

Hydrogen 

Germanium 

0.44976 

Europium 

1.00000 

Helium 

Arsenic 

mm 

Gadolinium 

0.30158 

Lithium 

Selenium 

0.07381 

Terbium 

1.00000 

Beryllium 

Bromine 

mm 

Dysprosium 

0.43735 

Boron 

Krypton 

0.11427 

Holmium 

1.00000 

Carbon 

: 

Rubidium 

mm 

Erbium 

0.22896 

Nitrogen 

m 

Strontium 

0.06963 

Thulium 

mm 

Oxygen 

Yttrium 

1.00000 

Ytterbium 

0.30259 

Fluorine 

1.00000 

Zirconium 

0.11191 

Lutetium 

mm 

Neon 

0.00267 

Niobium 

mm 

Hafnium 

0.32121 

Sodium 

mm 

Molybdenum 

0.25119 

Tantalum 

mm 

Magnesium 

0.10412 

Technetium* 

mm 

Tungsten 

0.14331 

Aluminum 

Ruthenium 

0.29497 

Rhenium 

mm 

Silicon 

0.04849 

Rhodium 

mm 

Osmium 

0.17595 

Phosphorus 

1.00000 

Palladium 

0.21909 

Iridium 

mm 

Sulfur 

0.00782 

Silver 

mm 

Platinum 

0.33779 

Chlorine 

Cadmium 

0.24890 

Gold 

1.00000 

Argon 

0.00000 

Indium 

Mercury 

0.29946 

Potassium 

1.00000 

Tin 

0.16424 

Thallium 

1.00000 

Calcium 

0.00155 

Antimony 

mm 

Lead 

0.22578 

Scandium 

1.00000 

Tellurium 

0.07679 

Bismuth 

Titanium 

0.12772 

Iodine 

mm 

Polonium* 

Vanadium 

1.00000 

Xenon 

0.47074 

Astatine* 

Chromium 

0.09683 

Cesium 

1.00000 

Radon* 

WliffliS 

Manganese 

1.00000 

Barium 

0.17759 

Francium* 

0.50106 

Iron 

0.02233 

Lanthanum 

1.00000 

Radium* 

0.00000 

Cobalt 

1.00000 

Cerium 

0.00000 

Actinium* 

1.00000 

Nickel 

0.01235 

Praseodymium 

mm 

Thorium* 

liiiMi'i'il 

Copper 

mm 

Neodymium 

0.20396 

Protactinium* 

EE23 

Zinc 

Gallium 

0.04207 

1.00000 

Promethium* 

Samarium 

1.00000 

0.28324 

Uranium 

0.00711 

*No  stable  isotopes 


{Q/n)i  -  These  are  plotted  against  the  acceleration  differences  (A/c  in 

the  EPF  notation)  in  Fig.  1.  Fitting  these  data  to  a  straight  line, 

Ak  =  7  MQ/p)  +  (6) 
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Table  II:  Calculated  values  of  A (Q/p)  using  Eq.  (4)  for  the  EPF  samples  versus 
the  quoted  EPF  values  of  Ak.  The  misprint  in  the  sign  quoted  by  EPF  for  the 
RaBrj-Pt  datum  has  been  corrected,  as  discussed  in  Ref.  2. 
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Samples 

Legend 

MQ/ri 

10®  Ak 

Asbestos-Cu 

a 

-0.568 

-2±1 

HjO-Cu 

b 

-0.611 

— 5±1 

Tallow-Cu 

c 

-0.493 

-3±1 

CuS04(sol’n)-Cu 

d 

-0.505 

— 4±1 

CuS04  •  5HjO-Cu 

e 

•0.441 

-3±1 

Snakewood-Pt 

f 

-0.274 

— 1±2 

Ag-Fe-S04 

S 

+0.000 

+0±1 

RaBrj-Pt 

h 

+0.343 

+1±2 

Magnalium-Pt 

i 

+0.578 

+4±1 

Cu-Pt 

j 

+0.668 

+4±2 

Fig.  1  Plot  of  A (Q/n)  obtained  from  Eq.  (4)  versus  Ak, 
using  the  data  from  Table  II.  The  straight  line  is  the  result 
of  a  least  squares  fit  to  these  data,  as  discussed  in  the  text. 
The  labels  a-j  are  defined  in  Table  II. 


we  find 

7  =  (6.5  ±  0.8)  x  10~9,  6  =  (-0.1  ±  4.1)  x  10~10, 

XJ  =  5.1  (8  degrees  of  freedom). 


The  quality  of  this  fit  is  very  similar  to  that  obtained  originally  in  Ref.  1  for  a 
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coupling  to  B:  Using  the  results  of  Table  V  in  Ref.  2,  the  corresponding  values 
of  7  and  6  are: 

7  =  (5.08  ±  0.61)  x  10"*,  6  =  (0.87  ±  0.59)  x  10'®, 

(8) 

X2  =  5.2  (8  degrees  of  freedom). 

HA  note  in  Ref.  16  that  the  effect  of  J  =  1/2  nuclei  may  be  removed  (by 
setting  6  =  0  for  J  =  1/2  nuclei  in  Eq.  (4))  without  significantly  affecting  the 
correlation  in  Fig.  1.  Further,  the  EPF  data  do  not  discriminate  between  other 
variants  of  the  charge  Q  (such  as  replacing  6  by  [J(J  +  1)]1/J),  as  these  variants 
also  evidence  a  correlation  with  the  EPF  data.1*  Even  the  extreme  assumption 
that  Q  =  0  for  all  elements  but  Cu  and  A1  leads  to  a  correlation.16  While  these 
variants  may  or  may  not  be  significant,  we  consider  here  only  the  charge  Q  given 
in  Eq.  (4). 


MODELS  OF  THE  SPIN-DEPENDENT  CHARGE  Q 

In  this  section  we  discuss  preliminary  models  aimed  at  deriving  a  spin- 
dependent  interaction  from  a  more  fundamental  theory.  This  is  desirable  partly 
because  one  must  also  know  the  spatial-dependence  of  whatever  coupling  leads 
to  Eq.  (4),  in  order  to  design  appropriate  experiments.  This  is  especially  im¬ 
portant  for  composition-independent  experiments,  which  directly  measure  the 
spatial  variation  of  lfo(r),  as  we  discuss  below. 

One  of  the  challenges  in  constructing  a  model  of  the  EPF  data  based  on  the 
spin-dependent  charge  Q,  is  that  the  EPF  samples  were  presumably  unpolarized. 
It  is  reasonable  to  assume  that  (J)  a  0  for  the  samples  as  well  as  for  the  source, 
although  bulk  matter  may,  however,  have  a  small  net  polarization  due  to  the 
Earth’s  magnetic  field  At  room  temperature  (T  =  300° K),  however,  the 
polarization  V  is  expected  to  be  of  order17 

V  *  a  1  x  l<r10B®  (Gauss),  (9) 

kbT 

where  ns  i*  the  nuclear  magneton,  and  kB  is  the  Boltzmann  constant.  For 
B&  9 £  1  Gauss  this  gives  V  a  1  x  10~10.  Since  any  effect  would  depend  on  the 
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product  P,Pd  for  the  source  sod  detector,  experiments  sensitive  to  PtPd  would 
be  suppressed  by  a  factor  of  order  10-20  relative  to  those  with  P,  w  P4  a  1. 
This  means  that  the  effective  strength  of  the  “gravitational”  interaction  between 
spin-polarized  samples  (for  which  Pal)  should  be  «  1020  stronger  than  the 
strength  of  the  claimed  signal  in  the  EPF  experiment,  which  we  may  take  to  be 
of  order  (10~2Goo  -  10_sGoo)-  This  conclusion  is  significantly  at  variance  with 
the  results  of  recent  experiments  on  spin-polarized  samples18  which  set  upper 
limits  on  the  strength  of  a  spin-dependent  gravitational  coupling  to  electron  spin 
at  the  level  of  a  10~3G<X-  These  experiments  can  also  be  used  to  set  limits  on  a 
coupling  to  nuclear  spin,  because  the  polarized  electrons  produce  a  net  hyperfine 
field  at  the  nucleus,  which  in  turn  polarizes  the  nucleus.  Since  this  field  is  much 
larger  than  B9,  it  would  be  difficult  to  understand  how  an  effect  could  have 
shown  up  (incidentally)  in  the  EPF  experiment,  but  not  in  those  of  Ref.  18. 


The  implication  of  the  preceding  discussion  is  that  if  spin  is  indeed  relevant 
in  the  EPF  experiment,  it  must  be  so  in  a  way  that  does  not  depend  on  the  net 
polarization  of  the  nuclei  in  the  test  masses.  One  way  this  can  happen  is  if  the 
spin  were  to  contribute  a  term  to  the  nuclear  mass-energy,  and  that  term  coupled 
anomalously  to  gravity.  If  M/  and  Mg  denote  the  inertial  and  gravitational 
masses  of  a  nucleus,  then  the  content  of  the  WEP  is  that  Mi  =  Mg-  Suppose, 
however,  that  a  particular  contribution  ea  of  the  nuclear  mass-energy  coupled 
anomalously  to  gravity  and  violated  the  WEP  so  that  Mi  ^  Mq.  We  can  then 
define  a  parameter  r)a, 

Af>  "1+ ^  fe)  +  (10) 


which  measures  the  strength  of  the  anomalous  gravitational  coupling  of  ea.  The 
parametrization  in  (10)  has  been  widely  used  in  the  literature  to  set  limits  on 
the  coupling  of  various  possible  energy  terms  in  the  nuclei's  using  the  Eotvos 
experiment.10  Consider,  for  example,  the  acceleration  of  nucleus  1  towards  the 
Earth  in  the  presence  of  (10): 


M/iai  = 


°i  = 


GoqMqMgx  _  GooMqMh  /,  .  * 

~  d2  —  d2  l*T*al/ 

- +  *al)  =  -0©(1  +  *al)- 


(ID 
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It  follows  that 


"V’ ~ K"-i° [(«/),  (m,)J-  <i2) 

Since  ( ea/Mi )  is  in  general  different  for  dissimilar  nuclei,  it  follows  that  Eq.  (12) 
can  be  used  to  set  constraints  on  qa.  In  terms  of  this  picture,  the  problem  of 
justifying  the  form  of  Q  in  Eq.  (4)  reduces  to  finding  an  energy  term  in  the 
nucleus  which  depends  on  the  product  MS,  or  some  expression  similar  to  it. 


To  illustrate  how  a  dependence  of  some  energy  ea  on  MS  could  come  about, 
we  consider  for  illustrative  purposes  the  gravitational  contribution  to  the  spin- 
orbit  energy  £so  of  a  single  nucleon  outside  a  filled  or  partially-filled  shell.  This 
is  given  by20 


£S0  =  r 


3  Goo  M  ~  - 


2  mtPr3 


LS, 


(13) 


where  M  is  the  mass  of  the  shell,  m  is  the  nucleon  mass,  L  and  §  are  the  nucleon 
orbital  and  spin  angular  momenta,  and  r  is  the  distance  from  the  nucleon  to  the 
center  of  the  shell.  The  contribution  from  eso  appears  to  have  the  property  of 
being  proportional  to  ML  •  S  -  ■  |M[J2  -  I?  —  S2],  and  hence  vanishes  when 
there  are  no  unpaired  nucleons,  as  does  the  charge  MS  ia  (4).  Since  nuclei  have 
an  approximately  constant  density  p,  we  can  set  M  3  *»r3p,  however,  so  that 


T.'n pGooL  •  S 

£so  — - - 

mr 


(14) 


This  contribution  to  the  nuclear  energy  does  not  grow  with  M  as  desired,  and 
thus  is  not  a  candidate  for  the  charge  of  Eq.  (4).  A  priori,  there  is  an  even  more 
serious  objection  to  this  toy  model,  quite  apart  from  the  obvious  problem  of  the 
magnitude  of  eso-  This  has  to  do  with  the  fact  that  the  same  (gravitational) 
interaction  which  couples  the  valence  nucleon  to  the  filled-shell  (core),  also  cou¬ 
ples  the  core  nucleons  to  one  another.  Thus  if  the  spin-orbit  gravitational  energy 
coupled  in  some  anomalous  way  to  gravity,  the  spin-independent  energy  would 
as  well.  This  might  well  lead  to  a  violation  of  the  WEP,  but  it  would  be  predom¬ 
inantly  a  spin-independent  effect,  which  is  not  what  we  are  seeking.  Evidently, 
similar  arguments  apply  mutatis  mutandis  for  other  contributions  to  the  total 
spin-orbit  interaction.21  Since  any  contribution  to  the  energy  from  the  exchange 
of  bosons  in  the  normal  spin-parity  series  ( Jp  =  0+,  l-,  2+ , . . .)  will  lead  to  the 
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same  problem,  the  best  hope  for  constructing  a  model  of  the  spin-dependent 
charge  along  these  lines  appears  to  arise  from  the  exchange  of  bosons  in  the 
abnormal  parity  series,  Jp  =  0“,  1+,  2~, . . . 

Consider,  for  example,  the  potential  arising  from  the  exchange  of  a  massive 
axial-vector  ( Jp  =  1+)  field,12 

e-r/*A 

VA(r)  =  9\d1.d2——,  (15) 

where  3\  and  &2  are  the  spins  of  the  interacting  nucleons,  \A  =  h/mAc  is  the 
Compton  wavelength  of  the  1+  quantum,  and  gA  ia  &n  appropriate  coupling 
constant.  The  interaction  in  Eq.  (15)  describes  not  only  the  coupling  of  the 
valence  nucleon  to  a  nucleon  in  the  core,  but  also  the  interaction  of  the  core 
nucleons  among  themselves.  In  the  latter  case  £  -fc  may  average  to  zero 

for  the  core  nucleons,  depending  on  the  details  of  the  interaction  and  on  the  form 
of  the  core  wavefunction.  Due  to  exchange  (Pauli)  effects,  this  is  not  necessarily 
the  case  for  3\  ■  £  <?*,  where  1  denotes  a  valence  nucleon  and  k  a  nucleon  in 
the  core.  What  this  means  is  that  to  0(g\)  there  may  be  no  contribution  from 
VA  to  the  energy  of  a  nucleus,  except  for  possible  contributions  from  valence 
nucleons  interacting  with  the  core.  If  we  then  suppose  that  this  term  has  an 
anomalous  coupling  to  gravity,  or  perhaps  to  another  long-range  gravity-like 
field,  then  we  have  at  least  the  outline  of  a  model  which  could  give  rise  to  some 
of  the  qualitative  features  of  the  spin-dependent  charge  in  Eq.  (4). 

To  summarize,  the  proposed  mechanism  for  modeling  the  spin-dependent 
charge  as  a  fifth  force  is  based  on  an  anomalous  coupling  of  the  spin-dependent 
energy  arising  from  VA  in  Eq.  (15).  This  energy  term  could  have  the  property 
of  vanishing  for  J  =  0  nuclei,  while  at  the  same  time  being  nonzero  for  nuclei 
with  a  valence  nucleon  outside  a  core.  In  this  case,  the  total  energy  should 
scale  with  the  size  of  the  core,  and  these  features  may  qualitatively  simulate  the 
behavior  of  Q  =  MS  (or  one  of  the  alternatives  to  Q  which  also  explains  the 
EPF  data).  We  emphasize  that  considerable  effort  will  be  required  before  we 
will  know  whether  this  or  any  other  specific  model  actually  works.  Nonetheless, 
the  preceding  arguments  indicate  that  a  spin-dependent  charge,  such  as  that 
proposed  in  Ref.  16,  may  in  fact  have  a  field- theoretic  justification. 
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IMPLICATIONS  OF  A  SPIN-DEPENDENT  CHARGE 


We  outline  in  this  section  some  of  the  phenomenological  implications  of  the 
spin-dependent  charge  in  Eq.  (4).  Eq.  (3)  can  be  generalized  by  noting  that 
whatever  the  specific  form  of  the  charge  Q  that  determines  the  strength  of  Vjs, 
the  composition-dependent  acceleration  difference  will  always  be  proportional 
to  the  product 

<(?)  (is) 

V/*/ source  V^/ detector 

Here  “detector”  refers  to  the  two  masses  whose  accelerations  are  being  compared 
as  they  accelerate  toward  a  third  object  which  is  defined  as  the  “source.”  The 
product  S  =  9,A<7a,  termed  the  “sensitivity  function”  by  Adelberger,  et  al.23  is 
one  of  the  factors  that  determines  the  sensitivity  of  a  given  experiment  to  the 
putative  fifth  force.  Another  factor  is  the  source  integral13  f,  which  depends 
in  turn  on  the  the  r-dependence  of  V55.  Since  we  do  not  as  yet  have  even  a 
rudimentary  phenomenological  theory  which  predicts  the  r-dependence  of  V55, 
we  will  limit  ourselves  to  what  can  be  learned  from  Q  alone. 

Experiments  with  Large  Geophysical  Sources 

This  category  includes  the  experiments  of  Adelberger,  et  aL,23  Bizzeti,  et 
al.,24  Boynton,  et  al.,25  Fitch,  et  al.,1*  Kuroda  and  Mio,15  Niebauer,  et  al.14  and 
Thieberger.27  These  experiments  have  in  common  the  feature  that  their  source 
is  predominantly  SiOj  for  which  q„  3  0.  This  means  that  for  the  class  of  models 
we  are  considering,  the  relative  sensitivity  of  the  different  experiments  would 
depend  crucially  on  the  presence  at  different  sites  of  “impurities”  with  q,  ^  0. 
The  same  remark  would  apply  to  experiments  carried  out  at  sites  rich  in  CaCC>3 
(e.g.,  limestone),  which  is  another  common  mineral.  Further,  for  variants  of  the 
model  in  which  the  spin  charge  Q  =  0  for  hydrogen  (as  in  the  model  discussed  in 
the  previous  section),  the  presence  or  absence  of  water  in  the  source  is  irrelevant. 
In  this  context  it  is  interesting  to  note  that  if  Q  =  0  for  hydrogen  then  the 
constraints  imposed  on  V$s  by  the  solar  Eotvos  experiments25-30  are  reduced 
by  a  factor  of  ~  103  from  their  quoted  values. 
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Experiments  with  Laboratory  Sources 

A  number  of  experiments  with  laboratory  sources31  have  been  carried  out 
which  may  set  interesting  limits  on  Vjs,  once  the  distance  dependence  is  known. 
Laboratory  experiments  with  Pb  sources,  however,  would  have  very  limited  sen¬ 
sitivity  to  the  spin  charge  discussed  here,  as  qpb  —  0  using  the  charge  of  Gq.  (4). 
It  is  interesting  to  note  that  the  implications  of  the  spin  charge  model  are  similar 
to  those  of  the  isospin  model  (related  to  V5)  of  Boynton,  et  a \LM  and  Adelberger, 
et  a/.,13  with  one  critical  difference.  In  the  latter  model  Q  —  N  -  Z  =  where 
N  and  Z  denote  the  numbers  of  neutrons  and  protons  in  the  sample.  Since 
is  relatively  large  t or  Pb,  experiments  with  laboratory  sources  of  Pb  can 
set  stringent  limits  on  the  strength  of  a  coupling  to  It.  For  a  summary  of  recent 
laboratory  results,  see  Fischbach  and  Talmadge.12  Unlike  the  case  of  isaspin,  the 
spin-dependent  charge  of  Eq.  (4)  is  non-negative.  This  means  that  there  can  be 
no  cancellations  as  occur  for  /,  between  the  contributions  from  water  (for  which 
1,1  p  a  —0.112)  and  the  slightly  positive  Iz/p  of  mineral  impurities. 

Pumped  Lake  Eotvos  Experiments 

Bennett33  has  carried  out  an  Eotvos  experiment  using  a  pumped  water  fa¬ 
cility  as  his  source.  This  experiment  can  set  interesting  limits  on  a  coupling 
to  /,  or  to  Q  in  (4).  The  EPF  data,  however,  are  also  consistent  with  the 
assignment18  6  =  0  for  J  <  1  (and  the  toy  model  above  does  not  contain  a 
coupling  to  hydrogen).  In  this  case  the  Bennett  experiment  would  be  insensitive 
to  Q,  which  is  then  compatible  with  his  null  result. 

Tower,  Mine/Borehole  and  Lake  Experiments 

In  the  notation  of  Eq.  (9)  these  composition-independent  experiments4,33-33 
are  sensitive  to  the  product  qtqd  where  q 4  is  the  charge  of  the  standard  mass  in 
the  gravimeter.  Typically  measurements  are  carried  out  with  similar  gravime¬ 
ters,  so  any  substantial  differences  among  experiments  would  arise  from  the 
source  charge  qt.  At  present  there  is  no  compelling  evidence  for  deviations  from 
Newtonian  gravity  in  any  of  these  experiments,  but  since  lingering  anomalies 
remain,  measurements  are  continuing.  It  would  be  difficult  at  this  point  to  use 
the  spin-charge  model  to  analyze  these  experiments,  since  at  the  present  stage 
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the  model  makes  no  statement  about  the  spatial  variation  of  Vjs(r)  to  which 
these  experiments  are  sensitive. 


SUMMARY  AND  CONCLUSIONS 


The  phenomenological  spin-dependent  charge  in  Eq.  (4)  introduced  by  Hall 
and  Armbruster16  provides  an  alternative  to  the  original  hypercharge  model  in 
Refs.  1  and  2  as  an  explanation  of  the  EPF  results.  In  this  paper,  we  have  ex¬ 
plored  the  specific  possibility  of  modeling  this  charge  as  that  of  a  spin-dependent 
“fifth  force”  V55.  Although  the  spin-dependent  charge  of  Eq.  (4)  as  originally 
proposed  does  not  appear  to  arise  in  this  context  from  a  more  fundamental 
theory,  some  variants  of  a  model  based  on  the  charge  of  Eq.  (4)  may  have  a 
field-theoretic  basis,  as  we  have  discussed.  The  success  of  a  spin-dependent 
charge  in  accounting  for  the  EPF  data  raises  the  question  of  whether  there  are 
variables  other  than  B  or  Q  in  (4)  which  could  explain  the  EPF  data.  If  so,  then 
we  might  be  tempted  to  ask  whether  these  data  have  any  deep  physical  signifi¬ 
cance,  if  they  can  be  accounted  for  by  such  seemingly  different  variables.  On  the 
other  hand,  among  the  many  physical  variables  which  have  been  tried,3,36  B  and 
Q  are  the  only  two  discovered  to  date  which  give  the  correct  correlation.  Since 
both  involve  non-classical  variables,  and  may  be  related  to  each  other  at  a  deeper 
level,  it  is  possible  that  they  may  be  the  signal  for  new  physics.  Although  it  is 
still  too  early  to  assess  the  physical  significance  of  the  spin-dependent  charge,  or 
the  success  in  modeling  it  as  a  charge  of  a  spin-dependent  fifth  force,  it  is  clear 
that  the  correlation  noted  in  Ref.  16  raises  new  and  interesting  possibilities  that 
are  worth  pursuing. 
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We  show  that  in  certain  classes  of  theories  the  spatial  variation  of  the  non-Newtonian  potential 
would  be  dominantly  an  exponential  rather  than  a  Yukawa  potential,  and  we  compare  the  phenom¬ 
enological  interpretation  of  the  existing  data  in  the  exponential  and  Yukawa  models.  We  also  show 
that  generalized  forms  of  the  exponential  potential  can  arise  naturally  from  simple  mass  spectra. 
Although  such  models  cannot  reconcile  all  of  the  existing  data  on  non-Newtonian  gravity,  they  have 
novel  properties  that  can  be  directly  studied  experimentally. 


I.  INTRODUCTION 

A  number  of  experiments  are  presently  underway 
to  search  for  both  composition-independent  and 
composition-dependent  deviations  from  Newtonian  gravi¬ 
ty.1-8  To  date  there  is  no  compelling  evidence  for  the 
existence  of  such  deviations,  although  some  unexplained 
experimental  anomalies  remain  to  be  fully  understood. 5,8 
The  accumulation  of  a  large  number  of  null  results  can  be 
used  to  infer  stringent  constraints  on  various  theories 
which  would  attribute  any  departures  from  Newtonian 
gravity  to  the  existence  of  new  forces.  The  non- 
Newtonian  effects  expected  in  such  theories  are  conven¬ 
tionally  described  in  terms  of  a  modified  expression  for 
the  potential  energy  V[r)  of  two  point  masses  mu 
separated  by  a  distance  r. 

K(r)=— G— -1— ( 1  +ae~r/K)=  V„(r)+  V'(r)  . 


Here  Gx  is  the  Newtonian  constant  of  gravity,  and  the 
parameters  k  and  a,  respectively,  give  the  range  of  the 
new  force  and  its  strength  relative  to  gravity.  Also,  K'(  r) 
describes  the  correction  to  the  effective  gravitational  po¬ 
tential  arising  from  the  particular  non-Newtonian  in¬ 
teraction  we  are  considering  (which  in  this  case  is  a  Yu¬ 
kawa).  The  functional  form  of  the  Yukawa  contribution 
V'(r)=  VY(r)  in  Eq.  (1.1)  is  suggested  by  models  in  which 
Vyir)  arises  from  the  exchange  of  a  single  new  quantum 
with  mass  my  =  A.-1.  If  the  source  of  this  quantum  is  a 
charge  Q  =B  cos0'+/zsin0',  where  B  =N  +Z  is  baryon 
number,  and  IZ  =  N  ~Z  is  isospin,  then 

«*-$<Ci/iMe24*2>»  Mi.2="»i.2 /"»//•  <l-2> 

Here  £=f2/G xmj,  gives  the  coupling  strength  in  terms 
of  the  unit  of  charge  /,  and  mH=m(lHl).  For  later  pur¬ 
poses  it  is  instructive  to  exhibit  the  expression  for  the 
force  F(r)  implied  by  (1.1): 


F  <r)=-VF(r)= 


■Gxm]m2T 


[l+a{l+r/k)e 


-G{r)mtm2T 


G(r)  =  Gy(r)=G<r[l+a(l+r/A.)e"rA]  .  (1.3b) 

Each  experimental  limit  defines  a  contour  in  the  a-k 
plane  which  specifies  the  region  in  the  plane  excluded  by 
that  experiment.  A  summary  of  recent  limits  can  be 
found  in  Refs.  9-11. 

To  date  the  Yukawa  model  in  Eqs.  (1.1)— (1.3)  has  been 
the  most  widely  studied  framework  for  introducing  non- 
Newtonian  effects.  However,  there  are  no  compelling 
reasons  to  believe  that  the  phenomenological  non- 
Newtonian  coupling  must  have  the  simple  form  suggested 
by  Eq.  (1.3),  and  since  the  theoretical  implications  of  the 
experimental  data  are  quite  different  in  other  models,  a 
number  of  alternatives  to  Eq.  (1.3)  have  also  been  exam¬ 
ined.  For  example,  Moffat 12  has  studied  the  experimental 
consequences  of  his  nonsymmetric  gravity  theory,  in 
which  the  non-Newtonian  contribution  to  the  force  varies 
as  r~i,  rather  than  that  expected  from  a  Yukawa  poten¬ 
tial.  In  addition,  various  authors  have  considered  some 
of  the  implications  of  a  model  with  two  (nearly)  canceling 
Yukawa  potentials. 1,-15  The  latter  model  is  the  starting 
point  of  the  present  paper,  whose  objective  is  to  demon¬ 
strate  that  there  is  an  interesting  (and  hitherto  unex¬ 
plored)  limiting  case  of  two  canceling  Yukawa  potentials, 
in  which  the  resulting  potential  can  be  represented  as  an 
approximate  exponential. 15  ,6  By  analyzing  this  regime 
in  terms  of  a  single  exponential,  we  demonstrate  that 
such  a  coupling  leads  to  novel  phenomenology.  More¬ 
over,  by  representing  this  interaction  directly  as  an  ex¬ 
ponential  we  avoid  the  inevitable  computational  errors 
which  arise  when  the  contributions  from  two  (nearly) 
canceling  Yukawa  potentials  are  evaluated  numerically, 
as  is  discussed  below.  One  of  the  purposes  of  this 
analysis  is  to  compare  in  detail  the  phenomenological  im¬ 
plications  of  the  existing  data  in  the  Yukawa  and  ex¬ 
ponential  models,  in  order  to  establish  the  extent  to 
which  the  relative  sensitivities  of  different  experiments 
change  in  the  two  models. 

In  Sec.  II  below  we  obtain  the  exponential  potential  as 
the  limiting  case  of  two  canceling  Yukawa  potentials,  and 
discuss  models  in  which  the  necessary  cancellations  can 
take  place.  The  phenomenological  implications  of  the  ex¬ 
ponential  model  are  discussed  in  Sec.  Ill,  and  our  con¬ 
clusions  are  summarized  in  Sec.  IV.  Some  mathematical 
details  of  the  exponential  model  are  given  in  the  Appen¬ 
dixes. 
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II.  ORIGIN  OF  THE  EXPONENTIAL  POTENTIAL 


As  we  have  noted  in  Sec.  I,  the  exponential  potential 
arises  as  the  limiting  case  of  two  (nearly)  canceling  Yu¬ 
kawa  potentials.  In  this  section,  we  demonstrate  this  ex¬ 
plicitly,  and  discuss  at  the  same  time  models  in  which  the 
necessary  cancellations  could  come  about. 

Consider  the  interaction  of  two  nucleons  arising  from 
the  exchange  of  fields  <f>a  and  4>b  which  interfere  destruc¬ 
tively,  so  that  V‘ir)=  Va(r)—  Yb(r).  Denote  the  masses 
of  these  fields  by  mab  and  their  couplings  to  Q  by  fa  b.  If 
the  masses  and  coupling  constants  are  related  in  such  a 
way  that17 

ma=mfc[l+0(e)],  /„=/*[  l+0(5e>]  ,  (2.1) 

where  e,8  « 1,  then  it  is  straightforward  to  show  that  the 
leading  contribution  to  the  potential  V‘(r)=  Va(r)—  Vb(r) 
depends  exponentially  on  r  (see  Appendix  A): 


V'^V^r^f'Q^  ~ 

F  £ir)=-VV£(r)=f2Q,Q2 

-Ell rl-t  ^ A.  Q\ 

m\  ~  A 


-r/K 

c 

(2.2) 

A 

AA  e~"x ^ 

A  1  A2  '  ’ 

(2.3) 

Qk 

—  ?£(r.  A). 

(2.4) 

B-2 

Here  AA=Aa  — A6  =  l/ma  — 1  /mb,  f  =fb,  and  3E  is  the 
field  strength  for  a  point  source  m2  in  the  exponential 
model.  It  follows  from  the  preceding  discussion  that  we 
can  view  VEir)  as  a  simple  parametrization  of  the  limit¬ 
ing  case  of  two  nearly  canceling  Yukawa  potentials, 
which  has  novel  phenomenological  implications,  as  we 
discuss  below. 

We  turn  next  to  discuss  some  representative  models  of 
non-Newtonian  forces  which  illustrate  how  an  exponen¬ 
tial  potential  could  arise  in  specific  theories.  As  we  have 
already  noted,  the  possibility  of  two  nearly  canceling  Yu¬ 
kawa  potentials  has  been  investigated  by  a  number  of  au¬ 
thors,  including  Goldman,  Hughes,  and  Nieto13  (GHN). 
These  authors  and  others18  note  that  in  simple  phenome¬ 
nological  theories  the  exchange  of  a  vector  (Jp=\~) 
leads  to  a  repulsive  force,  whereas  scalar  Up= 0+)  ex¬ 
change  gives  an  attractive  force.  If  the  vector  and  scalar 
fields  were  in  turn  related  by  some  higher  symmetry,  such 
as  supersymmetry,  then  the  masses  and  couplings  of  these 
fields  might  be  sufficiently  close  for  a  substantial  cancella¬ 
tion  to  take  place.  However,  in  order  to  arrive  at  an  ex¬ 
ponential  potential  the  masses  and  coupling  constants 
must  satisfy  Eq.  (2.1)  and  whether  this  can  happen  in 
such  theories  is  not  known. 

A  less  obvious,  but  potentially  more  interesting,  mech¬ 
anism  for  producing  an  exponential  potential  is  to  consid¬ 
er  the  cancellation  between  two  scalar  fields,  rather  than 
between  a  scalar  and  a  vector  field  as  above.  Ordinarily 
the  exchange  of  a  scalar  field  coupling  to  a  simple  charge 
such  as  baryon  number  or  lepton  number  gives  rise  to  an 
attractive  force  as  we  have  noted  previously.  However, 
several  models  of  weak  gravitylike  forces  have  been  pro¬ 
posed  recently  in  which  a  scalar  field  couples  to  a  “com¬ 


posite”  charge,  i.e.,  one  which  is  a  linear  combination  of 
more  elementary  charges.  For  example,  in  the  theory  of 
Peccei,  Sola,  and  Wetterich19  (PSW),  scalar  exchange 
produces  a  coupling  between  nuclei  1  and  2  of  the  form 

-r/K 

r(r)=-G0C/2e1e2^—  ,  (2.5a) 

0=  1-—  M-{l-x)aB~\bL  .  (2.5b) 

Here  /  is  a  coupling  constant,  and  a,  8,  and  x  are  dynam¬ 
ical  parameters  which  we  define  below.  The  novel  feature 
of  Q  in  (2.5b)  is  that  the  charge  in  the  PSW  model  is  a 
linear  combination  of  the  mass  M,  baryon  number  B,  and 
isospin  Iz  of  a  nucleus,  and  hence  the  product  Q\Q2  in 
Eq.  (2.5a)  contains  cross  terms  among  M ,  B,  and  Iz.  For 
later  purposes  it  is  instructive  to  understand  how  the  ex¬ 
pression  for  Q  in  Eq.  (2.5b)  comes  about.  A  scalar  field 
can  couple  to  the  nucleons  in  the  nucleus  through  two 
natural  scalar  operators.  These  are  T'vix),  which  is  the 
trace  of  the  full  energy-momentum  tensor,  and  0^(jc), 
which  is  the  anomalous  trace  of  the  energy-momentum 
tensor.  The  latter  operator  is  proportional  to  the  diver¬ 
gence  of  the  dilatation  current  J  (x),  which  is  not  con¬ 
served  at  the  quantum  level  due  to  the  presence  of 
anomalies,  i.e., 

0,y‘(x)=v/llgap(x)]l0;u) ,  (2.6) 

where  g^ix)  is  the  metric  tensor.  The  matrix  elements 
of  T'y(x )  and  ©','.( x)  for  a  nuclear  state  I N )  are  related  as 
follows: 

<jv|r;:|tf>=M ,  (2.7a) 

(N\Qy\N)=M  —  (N\muiiu  +  mddd\N > 

s*Af  —  aB  —  \8IZ  +  ■  (2.7b) 

m„ 

Here  mu  (md)  is  the  mass  of  the  u  (d)  quark,  and  x  « 1 
is  a  parameter  which  represents  the  fraction  of  the  time 
that  the  operator  ( muiiu  +mddd )  contributes  to  the 
binding  energy  eB  rather  than  to  a  or  6.  Theoretical  ar¬ 
guments  suggest  that  the  term  proportional  to  x  can  be 
neglected,  and  if  we  temporarily  assume  that  5=0  as 
well,  then  { 0J!)  in  Eq.  (2.7b)  can  be  written  as 

(N\®';.\N)3*M  1-—  —  ,  (2.8) 

where  we  have  written  M=y.mH.  One  can  form  an 
infinite  number  of  charges  by  taking  linear  combinations 
of  (©!.}  and  <r;:>,  and  from  Eqs.  (2.7)  and  (2.8)  these 
charges  will  have  the  form 

Q{y)  =  M  1  -y  —  ,  (2.9a) 

where 

e=ow<0;:>+a7<r;:> ,  (2.9b) 
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and  where  the  “mixing  parameter"  y  is  given  by 


y~ 


a e 

oe+aT 


a 

mH 


Consider  now  the  interaction  of  test  masses  i  and  2  in  the 
presence  of  the  potential  that  would  arise  from  charges 
(2.9c)  Qa  ai>d  Qb  >n  E<l-  <2.9)  corresponding  to  mixing  parame¬ 
ters^  andy*: 


V(r)=-Gx$e 


-,/A.  m,m2 


i-ya 

B 

l  -y. 

B 

~~~ 

— ' 

M 

1 

l* 

-r/k.  m  iKI; 


*  — -V* 


1  -»> 


(2.10) 

At  this  stage,  we  can  reinstate  the  contribution  from  lz  by  letting  B/p— *B  /p  +  {6/2o)Iz/p  =  q.  The  potential  energy 
in  Eq.  (2.10)  can  then  be  rewritten  in  the  form 


V(r)=-Gx£, 


)-6^' 


m ,  m-, 


r  /k  . 


-r/k. 


+‘}\‘hlyae 


,  r/k  ,  r/A,  , 

2-  “+yte  *)]. 


(2.11) 


The  first  term  in  Eq.  (2.11)  is  composition  independent, 
and  hence  leads  to  the  usual  phenomenology  for  mine 
and/or  borehole  and  tower  measurements,  but  with  a 
modified  spatial  dependence.  The  remaining  terms  in 
(2.11)  will  also  contribute  to  composition-independent 
effects,  but  they  are  of  particular  interest  here  because  of 
the  composition-dependent  effects  to  which  they  lead. 
The  term  proportional  to  (q,  +q2 ),  which  arises  from  the 
interference  of  the  two  components  of  Q{y)  in  Eq.  (2.9a), 
will  lead  to  an  exponential  potential  whenever  ya  »  —yb. 
Since  ya  and  yb  can  be  fixed  by  an  appropriate  choice  of 
the  constants  ae  and  or,  a  model  based  on  the  PSW 
charge  can  lead  to  an  exponential  arising  from  scalar  ex¬ 
change. 

Another  example  of  how  an  exponential  potential 
could  arise  from  scalar  exchange  follows  from  the  model 
of  Halprin,  Barnhill,  and  Barr20  (HBB).  These  authors 
consider  a  weak  force  mediated  by  a  gauge-singlet  scalar 
field  x  which  has  no  direct  couplings  to  light  quarks  or 
leptons,  but  which  acquires  an  effective  coupling  by  mix¬ 
ing  with  the  usual  Higgs  field  <j>.  In  such  a  theory  the 
analog  of  the  charge  Q  in  Eq.  (2.S)  becomes 

Q«(B-omL)  ,  (2.12) 

where  B  and  L  denote  baryon  number  and  lepton  num¬ 
ber,  respectively.  This  theory  can  be  generalized  by  in¬ 
troducing  the  coupling  of  a  field  X  to  two  Higgs-doublet 
fields  $\  and  $2-  If  X  couples  to  the  combination 
<£,sina+$2cosa,  then  the  charge  Q  in  Eq.  (2.12)  becomes 
a  function  Q  =  Q(z[a])  of  the  mixing  parameter  a  such 
that 

Q(z)*B  +zL  ,  (2  |3) 

Q(z(a—0))<x B  —0.2L,  Q{z(a  =  n/2))^  B  +0.3L  . 

If  two  fields  xa  and  xb  exist,  which  to  lowest  order  couple 
to  and  $2  but  not  to  each  other,  then  the  resulting  po¬ 
tential  V'(r )  will  have  the  form 

~r/KU 

V(r)=f;Z-j—($+zJ)tf+zJ)2 

-  r/kh 

+fb~~r  (b+zbT){(b+z„t)2  ,  (2.14) 


where  b=B/p  and  t=L/p.  As  before  we  can  regroup 
the  terms  in  (2. 14)  to  give 


V'(r)  =  bxb2 


-  r/K  . 


fl 


■+/*■ 


-(-(Si T2  +  b2lt ) 


- r/kb 

r 

-  r/k , 


-r/kK 


fX^-—+flzbL- 


+/,/ 


1*2 


f2z2- 

J  a*a 


-  r/ktt 


-r/kk 


.  (2.15) 


Since  za  and  zb  can  have  opposite  signs,  the  term  propor¬ 
tional  to  (bj2  +  b2Tt )  can  give  rise  to  an  exponential  po¬ 
tential  for  appropriate  values  of  /2  and  fb.  We  note  that 
here,  as  in  the  PSW  model,  the  exponential  arises  from 
the  cross  terms  that  are  characteristic  of  scalar  exchange 
in  such  models. 


III.  PHENOMENOLOGY 
IN  THE  EXPONENTIAL  MODEL 


We  discuss  in  this  section  the  mechanisms  by  means  of 
which  the  exponential  potential  changes  the  phenomeno¬ 
logical  interpretation  of  the  existing  composition- 
dependent  experiments.  Since  the  present  status  of 
composition-independent  searches  for  non-Newtonian 
gravity  is  somewhat  unsettled,  the  implications  of  the  ex¬ 
ponential  model  for  such  experiments  will  be  considered 
elsewhere.  Starting  with  Eq.  (2.3)  the  sum  of  the 
Newtonian  and  non-Newtonian  contributions  can  be 
written  in  a  form  analogous  to  Eq.  (1.3),  but  with  G(r) 
now  having  the  form 


G(r)  =  Gjr(r)  =  G0 


A  A. 

Qj_ 

Qj 

k 

Mi 

m2 

-  r/k 


(3.1) 


Comparing  Eqs.  (1.3)  and  (3.1)  we  see  that  in  both  cases 
the  non-Newtonian  effects  vanish  for  r— «•  oo  as  expected. 
However,  for  the  exponential  force,  the  non-Newtonian 
effects  become  negligible  as  r~*0,  in  contrast  with  the 
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Yukawa  case  where  they  do  not.  It  follows  from  Eq.  (3.1) 
that  the  exponential  force  leads  to  a  suppression  of  the 
non-Newtonian  contributions  from  sources  located  at 
separations  r  «  X  from  the  detector.  Thus  a  laboratory 
experiment  utilizing  as  a  source  a  sphere  of  radius  R 
would  be  reduced  in  sensitivity  by  a  factor  of  order 
(R  /A.)2  compared  to  what  would  be  expected  from  a  Yu¬ 
kawa  potential.  For  typical  values  of  R  and  X,  e.g.,  R  =  1 
m  and  X  =  1000  m,  we  have  {R  /X)2  =  10'6.  Since  the 
O(eS)  correction  term  in  Eq.'  (2.1)  becomes  important 
when  (r/X)2  £8  (see  Appendix  A),  an  additional  Yukawa 
contribution  can  arise  whose  strength  is  of  order  /2(e8), 
and  which  does  not  vanish  as  r— »0.  A  reasonable  upper 
bound  on  8  is  aem/ir» 2  X  10“ \  which  is  the  largest  small 
parameter  that  naturally  appears  in  perturbative  theories. 
It  follows  that  laboratory  experiments  are  suppressed  by 
a  factor  of  either  (R  /X)2  or  6,  depending  on  the  values  of 
X  and  8  and  on  the  dimensions  of  the  source  (~~R).  Since 
constraints  on  a  possible  coupling  to  /,  come  primarily 
from  laboratory  experiments,21  there  are  no  significant 
constraints  on  the  strength  £,  of  such  a  coupling  for  X  5  10 
m  in  the  exponential  model. 

For  an  experiment  utilizing  a  large  source  {R/kZ  8), 
the  exponential  contribution  dominates.  Experiments 
with  large  man-made  sources  have  been  carried  out  by 
Thieberger22  and  by  Bennett,23  but  the  most  stringent 
tests  of  the  exponential  model  come  from  the  cliff  or  hill¬ 
side  experiments  of  Thieberger, 24  Adelberger  et  al.,25 
Boynton  et  al.,2t  Fitch  et  al.,21  and  Bizzeti  et  al. 28  To 
evaluate  the  intrinsic  strengths  7f(r,X)  of  the  various 
geophysical  sources  in  the  exponential  model,  one  can  ei¬ 
ther  directly  integrate  the  point  source  expression  given 
in  (2.4),  or  else  start  from  the  corresponding  results  for 
7y(r,X)  for  the  Yukawa  model  (if  known)  and  use 


3 F£(r,X)3*X 


d7y(r,X) 

ax 


(3.2) 


To  study  the  differences  between  the  exponential  and 
Yukawa  models  for  large  sources,  we  have  analyzed  the 
experiments  of  Adelberger  et  al. 25  (Seattle),  and  Boynton 
et  al. 26  (Mt.  Index),  which  have  the  greatest  intrinsic  sen¬ 
sitivities  among  the  hillside  experiments.  The  surface  ter¬ 
rain  surrounding  each  site  was  modeled  in  a  series  of 
grids  of  different  scales,  which  varied  in  size  from  SO 
m  X  50  m  to  500  m  X  500  m  depending  on  the  topography 
and  on  the  distance  from  the  experimental  site.  The  grids 
contained  5653  squares  at  the  Seattle  site,  and  5246 
squares  at  the  Index  site,  corresponding  to  a  total  area  of 
approximately  30  kmX30  km  surrounding  each  experi¬ 
ment.  For  each  square,  the  maximum  and  minimum 
elevations  were  recorded,  and  these  were  used  to  estimate 
the  error  arising  from  the  discretization  of  the  topogra¬ 
phy.  Our  results,  which  were  obtained  from  Eqs. 
(B15MB22),  are  shown  in  Figs.  1(a)  and  Kb).  These 
figures  exhibit  the  north  and  east  components  of  the 
source  strength  1 7|  =(  SfjLmh  +  &L* )' /2  as  a  function  of  X. 
In  the  Yukawa  model  | SFflndex)!/]  7(Seattle)|  varies  from 
2.6±0.7  (at  X=  100  m)  to  4.2±0.4  (at  X*  1060  m).  How¬ 
ever,  the  same  sources  appear  quite  different  in  the  ex¬ 
ponential  model,  as  we  see  from  Fig.  1(b).  For  this  model 


FIG.  !.  Field  strengths  for  the  Index  and  Seattle  sites,  in 
units  of  mgal  =10  'ms  2.  At  each  site,  the  north  and  east 
components  are  plotted  separately,  (a)  Yukawa  case  (7y(X)], 
and  (b)  exponential  case  [(7,.  (A.)]. 


|7(Index)|/|7(Seattle)|  is  larger,  varying  between  2.1±0.2 
(at  X  s*  2000  m)  and  6.2±0.2  (at  X»5970  m).  In  each 
model  the  direction  of  7  changes  as  a  function  of  X  for 
both  Seattle  and  Index,  but  this  variation  is  considerably 
greater  for  the  exponential  model  than  for  the  Yukawa. 

The  phenomenological  differences  between  the  Yukawa 
and  exponential  models  result  from  the  characteristic 
functional  forms  for  Gr{r)  and  GE(r)  in  Eqs.  (1.3)  and 
(3.1),  and  lead  to  different  weightings  for  the  contribu¬ 
tions  to  7  from  matter  at  a  given  distance  from  the  detec¬ 
tor.  We  note  specifically  from  (3.1)  that  the  maximum 
deviation  of  GE(r)  from  the  Newtonian  result  occurs  at 
r=2X,  whereas  for  the  Yukawa  model  it  is  at  r=0.  It 
follows  that  in  the  exponential  model  the  source  strength 
|7|  for  a  particular  site  with  a  large  matter  distribution 
at  ( r  >  «  2X,  may  be  enhanced  relative  to  that  at  another 
site  with  a  different  matter  distribution.  Moreover,  this 
effect  can  be  even  more  dramatic  in  the  variants  of  the 
exponential  model  which  we  now  discuss. 

Just  as  the  simple  exponential  model  is  an  interesting 
limiting  case  of  two  nearly  cancel  2  Yukawa  potentials, 
there  are  generalized  exponential  models  which  are  the 
limiting  cases  of  several  canceling  Yukawa  potentials. 
We  can  parametrize  these  in  the  form 


V'[n)(r)  =  (Gacmim2)^nqlq2 


.  -  r/K 


n  =0,1,2 .  (3.3) 
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G(r)  =  C‘"'(r) 

n  +  I 

=  GX  1+£,4i92  f  "~X  6  /X  ’  <3‘4) 

where  the  £„  are  appropriate  constants.  We  see  from 
(3.4)  that  for  n>l  the  deviation  of  Gul(r)  from  the 
Newtonian  value  Gx  vanishes  not  only  at  r=0  and 
r— *  qo,  but  also  at  r  =nk.  This  is  a  novel  feature  of  such 
an  interaction,  which  has  interesting  experimental  conse¬ 
quences^  Since  G'"'(r)  has  extrema  at  r=(l+/r 
± V  1  +  n  )A.,  and  changes  sign  in  going  from  r  <nk  to 
r  >nk,  it  follows  that  even  an  experiment  with  a  large 
geophysical  source  might  fail  to  detect  the  presence  of 
V'in\r),  due  to  fortuitous  cancellations  in  the  source  con¬ 
tributions  (see  below). 

It  is  clear  from  the  preceding  discussion  that  if  the 
non-Newtonian  force  were  described  by  the  potential  in 
Eq.  (3.3)  with  ti>  1,  our  view  of  the  current  experimental 
situation  could  be  quite  different  from  what  it  is  in  the 
conventional  Yukawa  formalism.  It  is  thus  interesting  to 
note  that  the  sum  of  a  small  number  of  primitive  Yukawa 
terms  can  in  fact  lead  to  an  expression  for  G(r )  having 
the  property  that  AC(r)=[6(r)- Gx  ]  vanishes  at  some 
intermediate  value  of  r.  As  an  example,  consider  the 
four-component  potential 

y'(r)—f2QiQ2 


r  r 

where  mu  =  m(  l— e/2),  mb=m(  1  +e/2), 

mc  =  m(l+»c)(l  +  e/2),  md  —  m  ( 1  +k)(  1  —  e /2),  with 
e,K  « 1.  Expanding  (3.5)  in  e  and  k,  and  setting  A  =  1  /m, 
we  find  for  the  force 

F(r)3-?/20,0j«(2-r/A.)e'rAA2  .  (3.6) 


10' 

10° 

io-' 
y 

10-2 

10"3 
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X 

FIG.  2.  Comparison  of  normalized  field  strengths 
Y  5  ly(A)|  =  \  Hk)\  /2-rrpG  y  k  of  the  Newtonian  (inverse-square), 
Y ukawa,  and  exponential  models,  along  with  the  n  =  I  model  in 
Eq.  (10)  and  the  four-component  model  of  Eq.  (12).  The 
strengths  are  plotted  versus  x=R/k  for  an  observer  on  the  axis 
of  a  cylindrical  source  (radius  =  length  =  R )  at  one  end  cap. 


It  follows  from  Eqs.  (3.5)  and  (3.6)  that  since 
F*(r  =  2A.)=0,  the  vanishing  of  F‘  at  intermediate  r  could 
conceivably  arise  from  a  relatively  simple  mass  spectrum. 
To  quantitatively  study  the  phenomenological  implica¬ 
tions  of  forces  which  vanish  at  intermediate  r,  consider 
V{  1 '  in  Eq.  (3.3),  and  the  four-component  potential  in  Eq. 
(3.5).  We  calculate  |7|  for  an  observer  on  the  axis  of  a 
cylindrical  source  (length  =  radius =/? )  at  one  end  cap, 
and  compare  these  resul's  to  the  corresponding  Newtoni¬ 
an,  Yukawa,  and  exponential  contributions  (see  Appen¬ 
dix  B  and  Fig.  2).  We  see  that  for  the  cylindrical  source 
all  of  the  generalized  exponential  models  are  suppressed 
at  small  x  =R  /k.  At  large  x,  however,  all  finite-range 
models  (including  the  Yukawa)  converge  to  the  same  lim¬ 
it.  This  implies  that  when  the  characteristic  size  is  much 
greater  than  k,  the  interpretation  of  large  scale  experi¬ 
ments,  such  as  the  Galileo  free-fall  experiments, 2’- 30  is 
essentially  model  independent.  Most  interestingly,  the 
n  —  1  results  show  that  in  some  models  an  observer  adja¬ 
cent  to  a  macroscopic  source  could  fail  to  detect  the  pres¬ 
ence  of  a  non-Newtonian  force  (in  this  case  for  A  se  2.9), 
even  though  F'"’(r)# 0.  The  fact  that  this  does  not 
occur  in  the  four-component  model  further  underscores 
the  sensitivity  of  the  final  results  to  both  the  nature  of  the 
source  and  to  the  detailed  model  of  F'(r). 

IV.  SUMMARY 

To  summarize,  we  have  shown  that  the  exponential  po¬ 
tentials  in  Eqs.  (2.3)  and  (3.3)  can  arise  from  simple  mod¬ 
els,  and  that  the  phenomenological  implications  of  such 
couplings  can  be  quite  different  from  those  of  the  conven¬ 
tional  Yukawa.  Although  the  exponential  model  cannot 
reconcile  the  existing  (conflicting)  data,  it  is  clear  that  the 
comparison  of  experimental  results  using  different 
sources  (both  geophysical  and  laboratory)  is  highly  model 
dependent.  Hence  by  analyzing  experimental  data  in 
terms  of  both  a  Yukawa  and  an  exponential,  we  are  in 
effect  exploring  the  implications  of  a  much  broader  class 
of  interesting  models. 
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APPENDIX  A:  EXPANSION 
OF  TWO  YUKAWA  POTENTIALS 

We  present  in  this  appendix  some  additional  details  of 
the  exponential  model.  1  As  noted  in  Sec.  II,  the  ex¬ 
ponential  model  arises  as  the  limiting  case  of  two  nearly 
canceling  Yukawa  potentials,  when  the  masses  and  cou¬ 
pling  constants  satisfy  Eq.  (2.1).  To  clarify  the  relation 
between  Eq.  (2.1)  and  the  exponential  model,  we  write  the 
non-Newtonian  potential  energy  V'(r)  arising  from  the 
sum  of  two  Yukawa  contributions  in  the  form 
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V\r)  = 


Gxm,m2  Q , 

Qj 

£a=$(l+e6/2), 

$=‘$.+£*>/2  . 

'  Mi 

Pi 

£*  =  £(1-68/2), 

X(-£ue  “+&«  *>. 


A.a=A.(l+e/2),  A  =  (Aa  +  A.t)/2  , 
kh  =  Ml-e/2),  e  =  (A.a  —kb  )/k  , 
9,=C?,/M,.  *'*1.2. 


rhere  by  hypothesis  the  two  contributions  are  assumed  g  —q  ,  =  (  2 

o  interfere  destructively.  To  incorporate  the  content  of 

!q.  (2.1)  we  first  simplify  our  notation  by  writing  Assuming  8,e  «  1  as  before,  we  find,  upon  expanding  Eq. 
)(£)—►«  and  0(c8)-*e6,  and  we  then  define  (Al), 


$Gxmtm2qtq2)  ":K*(r)  =-€fi- 


e",/A  .  cJ8 

+~r 


_ r_  e  r/K  e3  .  r  r2  e  r/k 

2k  k  4  A.  6A.2  A. 


+  - r— 

2k  2  k2  24k3 


e  'M  __  e5  _  2r  r2  ri  r 4 

_X  76  1  T+aT  6A.J  +  120A4 


This  expansion  can  be  reexpressed  in  the  form 

x  „  -  r /X  at  [ 


=  2  s„{k)=  — 

»i  —I 


7  •  <A4) 

a  —  - 1  A 


Since  6  and  e  are  assumed  small,  we  need  only  consider 
the  dominant  contribution  to  each  an  in  Eq.  (A4)  above. 
For  n  >  0,  we  then  obtain 


n  even  , 


A*3 


2"(n  +1)! 

+  2 

— — — ,  n  odd  . 
2"  +  « ! 


APPENDIX  B:  CLOSED-FORM  SOLUTIONS 

In  this  appendix,  we  present  closed-form  solutions  to 
the  non-Newtonian  fields  for  both  the  Yukawa  and  ex¬ 
ponential  models.  '1  We  begin  by  defining  these  quanti¬ 
ties  in  terms  of  the  differential  acceleration  between  two 
test  masses.  From  Eq.  (Al),  we  have 


Aa=-— 7P(r) 
m,m2 


-£eqsAq0[8Sfy(r,A.)  +  :7£(r,A.)]  , 


where  S  and  D  refer  to  the  source  and  detector,  respec¬ 
tively,  and’1 


We  can  show  that  only  the  n  =  —  1  and  n=0  contribu¬ 
tions  need  be  considered  for  present  phenomenological 
purposes.  We  do  this  by  noting  that  for  n=even  a„  will 
become  comparable  to  a0  when 
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*<>  =  *„  =  -^-e*P  - - — 1 —  •  <A7) 

Even  for  e=  ^  and  n= 2,  at  the  value  of  r/k  where 
r0=s„,  exp(  )  is  on  the  order  of  10-22;  for  e—  ^  exp(  )  is 
an  the  order  of  10  22\  The  analogous  terms  in  the  ex¬ 
ponential  decrease  even  more  rapidly  for  w=odd.  We 
conclude  that  the  phenomenology  will  be  dominated  by 
the  n  =  —  1  term  for  r/k  £8,  and  by  the  n=0  term  for  all 
other  values  of  r/k  for  which  V(r)  is  non-negligible. 
V(r)  is  then  given  to  the  required  accuracy  by 

-  r/X  r/X 

V’(r)  =  —  $eG, m\m1q\q2  8 - H - - —  (A8) 
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=  f  dV— - T  e ~!  .  (B3) 

Jbody  X2  I r — r  | 

Using  Eq.  (B3),  we  can  obtain  closed-form  expressions  for 
7K(r,A.)  when  closed-form  expressions  for  S7y(r,A.)  al¬ 
ready  exist.  In  the  examples  we  consider  below,  the  mass 
distribution  is  assumed  to  be  chemically  homogeneous 
and  of  constant  density  p. 

(i)  Point  mass.  We  have 


SUr,A.)=rGaM  1  +  -J1  — 


where  M  is  the  source  mass,  and  r  is  the  distance  from 
the  source  to  the  detector.  Hence 
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yE(r,k)=iGxM- 


—  r/X 
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A2 


(BS) 


in  agreement  with  Eq.  (2.3).  In  comparing  Eqs.  (B4)  and 
(B5),  we  note  that  for  r«k  the  exponential  force  is 
suppressed  relative  to  the  Yukawa  force,  since  the  point 
mass  is  always  (effectively)  a  distance  A  away  from  the 
source.  This  is  the  origin  of  the  short -distance  suppres¬ 
sion  of  the  exponential  force  which  we  have  discussed 
previously. 

00  Spherical  mass  distribution.  We  consider  the  case 
of  a  test  object  located  a  distance  r  from  the  center  of 
mass  of  a  spherical  mass  distribution  of  uniform  density 
with  mass  M  and  radius  R.  For  r  >  R  we  have1 
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(Bll) 

(Hi)  Cylindrical  mass  distribution.  We  consider  here 
the  case  of  a  test  object  located  along  the  symmetry  axis 
of  a  cylinder  of  uniform  density  p,  radius  R,  and  length 
L.  We  also  choose  z=0  to  occur  at  one  end  and  define  z 
to  increase  away  from  that  end.  We  then  have 


<b(x)=  — r(x  coshx  —  smhx)  . 
x 


Since 


A  ^/-k)  .  34><  R  /A )  -  3  sinh(  R  /k  > 
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we  find 
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For  purposes  of  comparison,  the  gravitational  accelera¬ 
tion  g  for  this  configuration  is  given  by 


~tGxM 
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(BIO) 


g(r)~  -~z27rpGec(V/R2+z'2  —  lz'l)lyZt  +  l'  .  (B14) 

(iv)  Parallelipiped  mass  distribution.  We  finally  consid¬ 
er  the  case  for  a  parallelipiped  whose  edges  are  aligned 
along  the  x,  y,  and  z  axes,  and  whose  extent  is  given  by 
vertices  located  at  (x0,y0,z0)  and  (xx,yx,Z\),  with 
x0  <x,,  yn  <yx,  z0  <z,.  (The  more  general  case  can  be 
obtained  from  the  results  presented  here  by  three- 
dimensional  rotations  about  the  appropriate  axes.)  Then 


5MA)-  ”pGk  2  (  - 1 )'  + ;  +  Ax, , Ay„  Az, ; A ) + yd^-l Ay,, Az,, Ax,; A ) +zd>os( Az, , Ax, , Ay, ; A ) ]  , 

ij,  k  ~  0  v 


where  Ax,  — x  — x,, . . . ,  and  4»gS(x,y,z;A)  is  given  by 

^QS^x,y,z  -,k)  =  <f>5(x  /k,y  /k\eo)  +  ib/x /k,z /k;6'0)  (y  <0,  z  <0) 

—  /k,y /k;d0)  —  <l>;(x /k,z /k;6'0)  ( y  <0,  z  >0) 
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and  where 


do=f~Oo<  (B18) 
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<Ma,/3)  =  <Ma,0;O)  .  (B20) 

*;(a.ft0)  =  2*(<a,0)-*,<a,/3;0)  ,  (B21) 

The  solution  for  7f(X)  for  a  parallelipiped  can  be  ob¬ 
tained  trivially  from  this  equation,  by  defining 


^f(a,/3;0)=  —  f  w/2J0'(  1  +  V'a2+/3^scc2<n 

IT  J  6 

X  exp(  -  vV+02 secJ0* )  ,  (B22) 

and  letting  ^s(a,/3;0)-»^f(a,/3;0)  everywhere  in  Eqs. 
(B16MB21). 
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Non-Newtonian  Gravity  and  New  Weak  Forces: 
an  Index  of  Measurements  and  Theory 


E.  Fischbach,  G.  T.  Gillies,  D.  E.  Krause,  J.  G.  Schwan  and  C.  Talmadge 


Abstract.  The  precise  measurement  of  weak  effects  plays  a  pivotal  role  in  metrology  and  in  the  determination  of  the 
fundamental  constants.  Hence,  the  possibility  of  new  weak  forces,  and  the  related  question  of  non-Newtonian  beha¬ 
viour  of  the  gravitational  force,  have  been  of  special  interest  to  both  measurement  scientists  and  those  involved  in  pre¬ 
cise  tests  of  physical  laws.  To  date  there  is  no  compelling  evidence  for  any  deviations  from  the  predictions  of 
Newtonian  gravity  in  die  nonrelativistic  weak-field  regime.  A  significant  literature  on  this  question  has  developed  over 
the  past  few  years,  and  a  host  of  experiments  and  theoretical  scenarios  have  been  discussed.  Moreover,  a  very  close 
relationship  exists  between  the  experimental  methodologies  used  to  determine  the  absolute  value  of  the  Newtonian  gra¬ 
vitational  constant  G,  and  those  employed  in  searches  for  new  weak  forces  and  fra  breakdowns  in  the  inverse-square 
law  of  gravity.  We  have  therefore  prepared  a  new  index  of  measurements  of  such  effects,  using  the  original  bibliogra¬ 
phic  work  of  Gillies  as  a  starting  point,  but  also  including  citations  to  the  appropriate  theoretical  papers  in  the  field. 
The  focus  of  the  present  version  of  the  index  is  then  studies  of  the  "fifth  force",  measurements  of  gravitational  effects 
on  antimatter,  searches  for  a  spin-component  in  the  gravitational  force,  and  related  phenomena. 


1.  Introduction 


During  the  past  several  years  interest  in  die  possibility  of 
deviations  from  die  predictions  of  Newtonian  gravity  has 
increased  significantly,  motivated  by  both  theoretical  and 
experimental  considerations.  On  die  theoretical  side,  a 
number  of  workers,  most  notably  Fujii  [1]  and  Scherk  [2], 
stimulated  interest  in  this  question  by  demonstrating  that 
various  models  suggest  the  existence  of  new  weak  inter¬ 
mediate-range  forces  coexisting  with  gravity.  These  forces 
share  with  gravity  the  property  of  acting  over  macroscopic 
distances  but,  unlike  gravity,  their  effects  are  negligible 
beyond  some  characteristic  distance  k  (called  the  range) 
from  a  source.  As  we  discuss  in  more  detail  below,  the 
combined  effect  of  such  a  finite-range  force  and  gravity  is 
a  net  gravity-like  force  whose  overall  strength  depends  on 
the  separation  r  -  )#■ ;  -  r  I  of  two  test  masses  i  and  j,  and 
thus  leads  to  deviations  from  the  expected  1/r2  force  law. 

Stimulated  by  these  ideas,  experimental  studies  were 
undertaken  in  the  1970s  and  1980s  to  search  for  deviations 
from  the  1/r3  law.  A  detailed  history  of  these  efforts  is 
contained  in  a  number  of  recent  reviews  [3, 4, 5]  to  which 
we  refer  the  interested  reader.  The  accumulation  of  a  num- 
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her  of  experimental  and  theoretical  results  eventually  led 
Fischbach  et  al.  [6]  to  reanalyze  the  classic  experiment  of 
Edtvfts,  Pekir,  and  Fekete  [7]  (EPF)  which  compared  the 
accelerations  of  different  pairs  of  materials  to  the  Earth. 
Fischbach  et  al.  observed  that  the  experimental  results  of 
EPF  could  be  interpreted  as  suggesting  the  existence  of  a 
new  weak  intermediate-range  force,  whose  presence 
would  lead  to  apparent  deviations  horn  the  predictions  of 
Newtonian  gravity. 

Since  the  publication  of  [6]  a  large  number  of  experi¬ 
mental  and  theoretical  papers  have  been  written  which 
explore  in  various  ways  the  possible  existence  of  new 
weak  gravity-like  forces.  At  present  there  is  no  compelling 
evidence  for  any  deviations  from  the  predictions  of 
Newtonian  gravity  in  the  nonrelativistic  regime  (i.e.  for 
systems  where  general  relativistic  effects  are  negligible). 
The  present  bibliography  is  an  attempt  to  collect  in  a 
single  place  all  the  relevant  literature  on  such  forces.  To 
understand  the  relationships  among  the  various  papers,  as 
well  as  the  criteria  we  have  used  in  compiling  this  biblio¬ 
graphy,  the  following  simplified  phenomenological  intro¬ 
duction  may  be  useful.  The  presence  of  a  single  new  field 
of  mass  m  gives  rise  to  an  additional  Yukawa  potential 
proportional  to  exp  (-  rfk)  /r,  where  X  -  h/mc.  This  poten¬ 
tial  modifies  the  Newtonian  interaction  so  that  the  new 
potential  describing  the  interaction  of  masses  m,  and  m. 
assumes  the  form 

Y(r)=-Gn  m‘  m>  (l+ae-'M).  (I) 
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Here  is  the  Newtonian  constant  of  gravity  (in  die  limit 
r  — »  oo),  and  a  is  a  constant  which  characterizes  the 
strength  of  the  new  interaction  relative  to  gravity. 
Differentiating  (1)  we  find  for  the  corresponding  force, 


P(r)=  K<r)  =  -  "•  ^  '  [l+a(l+rM)e-^] 

r* 


r 

G(r)=G„[l  +  a(l+r/A)c-"*]. 


We  see  from  (2)  dial  the  effect  of  the  new  interaction  is  to 
replace  the  Newtonian  constant  G^  by  the  function  G(r) 
in  (2b).  It  follows  that  tests  of  the  constancy  of  G(r)  as  a 
function  of  r,  i.e.  tests  of  deviations  from  a  simple  1/r2  force 
law,  are  also  direct  tests  of  die  presence  of  new  couplings. 

In  the  form  of  (2),  which  is  that  suggested  by  Fujii  [1] 
and  others,  a  is  a  universal  constant  determined  by  the 
coupling  strength  of  some  new  quantum  to  matter.  The 
work  of  Fischbach  et  al.  [6]  focuses  attention  on  the  fact 
that  in  many  theories  a  is  not  a  universal  constant,  but 
depends  instead  on  the  compositions  of  the  test  masses  i 
and  j.  In  the  original  "fifth  force"  model  of  [6],  a  is 
given  by 


where  ^  -F/Goa  mf,  expresses  the  strength  of  the  new 
interaction  in  terms  of  a  new  constant  /,  which  is  the  ana¬ 
log  for  this  interaction  of  what  the  electric  charge  e  is  for 
electromagnetism.  Here  B,  denotes  the  baryon  number 
(the  number  of  neutrons  and  protons)  in  the  sample  i,  and 
Hi  =  mJmH ,  where  mH  a  m  (,//')  is  the  mass  of  atomic 
hydrogen.  Since  (/?,//!,)  varies  from  one  material  to 
another,  a (  is  a  function  of  the  compositions  of  the  inter¬ 
acting  materials.  It  follows  that  the  accelerations  of  two 
samples  j  and  j'  towards  a  common  source  i  depend  on  the 
(generally  unequal)  constants  a y  and  a,/.  If  i  denotes  the 
Earth,  then  the  accelerations  of  j  and  /  towards  the  Earth 
depend  on  their  compositions,  and  this  is  what  the  Eotvos 
experiments  set  out  to  measure. 

It  follows  from  the  preceding  discussion  that  the  pre¬ 
sence  of  a  new  intermediate-range  weak  force  can  be 
detected  either  through  the  modification  of  the  usual 
inverse-square  law  for  the  net  force,  and/or  through  a 
composition-dependence  of  the  net  acceleration.  (The  lat¬ 
ter  effect  is  often  referred  to  as  a  violation  of  the  Weak 
Equivalence  Principle  or  of  the  universality  of  free-fall). 
Although  in  principle  both  of  these  effects  are  always  pre¬ 
sent,  in  practice  experiments  are  usually  designed  to  iso¬ 
late  one  or  the  other  of  these  manifestations  of  non- 
Newtonian  gravity.  When  appropriate  we  have  therefore 
classified  the  papers  in  this  bibliography  according  to 
whether  they  deal  primarily  with  tests  for  composition- 
dependent  effects  (denoted  by  CD),  or  with  tests  of  the 
inverse-square  law  (IS).  Forces  which  give  rise  to  compo¬ 
sition-dependent  effects  in  bulk  matter  will  generally 
affect  matter  and  antimatter  differently,  and  references  on 
antimatter  are  denoted  by  AM.  In  this  category  we  have 
also  included,  for  obvious  reasons,  papers  dealing  with 


the  1C-  K°  system.  New  forces  can  also  manifest  them¬ 
selves  via  spin-dependent  couplings,  and  papers  on  this 
topic  are  denoted  by  SD.  We  have  classified  papers 
depending  on  whether  they  are  primarily  experimental 
(E),  theoretical  (T),  or  phenomenological  (P),  the  final 
label  being  one  we  use  to  characterize  theoretical  papers 
analyzing  data.  Where  a  paper  quotes  a  new  experimen¬ 
tal  result  or  a  theoretical  limit,  we  have  denoted  this  by  a 
+  sign.  Summaries  of  the  existing  data,  along  with  des¬ 
criptions  of  the  techniques  that  have  been  used,  are  given 
in  [4]  and  [3]. 

In  addition  to  these  primary  classifications,  we  have 
also  noted  whether  an  experiment  utilizes  a  laboratory 
source  (L),  a  geophysical  source  -  e.g.  a  cliff  -  (G),  a  lake 
(LK),  or  an  astrophysical  source  (A).  For  theoretical 
papers  presenting  models  of  non-Newtonian  effects,  we 
have  denoted  by  AG  those  in  which  the  model  is  essential¬ 
ly  an  alternative  theory  of  gravity,  and  by  FI  those  for 
which  the  starting  point  is  a  new  fundamental  interaction 
such  as  the  "fifth  force".  In  addition,  theoretical  or  experi¬ 
mental  papers  dealing  with  high-energy  or  elementary  par¬ 
ticle  systems  are  denoted  by  H.  Finally,  REV  denotes  a 
review  paper,  and  INT  indicates  an  introductory  or  ele¬ 
mentary  exposition.  We  recognize  that  these  distinctions 
are  sometimes  subtle  and  not  always  well-defined,  and 
hence  the  labels  attached  to  each  paper  are  only  meant  to 
supply  a  general  indication  of  its  contents. 

We  have  been  guided  by  a  number  of  criteria  in 
selecting  the  papers  included  in  this  bibliography.  As 
explained  above,  our  focus  has  been  on  the  search  for  new 
weak  forces,  particularly  the  hypothesized  "fifth"  force 
first  proposed  in  1986.  We  have  thus  included  virtually  all 
papers  published  on  this  subject  since  January  1986.  In 
selecting  papers  published  prior  to  1986,  we  have  limited 
ourselves  to  those  which  can  be  viewed  as  the  direct  ante¬ 
cedents  of  die  current  searches  for  non-Newtonian  gravity, 
such  as  the  seminal  papers  by  Fujii.  With  very  few  excep¬ 
tions  we  have  thus  excluded  papers  already  cited  in  the 
comprehensive  bibliography  by  Gillies  [8]  dealing  with 
the  dependence  of  G  on  various  external  influences  (tem¬ 
perature,  shielding,  etc.)  which  can  also  be  interpreted  as 
searches  for  non-Newtonian  gravity.  Moreover,  given  the 
focus  of  our  effort,  we  have  chosen  not  to  include  papers 
aimed  primarily  at  rede  terminations  of  the  absolute  value 
of  G,  or  the  other  classical  effects  discussed  by  Gillies,  as 
these  will  be  the  topic  of  a  separate  updated  report.  It  is 
thus  our  intention  that  the  present  work  be  viewed  as  com¬ 
plementing  and  updating  that  of  Gillies  which,  coinciden¬ 
tally,  was  completed  just  as  interest  in  the  "fifth  force" 
was  beginning.  We  have  also  excluded  the  area  of  gravita¬ 
tional  effects  on  quantum  systems,  most  notably  the 
important  work  of  Colella,  Overhauser  and  Wemer  [9], 
because  the  extensive  literature  that  has  grown  up  in  this 
area  genuinely  warrants  an  index  of  its  own.  Finally,  we 
have  not  included  papers  dealing  with  the  deviations  from 
Newtonian  gravity  implied  by  general  relativity,  which  are 
covered  in  excellent  reviews  by  Will  [10],  or  papers  rela¬ 
ting  to  the  time  variation  of  G. 

We  acknowledge  that  our  choices  for  the  papers  to  be 
included  in  this  bibliography  are  necessarily  somewhat 
arbitrary.  We  also  recognize  that  there  are  papers  which 
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should  properly  be  included,  but  of  which  we  are  simply 
unaware.  We  therefore  apologize  in  advance  to  all  our  col¬ 
leagues  whose  papers  we  have  inadvertently  omitted.  It  is 
our  intention  to  update  this  listing,  and  we  therefore 
request  that  any  suggestions  for  papers  to  be  included  be 
brought  to  the  attention  of  the  authors.  Finally  we  wish  to 
thank  Sam  Aronson  and  Lisa  Schwan  for  their  help  in 
compiling  this  bibliography. 
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Key  to  Paper  Categories 

General  categories  Qualifiers 


CD 

Composition-Dependent 

E 

Experimental 

SD 

Spin-Dependent 

P 

Phenomenological 

IS 

Inverse-Square 
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Theoretical 

AM 

Antimatter 
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Laboratory 

FI 

Fundamental  Interactions 
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Geophysical 

AG 

Alternative  Gravity 
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Astrophysica) 
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Introductory 
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High  energy 
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Addendum 

Since  submitting  the  original  version  of  our  biblio¬ 
graphy,  we  have  received  numerous  suggestions  for 
additional  papers  to  be  included,  and  most  of  these  have 
been  incorporated  into  the  preceding  text.  However,  we 
have  also  decided  to  include  in  proof  a  number  of  papers 
related  to  spin-dependent  effects,  and  these  are  listed 
below.  Although  they  could  not  be  cross-referenced  in 
the  index,  we  hope  that  their  inclusion  here  will  serve  to 
make  this  bibliography  more  useful. 

1.  Aleksandrov  E.  B„  Ansel’m  A.  A.,  Pavlov  Yu.  V., 
Umarkhodzhaev  R.  M.,  "A  restriction  on  the 
existence  of  a  new  type  of  fundamental  interaction 
(the  anon'  long-range  interaction)  in  an  experiment 
on  spin  precession  of  mercury  nuclei",  Soviet 
Physics  JETP ,  1983,  58,  No.  6.  1103-1107. 
[Translation  of  ZJi.  Eksp.  Teor.  Fiz.,  1983,  85,  1899- 
1906].  Topics:  SD,E,+ 

2.  Ansel' m  A.  A.,  Uraltsev  N.  G.,  "Long-range  an  on' 
field  in  the  radio  frequency  band".  Physics  Letters, 
1982, 116B,  No.  2, 3. 161-164.  Topics:  SD,P 

3.  Ansel’m  A.  A..  "Possible  new  long-range  inter¬ 
action  and  methods  for  detecting  it”,  JETP  Letters, 
1982,  36,  No.  2,  55-59.  [Translation  of  Pis'ma  Zh. 
Eksp.  Teor.  Fiz.,  1982,  36,  No.  2,  46-49.]  Topics: 
SDJ> 

4.  Ansel’m  A.  A.,  Neronov  Yu.  I.,  "Restrictions  on  the 
existence  of  spin-spin  coupling  of  nonelectro- 
magnetic  origin  in  experiments  on  the  measurement 
of  the  gyromagnetic  ratios  of  the  proton  and 
deuteron",  Soviet  Physics  JETP,  1985,  61,  No.  6, 
1154-1155.  [Translation  of  Zh.  Eksp.  Teor.  Fiz., 
1985, 88. 1946-1949.]  Topics:  SD,E,+ 

5.  Hari  Dass  N.  D.,  "Test  for  C,  P,  and  T 
nonconservation  in  gravitation".  Physical  Review 
Utters,  1976, 36,  No.  8. 393-395.  Topics:  SD,FI»P 

6.  Hari  Dass  N.  D.,  "A  new  spin  test  for  the 
equivalence  principle".  General  Relativity  and 
Gravitation,  1977,  8,  No.  2,  89-93.  Topics: 
SDW 

7.  Hari  Dass  N.  D.,  "Experimental  tests  for  some 
quantum  effects  in  gravitation".  Annals  of  Physics 
(New  York),  1977, 107,  337-359.  Topics:  SD,FI,T 

8.  Leitner  J.,  Okubo  S..  "Parity,  charge  conjugation, 
and  time  reversal  in  the  gravitational  interaction". 
Physical  Review,  1964.  136,  No.  5B,  B1542-B1546. 
Topics:  SDJFI,P 

9.  Morgan  T.  A.,  Peres  A.,  "Direct  test  for  the  strong 
equivalence  principle".  Physical  Review  Letters, 
1962, 9,  No.  2, 79-80.  Topics:  SD,AG,P 

10.  Naik  P.  C.,  Pradhan  T.,  "Long-range  interaction  be¬ 
tween  spins".  Journal  of  Physics  A,  1981,  14,  2795- 
2805.  Topics:  SD,FI,T 

11.  O’Connell  R.  F.,  "Spin,  rotation,  and  C,  P,  and  T 
effects  in  the  gravitational  interaction  and  related 
experiments".  In  Experimental  Gravitation,  Proceed¬ 
ings  of  the  International  School  of  Physics,  Enrico 


Fermi,  Course  LVI,  (Edited  by  B.  Bertotti),  New 
York,  Academic  Press,  1974,  496-514.  Topics: 
SDJT.T 

12.  Peres  A.,  "Test  of  the  equivalence  principle  for 
particles  with  spin".  Physical  Review  D,  1978,  18, 
No.  8, 2739-2740.  Topics:  SD,AG,P 

13.  Pradhan  T.,  Malik  R.  P.,  Naik  P.  C.,  "The  fifth  inter¬ 
action:  universal  long-range  force  between  spins”, 
Pramana,  1985,  24,  No.  1,  2,  77-94.  Topics: 
SDJFI.T 

14.  Stueckelberg  E.  C.  G.,  “A  possible  new  type  of  spin- 
spin  interaction”,  Physical  Review,  1948,  73,  No.  7, 
808.  Topics:  SDJ* 

15.  Vasil’ev  B.  V.,  "Concerning  the  gravitational 
moment  of  the  proton”,  JETP  Utters,  1969,  9,  175- 
177.  [Translation  of  ZhETF  Pis.  Red.,  1969, 9,  No.  5, 
299-301.]  Topics:  SD,E,+ 

16.  Velyukhov  G.  E.,  "Searches  for  the  gravitational 
moment  of  the  proton",  JETP  Utters,  1968,  8,  229- 
232.  [Translation  of  ZhETF  Pis.  Red.,  1968, 8,  No.  7. 
372-375.]  Topics:  SD,E,+ 

17.  Young  B.  A.,  "Search  for  a  gravity  shift  in  the  proton 
Larmor  frequency",  Physical  Review  Utters,  1969, 
22.  No.  26,  1445-1446.  Topics:  SD,E,+ 

In  addition,  there  are  several  very  recent  papers  which 

we  include  here  for  the  sake  of  completeness. 

1.  Atiya  M.  S.,  Chiang  I.-H.,  Frank  J.  S..  Haggerty  J.  S„ 
Ito  M.  M.,  Kycia  T.  F.,  Li  K.  K.,  Littenberg  L.  S.. 
Stevens  A.  J.,  Sambamurti  A.,  Strand  R.  C.,  Louis 
W.  C.,  Akerib  D.  S.,  Marlow  D.  R.,  Meyers  P.  D., 
Selen  M.  A.,  Shoemaker  F.  C.,  Smith  A.  J.  S., 
Blackmore  E.  W.,  Bryman  D.  A.,  Felawka  L., 
Kitching  P.,  Konaka  A.,  Kuno  Y.,  Macdonald  J.  A., 
Numao  T.,  Padley  P.,  Poutissou  J.-M.,  Poutissou  R., 
Roy  J.,  Turcot  A.  S.,  "Search  for  the  decay  it0  -*  y  + 
X",  Physical  Review  Utters,  1992,  69,  No.  5,  733- 
736.  Topics:  AM,HJE,+ 

2.  Hughes  R.  J.,  ’Tests  of  gravity  and  CPT  symmetry 
with  trapped  antimatter”.  In  Intersections  Between 
Particle  and  Nuclear  Physics,  AIP  Conference 
Proceedings  243,  (Edited  by  W.  T.  H.  van  Oers), 
New  York,  American  Institute  of  Physics,  1992,  31 1- 
313.  Topics:  AMJ* 

3.  Smith  G.,  Adelberger  E.  G.,  Heckel  B.  R.,  Su  Y.,  ”A 
test  of  the  equivalence  principle  for  ordinary  matter 
falling  toward  dark  matter",  1992,  unpublished, 
10  pp.  Topics:  CD,A,+ 

4.  Carusotto  S..  Cavasinni  V.,  Mordacci  A.,  Penone  F„ 
Polacco  E.,  lacopini  E.,  Stefanini  G.,  "Test  of  g 
universality  with  a  Galileo  type  experiment”. 
Physical  Review  Utters,  1992,  69,  No.  12,  1722- 
1725.  Topics:  CD^P 
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Roland  Eotvoe'  classic  experiment  concerning  the  proportionality  of  inertial  and 
gravitating  masses,  performed  at  first  in  1889,  has  again  become  the  focus  of  scientific 
interest  in  the  1980's,  due  to  the  possibility  of  the  existence  of  a  Fifth  Force,  as  proposed 
by  Fischbach  and  co workers.  The  publication  of  Eotvoe,  Pekir  and  Fekete  omitted  various 
details  of  their  experiment  which  may  be  relevant  for  the  re-interpretation  of  their  results. 
The  aim  of  this  report  is  to  fill  in  some  of  these  details,  and  to  discuss  the  impact  of  the 
Eotvoe  experiment  on  modern  research. 


“Are  longs,  vita  brevis” 

Inspired  by  the  beauty  of  the  Newtonian  system,  Baron  Roland  von  Eotvos 
experimentally  investigated  the  proportionality  of  inertial  and  gravitating  masses  in 
1889,  and  reported  his  results  in  the  Proceedings  of  the  Hungarian  Academy  in  1890 
[1],  In  this  work  he  improved  Bessell’s  accuracy  1/60  000  to  1/20  000  000.  This  was 
a  short  report  of  3  pages.  Inspired  by  this  achievement,  the  Royal  Scientific  Society 
of  Gottingen  in  1906  offered  a  prize  (see  Appendix  I)  for  the  following  task: 

“ A  very  sensitive  method  was  given  by  Eotvos  to  make  a  comparison 
between  the  inertia  and  gravity  of  matter.  Considering  this  and  the  new  develop¬ 
ment  of  electrodynamics  as  well  as  the  discovery  of  radioactive  substances,  Newton ’s 
law  concerning  the  proportionality  of  inertia  and  gravitation  is  to  be  proved  as  ex¬ 
tensively  as  possible.” 

Eotvos  began  a  series  of  investigations  with  his  co-workers  Pekar  and  Fekete  in 
the  years  1906-1909.  This  included  data  taking  through  approximately  4000  hours. 
Eotvos  personally  reported  his  results  at  the  16th  International  Geodesic  Conference 
in  London  in  1909  [2],  quoting  an  achieved  accuracy  of  1/100  000  000.  The  complete 
work  of  Eotvos,  Pekar,  Fekete  was  submitted  to  the  Beneke  Foundation  in  1909  [3]. 
Its  motto  was  “Ars  longs,  vita  brevis ”  (the  art  lasts  long,  life  lasts  short),  which 
is  indeed  a  true  characterization  of  the  fate  of  Eotvos’  work.  The  evaluation  of 
C.  Runge,  Dean  of  the  Faculty  of  Science  in  Gottingen  [4]  says  that  Eotvoe  has 
quoted  an  accuracy  of  1/200  000  000,  but  since  the  submitted  text  does  not  include 
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the  real  theoretical  discussion  of  the  observational  data,  the  Faculty  recommends 
only  a  reduced  prize  for  this  work  (3400  German  Marks  instead  of  4500  Marks).  In 
Appendix  I  we  reprint  the  text  of  this  evaluation. 

Shortly  thereafter  the  First  World  War  came.  Roland  Eotvos  died  in  1919,  and 
the  detailed  description  of  the  experiments  performed  in  1906-1909  was  published 
by  his  assistants  Pekar  and  Fekete  only  in  1922  [5].  This  is  the  text  known,  cited, 
and  translated  by  the  international  scientific  community  worldwide.  When  the 
collected  works  of  Roland  Eotvos  were  published  by  the  Hungarian  Academy  of 
Sciences  [6],  the  editor  (Eotvos’  former  student  P.  Seienyi)  included  some  additions 
in  parentheses  [...]  in  the  reprinted  Eotvos- Pekar-Fekete  paper  [5],  taken  from  the 
original  manuscript  [3].  The  original  Beneke-prize  manuscript  was  lost  somewhere 
in  the  hands  of  the  heirs  to  Pekar  and  Seienyi.  The  more  complete  text  taken  from 
the  Volume  [6]  has  been  reprinted  in  English  in  Budapest  in  1963  [5]. 

The  Eotvos  experiment  was  repeated  by  J.  Renner  (a  former  student  of 
Eotvos,  and  a  physics  teacher  in  the  famous  Lutheran  High  School  in  Budapest 
where  among  others  J.  von  Neumann  and  E.  P.  Wigner  studied).  The  results  of  Ren¬ 
ner’s  experiment  were  published  in  Hungarian  [7]  (with  a  German  abstract,  reprint¬ 
ed  in  Appendix  II).  Renner  claimed  an  empirical  accuracy  of  1/2  000  000  000  to 
1/5  000  000  000. 


About  Dicke’s  criticisms 

Acknowledging  the  basic  role  played  by  the  connection  between  inertial  mass 
and  gravitating  mass  in  General  Relativity,  P.  G.  Roll,  R.  Krotkov  and  R.  H.  Dicke 
carried  out  a  new  experiment,  using  modern  technology,  and  achieved  an  accuracy 
of  1/100  000  000  000  [8].  Dicke  and  co-workers  were  able  to  increase  the  sensitivity 
compared  to  Eotvos  in  part  by  measuring  the  accelerations  of  their  test  masses  to 
the  Sun,  rather  than  to  the  Earth  as  Eotvos  had  done.  In  such  an  experiment 
any  signal  arising  from  the  difference  between  gravitational  and  inertial  mass  would 
have  the  same  24-hour  periodicity  as  the  Earth’s  rotation.  The  advantage  of  such 
an  approach  from  an  experimental  point  of  view  is  that  it  allows  such  a  signal  to  be 
discriminated  from  background  perturbations,  without  disturbing  the  torsion  fibre. 
Of  course,  one  must  be  careful  to  exclude  other  perturbations  which  will  have  the 
same  24-hour  period.  In  fact  Eotvos,  Pekar  and  Fekete  were  the  first  to  compare  the 
accelerations  of  different  materials  to  the  Sun,  and  for  platinum  versus  mangalium 
they  quote  a  fractional  difference  of  6  x  10~9.  However,  since  no  error  is  quoted  and 
few  other  details  of  their  analysis  are  presented,  it  is  difficult  to  know  precisely  how 
the  sensitivity  of  this  part  of  their  experiment  compared  with  that  of  their  more 
extensive  work  measuring  accelerations  to  the  Earth. 

In  analyzing  the  Eotvos  results,  Dicke  expressed  his  polite  doubts  about  the 
accuracy  claimed  by  Eotvos’  assistants.  Among  his  concerns  were  the  following: 

1.  Dicke  was  suspicious  about  the  perturbation  of  air  motion  created  by  tem¬ 
perature  differences.  The  present  authors  think  that  Eotvos’  team  was  quite 
careful  in  this  respect.  The  observations  were  performed  in  a  shaded  closed 
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Fig.  2.  Sketch  of  the  toreion  balance 


room,  in  a  double-walled  tent  (the  space  between  the  two  walls  was  stuffed 
with  sea-weed),  with  a  torsion  balance  protected  by  triple  copper  coating  with 
air  space  between  (Figs  1,2  and  the  letter  of  J.  Renner  to  R.  H.  Dicke,  due 
to  the  encouragement  of  one  of  us  (G.M.)  as  reprinted  in  Appendix  III.)  Fur¬ 
ther  evidence  for  the  concern  of  Eotvos  and  co-workers  over  thermal  effects  is 
reflected  in  the  fact  that  they  affixed  thermometers  to  the  torsion  balance  at 
various  locations.  (The  double-arm  balance,  used  in  their  third  method,  had 
three  thermometers:  One  along  each  arm,  and  one  near  the  torsion  fibre.) 
The  question  of  possible  thermal  effects  was  raised  by  Dicke  as  an  alterna¬ 
tive  to  the  Fifth  Force  to  explain  the  correlations  in  the  Eotvos  data  noted 
by  Fischbach  and  co-workers.  This  question  has  been  discussed  in  detail  in 
[12],  where  it  is  noted  that  the  main  objection  to  such  an  explanation  of  the 
Eotvos  results  is  that  the  thermal  effects  would  have  to  be  of  constant  sign 
and  magnitude  during  approximately  4000  hours  of  data  taking,  spread  out 
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over  several  yews,  which  seems  unlikely. 

2.  Dicke  was  also  concerned  about  the  poorly  defined  gravitational  perturbation, 
caused  by  the  observer  himself.  This  potential  source  of  error  was  welt  known 
to  Eotvos.  While  the  vibrations  of  the  torsion  pendulum  were  damped,  the 
observer  was  far  away.  When  the  balance  had  come  to  rest,  the  observer  came 
running  and  made  the  reading,  before  the  pendulum  (of  period  of  40  minutes) 
had  time  to  swing  out  [9]. 
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R|.  3.  Dal*  sheet  from  the  R.  Eolvoa  bequest 

3.  Unfortunately,  the  statistical  evaluation  of  the  empirical  data  cannot  be  fully 
reconstructed  in  the  case  of  the  Eotvos- Pekar-Fekete  experiment,  as  the  de¬ 
tails  are  not  given  in  any  of  their  publications.  (In  the  bequest  of  Eotvos,  kept 
in  the  Library  of  the  Hungarian  Academy  of  Sciences,  sheets  with  laboratory 
data  readings  can  be  found,  e.g.  Fig.  3,  but  they  do  not  contain  sufficient 
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information  to  allow  the  statistical  evaluation  to  be  repeated.)  In  the  Hun¬ 
garian  publication  of  J.  Renner  [7]  more  details  can  be  found.  It  can  be  seen, 
as  pointed  out  by  Dicke  [8],  that  Renner  tried  to  eliminate  the  influence  of 
environments  changes  by  interpolation  in  time.  The  numbers  obtained  by 
interpolation  are  not  statistically  independent,  but  Renner  treated  them  as  if 
they  were.  For  this  reason  the  statistical  errors  might  be  larger  by  a  factor  of 
3  than  claimed  by  Renner.  It  should  be  noted  that  this  factor  is  2.4  as  shown 
by  Kiraly  [19].  Taking  this  source  of  error  into  account,  Renner’s  data  lay  sus¬ 
piciously  near  to  zero,  concerning  the  difference  of  the  two  masses.  Therefore 
Dicke  and  others  have  suggested  that  Renner’s  conclusions  cannot  be  relied 
upon.  Renner  learned  this  statistical  technique  from  Eotvbs’  team.  The  inac¬ 
curacy  quoted  by  Pekar  in  his  paper  [5]  originated  from  observational  errors 
and  statistical  errors  in  a  ratio  unknown  to  us;  it  could  be  that  the  accuracy 
of  10~*  [2],  given  by  Eotvos  personally  is  the  reliable  one.  This  is  in  any  case 
a  marvellous  achievement,  and  the  curious  trends,  noted  by  Fischbach  and 
co-workers  seem  to  survive.  (See  the  next  Section  ) 

It  should  be  emphasized  that  the  analyses  of  both  Dicke  and  Fischbach  agree 
that  the  errors  quoted  by  Eotvos,  Pekar  and  Fekete  are  consistent  with  the  statisti¬ 
cal  scatter  of  their  data.  Moreover,  the  confidence  level  of  the  best  fit  of  Fischbach  et 
al  to  the  Eotvos  data,  viewed  as  suggesting  a  Fifth  Force,  is  86%,  which  is  perfectly 
reasonable. 

Let  us  quote  Eotvos  himself  [5]:  “ Ars  longa,  vita  brevis.  The  admonition  of 
this  old  saying  motivates  the  authors  of  this  paper  to  compile  the  results  of  their 
investigation  and  to  submit  them  to  the  judgement  of  a  high  scientific  Aeropag. 
Methods  of  observation  refine  and  improve  naturally  in  the  course  of  observation, 
and  hence  no  mortal  could  close  his  work  if  without  cease  would  follow  the  otherwise 
laudable  impetus  to  replace  the  useful  by  the  even  better.” 


The  hypothesis  of  the  Fifth  Force 

Eotvos’  experiment  is  one  of  the  last  pearls  of  the  grand  epoch  of  classical 
physics.  At  the  end  of  their  investigations  [5]  Eotvos,  Pekar,  and  Fekete  studied  how 
far  the  proportionality  of  inertia  and  gravity  is  valid  in  case  of  radioactive  materials. 
(This  was  already  in  the  era  of  E  =  me3.)  The  proportionality  was  verified  for  a 
0.1  g  sample  of  RaBr2  with  an  accuracy  of  1/2  000  000. 

In  the  following  decades,  the  investigation  of  the  structure  of  matter  called 
attention  to  other  possible  forces  of  Nature  beyond  the  long-ranged  gravity  and 
electricity,  and  beyond  the  short-ranged  nuclear  and  weak  interactions.  According 
to  the  quantum  theory  the  range  of  the  force  (r)  is  related  to  the  mass  (m)  of  the 
quanta  of  the  transmitting  field  by  the  quantum  law  r  =  h/mc,  where  h  is  Planck’s 
constant  and  c  is  the  speed  of  light. 

The  infinite  length  of  the  gravitational  and  electric  field  lines  is  logically  con¬ 
nected  to  the  absolute  conservation  laws  of  mass  (energy)  and  charge.  If  there  is 
any  further  exact  or  approximate  conservation  law  (e.g.  the  conservation  of  the 
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baryonic  charge  B,  discovered  by  E.  P.  Wigner  [10]),  it  may  be  that  a  further  un¬ 
known  field  exists  which  may  transmit  a  Fifth  Force.  But  if  the  rest  mass  of  the 
field  quanta  were  exactly  zero  (as  in  the  case  of  photons),  hot  bodies  would  radiate 
these  quanta  as  well,  in  contradiction  to  thermodynamical  experience.  Thus  one 
may  hypothesize  that  the  field  quanta  should  have  a  small  but  nonvanishing  rest 
mass  (m),  consequently  the  transmitted  Fifth  Force  would  have  a  long  but  finite 
range  (xo  =  h/mc).  This  gives  rise  to  a  composition-dependent  “action  at  a  dis¬ 
tance”  between  two  massive  bodies,  where  the  interaction  energy  is  the  sum  of  the 
gravity  and  the  Fifth  Force  contributions: 

V(x)  =  -Gmm'/x  +  FBB'  exp(-x/zo)/x. 

For  astronomical  distances  only  Newtonian  gravity  contributes, 

V(x)  ~  —Gmm'/x  if  x  >■  xo. 

For  laboratory  distances,  however,  one  may  experience  an  “effective  gravity” 

V(x)  =  —Gejjmm'/x  if  x  C  x0, 
with  an  effective  gravitational  constant 

G?«//  —  G  —  F(B/m)(B' /m'), 

which  may  differ  from  the  (astronomical)  Newtonian  gravitational  constant  G.  If 
B  is  the  baryon  number  (protons  plus  neutrons)  in  the  atom,  and  if  M  is  the  mass 
of  the  atom  in  Hydrogen  atom  mass  units  mu,  that  is  M  =  m/mj*,  then 

GtJJ  =  G[1  -  a(B/M)(B'/M')} 

with  a  —  F/Gmjj.  In  the  case  of  hydrogen  B/M  1,  whereas  for  carbon  B/M  = 
1.00782,  for  copper  B/M  =  1.00895,  and  for  lead  B/M  =  1.00794.  Hence  the 
effective  gravitational  constant,  manifesting  itself  over  laboratory  scales  may  be 
composition  dependent.  This  idea  can  be  checked  by  comparing  the  empirical  value 
of  G  on  astronomical  and  laboratory  scale,  and  by  testing  its  composition-depen¬ 
dence.  Eotvos’  experiments  contributed  to  both.  By  plotting  (see  in  the  work  of 
P.  Kiraly  [19])  the  measured  ratio  inertial  mass/gravitating  mass  with  respect  to 
B/M  (Fig.  4),  from  the  results  published  by  Eotvos-Pekar-Fekete  [5],  Fischbach 
and  co-workers  concluded  that  the  slope  of  the  resulting  line  was  (5.65±0.7)  x  10~6, 
which  differs  from  the  expected  value  of  zero  by  several  standard  deviations  [11].  If 
the  Fifth  Force  really  exists  with  a  range  of  100  m,  say,  the  composition-dependence 
must  be  due  to  the  action  of  nearby  mass  distribution. 

The  Fifth  Force  hypothesis  has  had  a  double  effect:  It  has  encouraged  a  series 
of  modern  experiments,  and  it  has  increased  the  interest  in  details  of  the  original 
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Ebtvos-Pekar-Fekete  experiment  (5]  (which  seemed  to  indicate  a  positive  effect, 
full  lines  in  Fig.  4;  even  at  increased  statistical  error,  indicated  by  dashed  lines 
in  the  Figure),  and  in  the  environment  of  the  Renner  experiment  [7]  (apparently 
suggesting  a  zero  result). 


Fig.  4.  Plot  of  £x{a/g)  u  a  function  of  A(B/M)  (19).  In  EPF  measurement  (thick  lines): 
1.  tallow  -  Cu;  2.  water  -  Cu;  3.  CuS04  solution  -  Cu:  4.  CuSO  crystals  -  Cu;  5.  Asbestos  -  Cu; 
6.  snakewood  -Pt;  7.  AgS04  and  FeS04  (before  and  after  the  reaction);  8.  mangalium  -  Pt; 
9.  Cu  -  Pt.  In  Renner’s  measurement  (thin  lines):  a)  paraffin  -  brass;  b)  NH4F  -  Cu;  c)  Bi  - 
brass;  d)  Pt  -  brass;  e)  glass  -  brass;  f )  Mn  Cu  alloy  -  Cu.  (The  dotted  lines  include  the  increased 

statistical  scatter,  as  indicated  in  the  text.) 

The  authors  of  the  present  report  do  not  intend  to  enter  the  field  of  contro¬ 
versies  related  to  the  Fifth  Force,  the  interested  readers  may  find  references  in  the 
review  papers  [e.g.  12].  Our  only  goal  is  to  supply  information  about  the  environ¬ 
ments  of  the  two  classic  Hungarian  experiments.  But  when  doing  so,  let  us  keep  in 
mind  what  Nieto,  Goldman  and  Hughes  wrote  [13]  •  "neither  the  concept  of  baryon 
number,  nor  the  mass  defect  existed  at  that  time.  Without  these  concepts,  Eotvos 
could  have  spent  considerable  time  and  effort  in  a  fruitless  attempt  to  find  out  why 
the  scatter  in  his  data  points  was  larger  than  his  error  estimates.  We  can  easily 
sympathize  and  imagine  the  gnawing  feeling  that  something  was  wrong,  or  that 
something  very  important  was  being  missed.” 


Eotvos'  Laboratory  revisited 

The  Faculty  of  Science  of  the  University  of  Budapest  (carrying  the  name  of 
Roland  von  Eotvos  since  1950)  is  located  in  the  downtown  of  Budapest,  East  of  the 
Danube.  The  river  follows  a  geological  break:  Its  West  shore  abounds  in  steep  hills 
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Fig.  5.  Geological  cross-section  between  Gelfcrt  Hill  and  Eotvos  University  in  Budapest 
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(CaCC>3,  MgCOa),  its  East  shore  is  flat  (mostly  wet  sand  deposit  of  the  river).  The 
East-West  asymmetry  is  the  dominating  geological  feature  (Fig.  5). 

Eotvos  designed  and  built  the  building  of  the  Physics  Institute  in  the  1880’s, 
his  laboratory  is  now  the  Department  of  Atomic  Physics  (Puskin  utca  5).  Accord¬ 
ing  to  the  Eotvos-Pekar-Fekete  paper  [6,  page  328]  the  laboratory  room  where  the 
Eotvos-Pekar-Fekete  experiment  was  carried  out  looks  South  with  two  windows  on 
the  ground  floor;  opposite  to  it  there  are  tall  buildings  [6,  page  328].  (Recollections 
made  two  decades  later  [13]  contradict  this  original  paper  [5]  but  are  compatible 
with  the  site  of  Renner’s  experiment,  therefore  this  hint  should  be  probably  disre¬ 
garded).  J.  Barnothy  joined  the  Institute  5  years  after  the  departure  of  Eotvos,  and 
he  firmly  locates  the  site  of  the  Eotvos-Pekar-Fekete  experiment  is  a  small  annex 
at  the  SW  end  of  the  building  [14]  (E  on  Figs  6  and  7),  which  now  houses  neutron 
generators.  At  Eotvos’  time  there  was  no  building  to  the  West.  To  the  SW  there 
was  a  temporary  hole  that  was  dug  for  future  construction,  to  the  East  there  is 
the  huge  complex  of  the  Physics  Institute  with  a  strong  concrete  tower,  about  20 
meters  NE  (Fig.  8).  Below  the  experimental  room  there  was  no  cellar  but  only  soil, 
above  it  there  was  no  floor. 

In  contrast  to  the  highly  asymmetric  location  of  the  Eotvos-Pekar-Fekete 
site,  the  Renner  experiment  was  performed  25  years  later,  probably  in  the  geophys¬ 
ical  laboratory,  on  the  North  side,  in  the  middle  of  the  ground  floor  of  the  Institute 
of  Physics  building  (about  where  the  computer  room  is  now).  This  location  of  the 
geophysical  laboratory  is  given  G.  Barta  [15]  indicated  by  R  in  Fig.  7.  Below,  this 
room  there  was  a  cellar,  above  it  one  additional  floor,  which  means  that  the  Renner 
site  is  located  rather  symmetrically  (up-down,  E-W)  in  the  building.  According  to 
Talmadge  et  al  [16],  the  asymmetric  location  of  the  Eotvos  site  may  be  the  source 
of  a  Fifth  Force,  explaining  the  positive  (composition-dependent)  outcome  (Fig.  8). 
The  symmetric  position  of  the  Renner  site  R  (compensated  Fifth  Force)  may  ex¬ 
plain  the  zero  (composition-independent)  outcome  (Fig.  8).  The  explanation  works 
if  the  Fifth  Force  exists  with  a  range  of  10-50  meters.  These  conclusions  have  to  be 
checked  by  modern  experiments  by  observers  who  are  ready  to  learn  patience  from 
Baron  Roland  von  Eotvos.  (G.  Barta,  who  is  presently  repeating  the  original  Eotvos 
Experiment,  speaks  about  2  days  of  waiting  before  one  single  reading  of  the  equilib¬ 
rium  position  of  the  torsion  balance  [17]).  Eotvos  selected  the  best  (most  linear  and 
most  sensitive)  Tungsten  wires  which  hung  with  weight  for  several  years,  to  get  rid 
of  any  distorting  internal  tension.  Eotvos  demanded  thousands  of  hours  of  patient 
unbiased  observations  from  his  assistants.  Let  us  conclude  this  report  from  the  past 
with  Eotvos’  message  [5]:  “The  authors  bow  to  the  fate  of  human  limitations  and 
leave  it  to  future  times  and  future  workers  to  further  elaborate  those  observations 
which  they  themselves  believe  upon  mature  experience  to  be  able  to  still  improve." 


Impact  of  the  Eotvos  experiment 

The  recent  revival  of  interest  in  the  possibility  of  non-Newtonian  gravity, 
which  followed  the  reanalysis  of  the  Eotvos  experiment  by  Fischbach  and  co-work- 
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Fig.  7.  Ground  plan  of  the  building 


ers  [11,12]  owes  much  to  the  perception  that  the  experiments  of  Eotvos,  Pekar,  and 
Fekete  were  carefully  done  and  hence  deserve  to  be  taken  seriously.  The  widespread 
favorable  view  of  this  series  of  experiments  is  due  in  part  to  the  detailed  description 
of  their  experiment  contained  in  the  published  literature,  and  in  part  to  other  details 
of  their  experiment  which  we  know  of  from  personal  communications  [9,14,15,17], 
and  from  aspects  of  their  methodology  that  we  can  infer  from  various  sources.  The 
following  are  two  additional  examples  of  some  of  the  details  of  their  experiment 
which  were  not  described  in  their  paper.  The  torsion  balances  used  in  the  experi¬ 
ment  were  mounted  on  stone  piers  (approximately  one  meter  on  a  side)  which  were 
sunk  deep  into  the  ground.  The  purpose  of  these  piers  was  to  provide  a  stable 
shock-free  platform  for  the  sensitive  balances,  and  a  number  of  these  are  still  visible 
today  at  the  Atomic  Physics  Institute.  Another  interesting  example  deals  with  their 
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Fig.  8.  Campus 


comparison  of  the  reactants  before  and  after  the  chemical  reaction 
Ag2SO«  +  2FeS04  — ♦  2Ag  +  Fe2(S04)3. 

Since  the  Ag  produced  by  this  reaction  precipitates  out  of  the  liquid,  the 
center-of-mass  of  the  initial  reactants  would  not  coincide  with  that  of  the  final 
products.  If  the  difference  of  the  centers  of  mass  were  not  corrected  for,  then  it 
would  couple  to  local  gravity  gradients  and  produce  a  large  (but  spurious)  signal 
which  could  simulate  a  violation  of  Equivalence  Principle  [12].  In  fact  Eotvos  and  co¬ 
workers  found  that  the  accelerations  of  the  reactants  before  and  after  the  chemical 
reaction  were  the  same,  which  is  what  we  expect  in  all  theories.  This  indicates  that 
the  authors  were  evidently  careful  to  correct  for  this  effect,  although  the  details  of 
their  methodology  are  not  provided. 

It  has  now  been  approximately  five  years  since  the  classic  work  of  Eotvos, 
Pekar  and  Fekete  stimulated  interest  in  the  possibility  of  a  fifth  force.  During  this 
period  numerous  experiments  have  been  carried  out,  and  many  are  still  under  way. 
To  date  these  experiments  have  not  confirmed  the  original  suggestion  of  a  fifth 
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force,  as  inferred  from  the  Eotvos  data  by  Fischbach  and  co-workers  [12].  However, 
neither  has  any  g.  oup  pinpointed  an  error  in  the  Eotvos  experiment  which  could  be 
the  source  of  their  suggestive  data.  Since  all  of  the  recent  experiments  differ  from 
the  original  Eotvos  experiment  in  various  ways,  the  possibility  remains  that  there  is 
some  theoretical  model  in  which  a  subtle  aspect  of  the  original  experiment  which  we 
have  heretofore  overlooked  could  explain  why  those  authors  saw  an  effect  while  the 
more  recent  ones  do  not.  The  significance  of  the  Eotvos  experiment  is  that  it  will 
continue  to  be  a  stimulus  for  new  ideas,  such  as  the  recent  suggestion  [18]  that  spin 
may  have  played  a  role  in  the  original  work.  However  the  search  for  new  gravity-like 
forces  turns  out,  it  is  clear  that  the  Eotvos  experiment  has  played  a  fundamental 
role  in  shaping  our  understanding  of  gravity  and  other  possible  forces  in  Nature. 
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Appendix  I 

Nackrickten 

von  der 

Koniglichen  Gesellschaft  der  Wissenschaften 

zu  Gottingen. 


Gesch&ftliche  Uittcilungen. 

1909.  Heft  1. 

Benekeache  Preiastiftung. 

Anf  die  im  Jalure  1906  gestellte  Preisaofga.be : 

:,Von  Etvbs  wurde  eine  «ehr  empfindliche  Hethode  ange- 
geben,  Trigheit*und  Gravitit  der  Materie  an  vergleichen.  Hit 
Riicksicht  hieranf  usd  im  Hinblick  anf  die  neuere  Entwicklang 
der  Elektrodynamik  aowie  auf  die  Entdeckung  der  radioaktiven 
Snbstanzen  ist  daa  Kewtonsche  Gesetz  der  Proportionalitit  von 
Tragheit  and  Gravitit  moglichst  weitgehend  zu  prfifen" 
ist  eine  Bewerbungtschrift  mit  dem  Titel: 

jjBritr’-'ge  zum  Gesetz  der  Proportionalitit  von  Trigbeit  und 
Gravitit'1, 
und  dem  Uotto: 

nArt  longs,  vita  brevis" 

eingelaufen. 

Urn  zu  einem  Urteil  der  Bewerbungtschrift  zu  gelangen,  scheint 
es  wichtig,  die  Gesiebtzpunkte  zu  beachten,  welche  die  Gesellschaft 
bei  der  Stellung  der  Preisfrage  leiteten. 

Das  von  Newton  entdeckte  and  nacb  ihm  benannte  Gesetz 
der  Ailgemeinen  Gravitation,  welches  die  Erscheinongen  der  ma* 
teriellen  Welt  in  ihrer  Gesamtheit  nmfaSt,  spricbt  einige  Sitze 
ana,  die  Bberans  merkwBrdig  sind,  die  man  aber  trotzdem  nicht 
bervorzubeben  pflegt:  1)  Die  Anziehnng  wird  garnicbt  beeinflufit 
durcb  die  pbysikaliscbe  Bescbaffenbeit  der  Materie,  sondem  wird 
einzig  and  aliens  darcb  die  „TrSgheit"  bestimmt.  Mit  dieser  Trig* 
heit  sind  die  Fernwirkangen  proportional,  sodafi  man  knrz  den 
Satz  formnlieren  kann:  Das  Verh&ltnis  von  Gravitit  and  Tragheit 
ist  fUr  alle  materiellen  Teile  anverinderlich  and  fUr  alle  gleich 
grofi.  Da  alien  materiellen  Teilen  eine  nnverinderliche  Tragheit 
anznhaften  scheint,  so  wire  zu  folgern,  dad  aucb  die  Gravitit 
eine  nnverinderliche  Eigenschaft  der  Materie  ist.  2)  Die  Fern* 
wirkung  irgend  zweier  materieller  Teile  wird  dorch  die  Anwesen- 
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33  Denekesche  rrcisitiftung 

beit  der  ubrigen  materiellen  Teile  nicbt  beeinfluflt  —  Ein  Teilchen 
im  Innern  der  Erde  und  eines  inmitten  des  Sonnenb&lles,  ziehen 
biernach  einander  gerade  so  an,  als  ob  die  Ton  dem  Erdkorper 
and  dem  Sonnenkorper  gebotenen  materiellen' Mantel  garnicht  vor- 
handen  waren.  3)  Die  Fernwirkung  hangt  allein  von  der  jewei- 
ligen  gegenseitigen  Lage  der  materiellen  Korper,  nicbt  von  ihrem 
Bewegangszustand  ab.  Es  scbeint  biernach  die  Gravitation  sich 
in  nnendlicher  Geschwindigkeit  aoszabreiten.  — 

Man  hat  sicb  an  den  Gedanken  der  unbeschrankten  Giiltigkeit 
des  Xewtonscben  Gesetzes  so  sdbr  gewobnt,  dab  das  Gefiihl  fUr 
die  Merkwiirdigkeit  der  hervorgehobenen  Satze  fast  verloren  ge- 
gangen  ist.  — 

In  ein  neues  Licbt  ist  die  Frage  nacb  der  Giiltigkeit  der 
Xewtonschen  Gesetze  geriickt  worden  durch  die  Erfolge  der  tbeo- 
retiscben  Elektrodynamik.  Auch  in  den  elektrischen  and  magne- 
tischen  Kraften  hatte  man  anfanglicb  Einwirknngen  ganz  ahnlicher 
Art  wie  bei  der  Gravitation  gesehen.  Nun  lehrte  Maxwell,  dafl 
die  elektriscben  und  magnetiscben  Krafte  sicb  nicbt  mit  nnend¬ 
licher,  sondern  mit  der  Geschwindigkeit  des  Lichtes  ansbreiten. 
Er  zeigte  ferner,  daB  die  Wechselwirkang  bei  den  magnetischen 
and  elektriscben  Erscheinangen  sebr  wesentlich  je  nacb  der  Art 
des  Zwischenmediums  variiert.  Dadarcb  schon  warde  die  Pbysik 
von  neaem  angeregt,  die  Giiltigkeit  der  Satze  2)  und  3)  in  Zweifel 
zu  stellen.  Xocb  tiefere  theoretische  Bedeutang  gewann  der  Satz 
1),  der  die  ProportionalitSt  von  Gravitat  und  Tragbeit  ausspricht. 
Es  gelang  der  Elektrodynamik  za  zeigen,  daB  mit  einer  elektri¬ 
schen  Ladung  das  Besteben  einer  Tragheit  im  Sinne  der  Mechanik 
verbnnden  ist,  and  es  warde  fcstgestellt,  daB  im  Innern  eines  jeden 
materiellen  Korpers  eine  auBerordentlich  groBe  Zab  1  sehr  stark 
elektrisch  geladener  kleiner  Teilchen  vorhanden  ist.  Danach  nan 
schien  es  nicht  aasgeschlossen,  daB  ein  wesentlicber  Teil  der  be- 
obacbteten  Tragheit  der  Materie,  vielleicbt  die  Tragbeit  iiberhanpt, 
sich  elektrodynamisch  erklfire.  So  warde  die  Gravitation,  die  doch 
mit  der  Tragheit  znsammenbSngt,  jetzt  aacb  in  Verbindang  mit 
der  Elektrodynamik  gebracht,  and  es  riickte  der  Satz  von  der 
Proportionality  von  Gravitat  and  Tragbeit  in  eine  nene  fiber- 
raschende  Beleuchtnng. 

Diesen  Erwagangen  weiter  nachgehend,  welcbe  za  der  Frage 
ffihren,  wie  die  Materie  in  das  physikaliscbe  Weltbild  einzaordnen 
ist,  schien  es  der  Fakaltat  besonders  wicbtig,  dem  Satz  von  der 
Proportionality  der  Tragbeit  and  der  Gravitat  erneat  die  Aaf- 
merksamkeit  za  schenken,  and  insbesondere  erscbien  die  denkbar 


Acta  Phfiica  Hang  aria  69,  1991 


102 


EOTVOS  EXPERIMENT 


351 


Denekeiche  Preiutiftong.  39 

scharfste  experimentclle  Priifur.^  dr  intend  erwiiascht.  Dies  war 
der  AnlaB  zur  Stellung  der  Preisaafgabe  fUr  1909.  — 

Die  Bewerbungsschrift  mit  dem  Motto:  ,,Ar«  longa,  vita 
brevis1'  geht  aaf  die  theoretischen  Erwagangen,  an  die  erinnert 
wnrde,  garnicht  ein,  beruhrt  sie  nicbt  einmal.  Damit  ist  klar, 
da£  einem  wesentlicben  Toil  der  Wfinsche,  welcbe  die  Fakoltft 
bei  Stellung  der  Preisaufgabe  leiteten,  nicbt  Recbnong  getragen 
wird.  Dafiir  wird  die  ganze  Kraft  aaf  die  Ausffihrang  der  ex- 
perimentellen  Untersacbang  verwendet,  and  es  wird  aacb  gezeigt, 
wie  man  nnter  Verwertang  bekannter  Erfabrangen  fiber  die  Er- 
scbeinaag  bei  Ebbe  and  Flat  weitere  wertvolle  Folgerangen  fiir 
das  bier  zur  Bebandlang  stebende  Problem  zieben  kann.  — 

Es  werden  obne  eine  Aenderung  die  darch  Etvos  konstraierten 
Apparate  benutzt.  Dem  Studiam  der  Fehlerqaellen  wird  eine 
grofie  Aufmerksamkeit  geschenkt,  sodaB  die  Beobachtungen  einen 
bohen  Grad  der  Zoverliissigkeit  erhalten.  Die  gewonnenen  Resul- 
tate  sind  so  wertvoll,  daB  es  die  Fakoltat  mit  Genagtnnng  begruflen 
darf,  darch  Stellang  der  Preisaafgabe,  za  den  Beobacbtangen  An- 
laB  gegeben  zu  haben.  — 

Newton  fand  fur  eine  Reihe  von  antersuchten  Materialien. 
daB  Gravitat  and  Tragheit  jedenfalls  bis  auf  etwa  1/1000  ihrer 
GrBBe  mit  ein&nder  proportional  sind.  Bessel  zeigte,  bei  Ge- 
legenbeit  seiner  Pendelversucbe,  daB  die  etwaigen  Abweicbongen 
hocbstens  1:60000  erreicben  konnten.  EtvSs  hat  mitgeteilt,  daB 
es  ihm  mit  teinen  Apparaten  mdglicb  geworden  sei,  die  Propor- 
tionalitfit  bis  aaf  1/20000000  za  erweisen;  er  macht  aber  keine 
Angaben  fiber  das  Beobacbtangsmaterial.  In  der  Preisschrift  wird 
aaf  Grand  des  neaen  Beobachtaagsmaterials  nachgewiesen,  daB 
fiir  eine  ganze  Reihe  sebr  verschiedener  Materialien  (Platin,  Ma- 
gnalium,  Kupfer,  Wasser,  krystallisiertes  Kapfersolfat,  Kapfer- 
salfatlosang ,  Asbest  and  Talg)  die  Abweicbong  von  dem 
Gesetz  der  Proportionalitat  jedenfalls  nicbt  groBer  als  etwa 
1/200000000  ist.  Aacb  lieB  sicb  zeigen,  daB  mit  der  chemiscben 
Reaktion  Silbersulfat  -  Ferrosolfat,  die  seinerzeit  von  L  an  dolt 
verwertet  wnrde,  am  den  Satz  von  der  Konstanz  der  Masse  za 
priifen  and  mit  der  Aaflosang  von  Kapfersalfat  in  Wasser,  die 
Heydweiller  za  gleichem  Zwecke  verwendete,  Yeranderangen 
des  Yerbaltnisses  von  Gravitat  and  Tragbeit  in  dem  angegebenen 
Grenzbetrage  jedenfalls  nicbt  eintreten.  —  Yon  besonderem  In- 
teresse  ist  aacb,  daB  eine  sebr  stark  radioaktive  Sabstanz,  Radiam- 
bromid,  der  Untersacbang  aaterworfen  warde.  Die  experimentellen 
Schwierigkeiten  waren  bier  natargemaB  sehr  viel  groBer  nnd  d&ram 
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die  Scharfe  lies  Rcsnltates  erheblich  gcringer.  E»  ergab  sich,  dad 
cine  etwaige  Abweichung  des  Verbal  tnisses  von  Gravitiit  and 
Trigbeit  bei  dem  Praparat  jedenfalls  nicht  griSBer  war  als  etwa 
1/2000000.  — 

Auch  der  oben  unter  2)  anfgefubrte  Satx  wird  einer  Prufung 
nnterzogen,  nach  dem  die  Wecbaelwirknng  der  Gravitation  dorch 
Zwiacbenscbalten  von  Uaterie  nieht  beeinflofit  warden  toll.  Dabei 
warden  Vertuche  angefUhrt,  die  schon  vor  langeren  Jabren  (1902) 
angastellt  worden  tind,  denen  aber  nor  der  Cbarakter  alt  Vor- 
vertncb  beigelegt  wird.  Ei  11  fit  sich  folgern,  dafi  durch  eine 
zwischenliegende  Bleiplatte  von  einer  Dicke  gleich  dem  Dorch- 
metter  der  Erde  die  Gravitation  nm  nicht  mehr  ah  etwa  am  1/800 
ihres  Betrages  geandert  wiirde.  Hieran  kniipfen  die  Verfaster 
theoretische  Erorterangen  fiber  die  Ebbe-  and  Fluterscheinongen 
and  folgern,  —  insbetondere  auch  aus  den  Heckerscben  Beob- 
achtungen  der  Flaterzeagenden  Kraft  mittels  des  Horizontalpen- 
dels  —  dafi  die  Zwischenschaltung  der  ganzen  Erde  die  Anziebong 
der  Sonne  anf  ein  materielles  Teilcben  am  weniger  ah  den  10000. 
Teil  andert.  — 

Das  Endresoltat  der  ganzen  Arbeit  wird  so  ausgcsprocben : 
„Wir  babeq  eine  Reihe  von  Beobachtungen  angestellt,  die  an  Ge- 
naaigkeit  alle  vorangebenden  fibertrafen,  doch  konnten  wir  in 
keinem  Falle  eine  bemerkbare  Abweichung  von  dem  Gesetz  der 
Proportionalitat  von  Tragheit  and  Gravitiit  entdecken*. 

Die  Verfaster  bemerken  zu  Anfang  ihres  Berichtes :  Blit  Riick- 
sicht  anf  die  Kiirze  der  Zeit,  die  ons  fiir  die  genauere  Durchsicbt 
onserer  Arbeit  znr  VerfUgung  stand,  bitten  wir  fiir  eventuell  vor- 
kommende  Scbreibfehler  and  das  W'escn  der  Resultate  nicht  be- 
eintrachtigende  Recbenfehler  am  Kacbsicbt.  So  mag  denn  nicht 
viel  Gewicht  daranf  gelegt  werden,  dafi  in  der  Tat  bei  der  Be- 
arteilong  der  Flatwirkang  ibre  Darlegangen  mebrfacb  Verbesse- 
rongen  and  Vervolhtandigangen  bediirfen.  — 


Es  ist  gewifi,  dafi  die  Verfasser  der  Preisarbeit  in  sehr  we: 
sentlicben  Ponkten  den  Erwartnngen  der  FakoltSt  nicht  entsproeben 
baben,  and  es  mafi  each  bemerkt  werden,  dafi  in  Einzelheiten  die 
AnsfObrongen  nicht  anerkannt  twerden  konnen.  Trotzdem  aber 
bringt  die  Arbeit  hochst  wertvolle  Resnltate,  indem  aie  als  Grand- 
lage  fiir  alle  tbeoretischen  Speknlationen  den  anfierordentlieb 
weitgebenden  GBltigkeitsbereicb  der  N ewtonschen  Gesetze  zeigt. 
Die  FaknlUt  stebt  darnm  nicht  an,  der  Arbeit  den  vollen  Preis 
an  erteilen. 

Gfittingen,  den  1.  April  1909. 

Die  pliilosopbisclie  Fakultat. 

Der  Dekan: 

C.  Range. 
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Dear  Professor  R.  H.  Dicks 1 

I  read  with  great  interest  your  valuable  treatise  about  the  Botvos 
experiment  published  in  the  "Scientific  American"  in  the  number  of  December  1961. 
I  know  your  previous  studies  in  this  matter  too.  1  am  fully  convinced  of  the 
great  importance  of  your  experiments  carried  out  with  up  to  date  methods.  The 
accuracy  achieved  by  your  experiment  is  much  greater  than  the  accuracy  of  the 
previous  works  done  by  R.  Botvos  and  his  collaborators.  The  fundamental  question 
of  the  independence  of  gravitational  acceleration  from  the  quality  of  material 
can  be  considered  now  as  proved.  I  have  to  congratulate  for  your  experimental 
results  of  high  accuracy. 

I  was  student  and  later  on  collaborator  of  Botvos  and  I  repeated  his  famous 
experiments  in  the  years  1930-1935.  I  bag  you  to  allow  me  some  suggestions 
relating  your  remarks  about  the  original  experiments  of  Botvos.  I  have  to 
mention,  that  your  interesting  treatise  was  translated  in  Hungarian  and  published 
in  the  periodical  "Fizikai  Szemle."  My  suggestions  were  published  in  the  same 
periodical.  I  enclose  a  copy  of  it  in  Hungarian. 

As  the  Hungarian  experiments  were  made  by  visual  readings,  the  mass  of  the 
observer  could  have  exert  any  influence  to  the  position  of  equilibrium.  But  the 
observing  time  was  always  rather  short  and  the  observer  did  not  stay  for  long 
time  near  to  the  instrument,  as  you  already  assumed  it. 

Concerning  the  influence  of  sudden  changes  of  temperature,  Botvos  took  care 
to  eliminate  the  unfavourable  effects.  The  experiments  were  carried  out  in  a 
dark  room,  where  temperature  was  practically  steady.  In  the  room  there  was 
mounted  a  linen  tent,  the  walls  of  it  were  filled  with  isolating  material  and  the 
instrument  was  brought  inside  this  tent.  Besides  the  instrument  had  a  triple 
shell  of  metal  around  the  chamber  containing  the  torsion  balance,  as  Professor 
had  already  mentioned  in  the  publication.  Accordingly  there  was  no  reason  of 
taking  place  any  convection  currents  in  the  inside  of  the  balance.  Specifically 
prepared  torsion  wires  were  used,  which  had  no  any  drift  caused  by  elastic 
properties  or  by  change  of  temperature.  The  positions  of  equilibrium  appeared 
even  in  long  observation  sets  very  constant. 

You  have  correctly  suggested  the  disturbing  effect  of  even  quite  small 
magnetic  impurities  in  the  moving  parts  of  the  balance.  Botvos  and 
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his  collaborators  aliminatsd  tha  assumable  effect  by  compensating  the  geomagnetic 
field  with  permanent  magnets  and  electromagnetic  coils.  The  compensation  was 
always  controlled.  The  compensating  magnets  were  so  mounted,  that  they  should 
not  produce  any  translatoric  effect.  Besides  the  materials  used  in  the  balance 
were  controlled  concerning  their  magnetic  properties. 

There  is  another  question,  that  Botvos  and  his  collaborators  used  a 
horizontal  variometer,  in  which  one  of  the  masses  was  suspended  lower,  so  that 
the  horizontal  gradient  of  the  acceleration  had  influence  too.  This  caused  no 
any  disadvantage,  as  the  method  of  Eotvos  counted  on  the  effect  of  the  gradient. 
When  different  materials  were  hanged  on  the  balance,  the  centre  of  gravity  was 
always  in  the  same  height  and  the  balance  took  place  in  the  same  gravitational 
field.  The  optical  diffraction  in  the  telescope  caused  no  erroneous  readings. 
The  20-th  part  of  a  scale  division  could  be  precisely  read  by  a  skilled  observer. 

Concerning  my  own  experiments  carried  out  in  the  years  1930-1935,  I  have  to 
mention,  that  I  did  not  use  the  old  Eotvos  balance,  but  an  improved  new  one,  the 
thermical  and  magnetic  effects  were  perfectly  eliminated  and  very  reliable 
torsion  wires  were  used.  The  positions  of  equilibrium  were  kept  in  long 
observation  sets  constant.  In  this  way  it  succeeded  to  extend  the  accuracy. 

I  would  be  very  grateful,  if  you  would  be  kind  to  publish  a  remark  of  my 
suggestions  above  concerning  the  original  experiments  of  Eotvos  and  his 
collaborators,  if  possible  in  the  same  peeriodical,  in  which  your  interesting 
study  appeared. 


With  kindest  regards 


Budapest,  26.  VII.  1963 


Yours  faithfully 


/  Dr.  J.  Renner/ 

Budapest  VII.  Damjanich  u  28/b 
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Six  years  of  the  fifth  force 

Ephraim  Fischbach  &  Carrick  Talmadge 


The  enunciation  of  the  ‘fifth  force’  hypothesis  in  1986  spawned  a  generation  of  experiments  searching 
•  for  deviations  from  newtonian  gravity.  Although  no  compelling  evidence  for  any  new  weak  forces  has 
emerged  in  the  past  six  years,  the  searches  for  anomalous  gravitational  effects  have  produced  a  large 
number  of  important  experimental  and  theoretical  results. 


The  suggestion  roughly  six  years  ago  of  a  possible  gravity-like 
'fifth'  force1  of  nature  called  attention  to  the  fact  that  gravity, 
itself  the  ‘oldest*  of  the  known  forces,  was  in  some  ways  also 
the  least  well  understood.  As  we  discuss  below,  a  fifth  force 
coexisting  with  conventional  gravity  would  lead  to  a  net  interac¬ 
tion  between  macroscopic  objects  which  showed  small  devi¬ 
ations  from  the  behaviour  expected  from  the  classical  newtonian 
inverse-square  law  of  gravity  (see  equations  (2)  and  (3)  below). 
Evidence  for  such  non-newtonian  behaviour  in  an  appropriate 
system  could  thus  be  a  sign  of  a  new  fundamental  force,  and 
this  idea  has  provided  part  of  the  impetus  for  the  current  efforts2,1 
to  re-examine  the  experimental  support  for  newtonian  gravity 
over  various  distance  scales.  Although  no  compelling  experi¬ 
mental  evidence  for  any  anomalous  results  has  yet  surfaced, 
there  has  nonetheless  been  a  dividend  from  the  intense  effort 
of  the  past  few  years.  This  is  the  substantial  improvement  in 
the  precision  with  which  newtonian  gravity  can  be  claimed  to 
be  experimentally  supported.  Moreover,  the  search  for  non- 
newtonian  effects  continues,  spurred  on  both  by  remarkable 
advances  in  experimental  techniques  and  by  an  improved  under¬ 
standing  of  the  theoretical  basis  for  such  effects.  Specifically, 
we  now  recognize  that  many  models  predict  the  existence  of 
new  intermediate-range  forces,  and  relate  the  properties  of  these 
forces  to  the  behaviour  of  physics  at  the  Planck  (energy)  scale. 

It  is  important  to  bear  in  mind  that  although  we  will  focus 
here  primarily  on  what  has  been  learned  in  this  field  since  the 
fifth  force  was  proposed  in  1986,  this  proposal  was  itself  a 
natural  outgrowth  of  many  earlier  experimental  and  theoretical 
studies  of  non-newtonian  gravity2.  In  the  1970s  interest  was 
stimulated  by  the  work  of  Fujii4  and  others3,6  who  argued  that 
various  theoretical  models  led  naturally  to  a  new  intermediate- 
range  force,  whose  strength  was  comparable  to  (or  perhaps 
slightly  weaker  than)  gravity.  Although  Fujii’s  theory  was  based 
on  detailed  dynamical  considerations,  at  its  core  are  a  few 
fundamental  observations  which  have  helped  to  motivate  much 
of  the  subsequent  work.  We  begin  by  noting  that  there  are  two 
natural  mass  scales  in  physics,  defined  by  mNal  GeV/c2  and 
Mpi,nck  =  (hc/G)1/2ai01,GeV/c2,  where  mN  is  the  nucleon 
mass,  ft  is  Planck's  constant,  c  is  the  speed  of  light,  and  G  = 
(6.67259  ±  0.00085)  x  10'"  N  m2  kg-2  is  the  newtonian  gravita¬ 
tional  constant.  The  ratio  of  these  scales,/ *  mN/  MP,,nck  =  10-19, 
introduces  a  new  small  dimensionless  constant  into  physics. 
.  This  constant  may  have  some  dynamical  significance  if  it  deter¬ 
mines  the  coupling  strength  of  a  new  field  to  matter,  just  as  the 
electric  charge  e  (e2/ftc*  atm  a  1/137)  determines  the  strength 
of  the  coupling  of  photons  to  matter.  If  this  were  the  case,  then 
•  the  analogue  of  a.m  for  this  field  would  be  f2/  ftc  s  10-38,  which 
corresponds  to  a  force  of  gravitational  strength  (p/Gm2^  =  1). 
The  second  observation  is  that  in  various  theories  the  product 
•J(f2/  ftc)  mN»  fi=  1<T10  eV/c2  determines  the  mass  of  the  new 
field.  As  the  Compton  wavelength  (which  characterizes  the 
distance  over  which  a  force  acts)  associated  with  this  field  would 
be  A  =  h/fic  a  2,000  m,  this  combination  of  values  for  the  para¬ 
meters  /  and  fi  could  describe  a  new  intermediate-range  field 
of  gravitational  strength. 

To  understand  how  such  a  field  could  be  detected  experi¬ 
mentally,  we  note  that  the  potential  energy  Vs( r)  describing  the 


interaction  of  two  point  masses  m,  and  m(  through  th 
force  would  be  given  by 


(I) 


where  by  convention  the  sign  has  been  chosen  to  make  V}(r) 
correspond  to  an  attractive  force  for  a>  0.  Here  r  =  |r,-rJ|  is 
the  distance  between  the  masses,  and  a  is  a  dimensionless 
constant  proportional  to  / 2  which  determines  the  strength  of 
the  new  force  (relative  to  gravity).  Because  a  can  be  either 
positive  or  negative  in  various  theoretical  models  (corresponding 
to  a  force  which  is  attractive  or  repulsive),  the  sign  of  a  carries 
important  physical  information.  For  example,  a  repulsive  force 
typically  arises  from  the  exchange  of  a  vector  (spin-1)  meson 
between  nucleons  or  electrons.  Theories  of  such  forces  lead  us 
to  expect  that  the  magnitudes  of  the  effects  they  produce 
necessarily  depend  on  the  chemical  compositions  of  the  test 
masses.  By  contrast,  an  attractive  force  between  like  objects 
arises  from  the  exchange  of  scalar  (spin-0)  and/or  tensor  (spin- 
2)  fields,  and  these  can  lead  to  deviations  from  the  1/r2  law 
that  are  independent  of  the  compositions  of  the  test  masses.  It 
follows  that  an  experimental  determination  of  the  sign  of  a  can 
be  used  to  discriminate  among  different  theories  of  the  fifth 
force.  Returning  to  equation  (1),  we  note  that  the  sum  of  V5(r) 
and  the  usual  newtonian  potential  VN(r)  can  be  written  as 


— Saa<i+«.-«) 


(2) 


In  the  form  of  equation  (2),  the  presence  of  the  new  interaction 
described  by  V5(r)  manifests  itself  as  an  apparent  deviation 
from  the  usual  1/r  newtonian  potential.  In  practice  nearly  all 
experiments  measure  a  force  (or  acceleration)  rather  than  a 
potential,  and  from  equation  (2)  we  find  for  this  force 

F (,)  =  — VV(r)  =  — G(f)!® 

r  (3) 

G(r)  =  G41  +  a(l  /A)e-r/A] 

In  equation  (3)  we  have  written  the  newtonian  constant  as  G„ 
to  indicate  that  this  is  the  constant  that  determines  the  strength 
of  the  interaction  in  the  limit  r-»oo.  This  serves  to  distinguish 
G*  from  the  constant  G0=  G*(l  +a),  which  characterizes  the 
interaction  strength  when  r/A«  1.  (The  latter  regime  is  often 
referred  to  as  the  ‘laboratory  scale’  on  the  assumption  that  A  is 
larger  than  the  typical  separation  of  test  masses  in  laboratory 
experiments.)  It  follows  from  equation  (3)  that  any  deviation 
of  G(r)  from  G „  leads  to  a  breakdown  of  the  usual  1/r2  force 
law,  and  can  be  interpreted  as  evidence  for  a  new  force. 


Elementary  phenomenology 

The  work  of  Fujii  and  others,  which  led  to  the  phenomenological 
framework  given  in  equations  (l)-(3),  motivated  further  experi¬ 
mental  and  theoretical  work,  as  discussed  in  ref.  7.  We  see  from 
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equation  (3)  that  every  experimental  result  implies  some  con¬ 
straint  on  the  parameters  a  and  A  which  deline  the  putative 
non-newtonian  interaction.  For  technical  reasons,  the  limits  on 
a  =  <*(A)  implied  by  a  given  experiment  are  most  stringent  for 
values  of  A  that  are  comparable  to  the  characteristic  separation 
r  of  the  test  masses  in  that  experiment.  It  follows  that  no  single 
experiment  can  provide  useful  limits  on  <*(A)  for  all  values  of 
A.  Rather,  a  collection  of  experiments  is  needed,  each  optimized 
to  measure  a(A)  near  a  selected  value  of  A.  By  1981  a  number 
of  such  results  were  available,  and  these  were  collected  together 
in  an  important  paper  by  Gibbons  and  Whiting8  (hereafter 
referred  to  as  GW). 

The  constraints  on  a  and  A  implied  by  the  data  available  to 
GW  are  shown  by  the  dark  shaded  area  in  Fig.  1.  In  this  and 
similar  figures  (see,  for  example,  refs  9,  10),  the  shaded  region 
in  the  a-A  plane  is  excluded  by  the  data  at  the  two-standard- 
deviation  (2 cr)  level.  The  GW  paper  drew  attention  to  the  region 
(10m«  A«  JO3  m)  where  the  limits  on  a  (and  hence  on  the 
validity  of  newtonian  gravity)  were  very  poor,  as  can  be  seen 
from  Fig.  1.  This  gap  in  our  knowledge  simply  reflects  the 
difficulty  in  designing  experiments  in  which  the  separation  of 
the  test  masses  is  of  order  10-10’ m.  Typically  information  on 
gravity  for  such  mass  separations  comes  from  geophysical 
experiments,  and  hence  this  range  of  distances  is  often  referred 
to  as  the  ‘geophysical  window'. 

Stacey  and  collaborators11'14  undertook  to  explore  the 
‘geophysical  window’  by  reviving  the  Airy  method  of  determin¬ 
ing  the  newtonian  gravitational  constant  G.  A  schematic  outline 
of  the  Airy  method  is  shown  in  Fig.  2  for  an  idealized  model 
of  a  spherical  non-rotating  Earth.  Suppose  we  compare  the 
acceleration  g(0)  of  a  test  mass  at  the  surface  of  the  Earth  with 
its  acceleration  g(z)  at  a  depth  z  below  the  surface.  (The  test 
mass  in  each  case  is  simply  the  proof  mass  in  a  standard 
gravimeter,  usually  of  the  LaCoste-  Romberg  type.)  It  is  straight¬ 
forward  to  see  that  there  are  two  different  effects  that  contribute 
to  Ag(z)«g(z)-g(0),  and  that  these  have  opposite  signs.  The 
first  is  the  free-air  gradient,  which  represents  the  increase  in 
g(z)  arising  from  the  circumstance  that  the  test  mass  at  z  is 
only  a  distance  (R  -  z)  from  the  centre  of  the  inner  mass,  rather 
than  at  R  as  would  be  the  case  at  the  surface  (see  Fig.  2).  The 
second  effect,  known  as  the  double-Bouguer  term,  accounts  for 
the  fact  that  at  a  depth  z,  g(z)  will  decrease  because  the  test 
object  is  now  attracted  to  the  centre  of  the  Earth  by  a  smaller 
total  mass  (M-4M).  (Recall  that  if  Newton’s  1/ r2  law  of 
gravity  holds,  then  the  shell  of  matter  of  thickness  z  does  not 
contribute  to  g(z).)  From  Fig.  2  and  the  preceding  discussion 


we  see  that  in  this  simplified  model  Ag(z)  is  given  by 


Ag(z)  =  g(z)-g(0)  = 


(R-z)2  R2 


^2g(0  U 

R 


-  4  irGpz 


(4) 


where  p  is  the  density  of  the  material  in  the  shell,  and  g(0)  = 
GM/R2.  It  follows  from  equation  (4)  that  G  can  be  determined 
over  geophysical  scales  by  a  series  of  local  measurements  to 
determine  g(z),  g(0)  and  p,  and  this  constitutes  the  Airy  method. 
In  practice  equation  (4)  must  be  modified  to  include  the  effects 
of  the  Earth’s  rotation  and  ellipticity11"15.  The  original  analysis 
by  Stacey  and  coworkers1’  in  the  early  1980s  in  fact  found 
an  anomalously  high  value  of  G:  their  best  value,  obtained 
from  a  mine  at  Mount  Hilton  in  Australia,  was  G  = 
(6.720 ±  0.024)  x  I0-11  N  m2  kg'2  which  was  larger  than  the 
laboratory  value  G  =  (6.67259 ±  0.00085)  x  10~H  N  m2  kg'2  by 
roughly  twice  their  maximum  admitted  error. 

We  note  for  later  purposes  that  as  G(r)  always  appears  in 
combination  with  a  mass  such  as  m,,  G(r)  cannot  be  determined 
separately  unless  m,  is  independently  known.  In  the  Airy  method 
the  appropriate  combination  is  G(r)AM,  and  as  AM  can  be 
determined  from  measurements  of  the  local  density  p,  G(r)  can 
be  inferred  for  values  of  r  comparable  to  the  depth  of  the 
measurement.  But  in  some  situations,  such  as  tower  experiments 
and  tests  involving  planetary  motion  (see  below),  the  appropri¬ 
ate  masses  cannot  be  determined  independently.  In  these  other 
cases  G  cannot  be  obtained  in  absolute  terms,  and  only  the 
variation  of  G(r)  with  r  can  be  measured. 

During  the  same  period  of  the  early  1980s,  Aronson  el  aL 16 
analysed  high-energy  Fermilab  data  on  the  K°-  R°  system  to 
look  for  anomalous  effects  that  might  be  associated  with  the 
existence  of  a  new  force.  They  used  neutral  kaons  because  the 
small  mass  difference  Am  =  mL- ms  =  4x  10~6  eV/c2  of  the  weak 
interaction  eigenstates  KL  and  Ks  makes  this  system  extremely 
sensitive  to  small  external  influences1718.  The  presence  of  a  new 
weak  force  would  show  up  as  an  apparent  dependence  of  Am 
on  the  velocity  v  of  the  kaons  (or  their  energy)  in  the  laboratory, 
because  the  strength  of  an  external  field  as  seen  in  the  rest  frame 
of  the  kaons  depends  on  v  through  the  Lorentz  transformations. 
Aronson  et  al.16  also  noted  that  similar  velocity-dependent 
effects  could  show  up  in  the  other  fundamental  parameters 
of  the  K°-K°  system,  such  as  the  lifetimes  tl  and  rs  of  Kt 
and  Ks,  and  CP- violating  parameter  v*  *  amplitude  (KL-» 
ir*ir~)/ amplitude  {Ks-»  ir*ir~).  Motivated  by  these  theoretical 


FIG.  1  Constraints  on  the  coupling  constant  a  as  a  function 
of  the  range  A  from  composition-independent  experiments. 
The  dark  shaded  area  indicates  the  status  as  of  1981.  and 
the  lighter  region  the  current  limits.  For  references  to  the 
earlier  experiments  which  contribute  to  the  curves,  see  refs 
9.10. 


10'1  -r- 

102  -j 
10-3 
10-4  -J 
10-* 

10-* 

IO7 
10-* 

10> 
10-10 


1981 


1991 


LAGEOS 


LAGEOS  -  Lunar 


N 

Planetary 


103  10-'  10’  103  105  107  10*  10"  10’3  10'5 

Mm) 


110 


NATURE  •  VOL  356  ■  19  MARCH  1992 


REVIEW  ARTICLE 


considerations,  they  reanalysed  data  from  a  series  of  Fermilab 
experiments  which  determined  the  K°-K°  parameters  for  kaons 
in  the  energy  range  30-130GeV.  They  found  evidence  at  a 
marginally  interesting  (s3<r)  level  for  an  energy  variation  of 
the  kaon  parameters,  particularly  for  the  phase  <6,  of  i)±, 
and  this  provided  at  that  time  another  hint  for  a  possible  new 
force.  We  will  return  later  to  discuss  the  present  status  of  kaon 
experiments. 

The  ‘fifth  force'  hypothesis1'’  in  1986  arose  from  the  recog¬ 
nition  that  the  anomalies  reported  by  both  Stacey  ei  aL"  '*  and 
Aronson  et  al.' 14  could  be  explained  in  terms  of  the  same  new 
force,  but  only  if  it  had  properties  that  were  somewhat  different 
from  those  of  V,(r)  in  equation  (2).  Additionally,  the  requisite 
force  would  be  capable  of  producing  anomalies  in  the  classic 
experiment  of  Eotvos,  Pekar  and  Fekete  (EPF)19-20,  which  com¬ 
pared  the  accelerations  of  different  pairs  of  materials  to  the 
Earth.  This  motivated  a  reanalysis 1,7  of  the  Eotvos  experiment, 
from  which  two  surprising  observations  emerged.  First,  the 
results  of  EPF  indicated  that  for  some  pairs  of  materials,  the 
acceleration  differences  for  the  two  samples  were  relatively  large 
(So  in  the  case  of  copper  against  H20).  Second,  these  acceler¬ 
ation  differences  fitted  a  pattern  that  was  compatible  with  the 
existence  of  a  new  ‘fifth  force’.  This  observation,  coupled  with 
the  suggestion  that  the  same  force  might  also  account  for  the 
earlier  anomalies13'16  helped  to  revive  interest  in  non-newtonian 
gravity. 

The  proposed  fifth  force  departed  from  the  earlier  work  of 
Fujii  by  introducing  a  repulsive  interaction  proportional  to  the 
hypercharge  of  a  test  object  V  =  B  +  S  ( fl  is  baryon  number,  S 
is  strangeness),  rather  than  to  its  mass  as  in  equation  (2).  For 
ordinary  bulk  matter  B  =  N  +  Z  and  S  =  0,  where  N  and  Z 
respectively  denote  the  number  of  neutrons  and  protons  in  a 
sample,  and  hence  Y  and  B  can  be  used  interchangeably.  But 
for  certain  elementary  particles  such  as  kaons  (specifically  K° 
and  its  antiparticle  K°),  B  =  0  and  S  *  0.  Hence  by  selecting  Y 
as  the  charge  we  allow  interactions  not  only  between  samples 
of  ordinary  bulk  matter,  but  also  between  bulk  matter  and  kaons. 
Such  an  interaction  could  account  for  the  anomalies  reported 
by  Aronson  el  al '6,  and  for  this  reason  hypercharge  emerged 
as  the  natural  choice  for  the  charge  in  the  original  fifth  force 


FIG.  2  Sketch  of  the  Airy  method  for  determining  the  newtonian  constant. 
By  comparing  the  acceleration  g(z)  at  a  depth  z  with  the  value  g(0)  at  the 
surface,  the  shell  of  mass  AM  can  be  'weighed'.  See  text.  ref.  15  and 
equation  (4)  for  details. 
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theory.  As  our  primary  focus  here  is  on  macroscopic  experi¬ 
ments,  we  can  set  Y  =  B  and  write  Vs(r)  in  the  form 

D  D 

Vj(r)  =  +/2  — i— ^  e~r/A  (5) 

r 

where  /  is  a  new  fundamental  constant.  It  is  straightforward  to 
show  that  when  Vj(r)  in  equation  (S)  is  combined  with  the 
newtonian  potential  VN(r),  the  resulting  expression  for  V(r) 
has  the  form  given  in  equation  (2),  but  with  a  replaced  by  a,, 


where  f  =  /2/  G^m^  ,/*,  =  m,/  mH  and  mH  =  m( , H 1 ).  As  ( B./ /x, ) 
and  ( B,/m> )  vary  from  one  material  to  another,  ait  is  a  function 
of  the  compositions  of  the  interacting  materials.  It  then  follows 
from  equations  (3)  and  (6)  that  the  accelerations  of  two  samples 
j  and  /  towards  a  common  source  i  would  depend  on  the 
constants  a,j  and  atj-  respectively,  and  these  are  generally 
unequal.  If  i  denotes  the  Earth  then  the  accelerations  of  j  and 
/  towards  the  Earth  would  depend  on  their  compositions,  and 
this  could  be  tested  by  simultaneously  dropping  two  dissimilar 
objects  as  Galileo  is  supposed  to  have  done2123.  Such  an  experi¬ 
ment,  although  relatively  simple  in  principle,  is  extremely 
difficult  in  practice  and  has  been  carried  out  to  high  precision 
only  recently24'23,  as  we  discuss  below.  Eotvos  recognized, 
however,  that  an  equivalent  experiment  could  be  done  by  sus¬ 
pending  the  two  samples  j  and  /  from  a  common  bar,  and 
searching  for  a  torque  on  the  bar  arising  from  the  differential 
force  exerted  on  the  samples  by  the  Earth.  Subsequent 
experiments26'28  used  the  Sun  as  a  source,  and  these  set  limits 
on  possible  new  weak  forces  for  which  A  is  of  order  1  au. 

It  follows  that  in  the  model  considered  here,  the  acceleration 
difference  A«„-  of  j  and  j'  towards  the  Earth  is  given  by 

where  &  is  the  field  strength  (in  units  of  acceleration)  of  the 
source,  which  in  this  case  would  be  the  Earth  (denoted  by  ©). 
As  in  equation  (12)  below,  B  can  be  replaced  by  a  ino-e  general 
charge  Q,  so  that  (B/#x)@-»  (G/m)®  and  so  on.  T:t  then  follows 
that  the  experimentally  determined  sign  of  Aa^-  (that  is,  whether 
j  falls  faster  than  /  or  vice  versa)  depends  on  the  signs  of  the 
factors  and  [(Q/fi), -(Q/^), ■],  and  on  the  direction 

of  9.  It  should  be  emphasized  that  for  a  force  of  finite  range 
the  effective  values  of  (G/#x)®  and  &  are  actually  determined 
by  the  local  chemical  composition  and  the  local  matter  distribu¬ 
tion,  respectively,  and  can  vary  from  one  location  to  another. 
As  this  information  for  the  site  of  the  original  EPF  experiment 
is  not  accurately  known,  the  sign  and  magnitude  of  a  implied 
by  this  experiment  cannot  now  be  established  unambiguously. 

From  equation  (3)  we  see  that  with  a  -»  atJ,  an  intermediate- 
range  weak  force  can  be  detected  through  a  modification  of  the 
inverse-square  law  (G(r)  *  constant)  and/or  through  a  compo¬ 
sition-dependence  of  the  net  acceleration  (a0  *  constant).  (The 
latter  effect  is  often  referred  to  as  a  violation  of  the  weak 
equivalence  principle  or  of  the  universality  of  free-fall.)  It  is 
evident  from  equation  (3)  that  both  effects  are  generally  present 
in  a  given  experiment.  In  practice,  however,  experiments  are 
usually  designed  to  isolate  one  or  the  other  of  these  signals  for 
a  new  force,  and  we  will  therefore  focus  on  each  of  these  classes 
of  experiments  separately. 

Composition-Independent  experiments 

We  first  consider  the  recent  progress  made  in  tests  of  the  inverse- 
square  law.  Returning  to  Fig.  1,  we  see  that  the  limits  on  the 
strength  of  a  possible  non-newtonian  component  have 

improved  substantially  over  several  distance  scales.  Most  not¬ 
ably  the  geophysical  window  at  A  *»  100  m  has  been  partially 
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closed  by  a  series  of  experiments  using  geophysical  sources,  as 
we  now  discuss.  In  1988,  Eckhardt  and  collaborators29  intro¬ 
duced  a  new  technique  into  the  field,  by  measuring  the  gravita¬ 
tional  acceleration  g(z)  as  a  function  of  height  z  using  the  600-m 
tower  of  the  television  station  WTVD  in  Gamer,  North  Carolina. 
If  newtonian  gravity  is  correct  then  Laplace's  equation  can  be 
used  to  predict  g(z),  given  a  knowledge  of  g  on  the  surface  of 
the  Earth  near  the  tower.  Although  this  technique  is  in  principle 
straightforward,  it  had  never  been  attempted  before  1988, 
apparently  because  of  the  widespread  (and  erroneous)  belief 
that  g(2)  could  not  be  measured  with  sufficient  precision  on 
towers.  Although  measuring  g(z)  on  a  tower  is  indeed  more 
difficult  than  in  the  relatively  protected  environment  of  a  mine- 
shaft,  the  tower  method  has  an  important  compensating  advan¬ 
tage.  If  there  are  pockets  of  unusually  dense  material  below  the 
surface  of  the  Earth  near  the  experiment,  these  will  be  located 
on  average  farther  away  from  the  gravimeter  being  used  to 
determine  g(z)  in  a  tower  experiment  than  for  analogous 
measurements  in  a  mine  or  borehole.  The  distortion  by  such 
local  mass  concentrations  of  the  average  gravitational  field  near 
the  tower  will  have  less  effect  than  in  a  similar  mine  experiment. 
This  ‘smoothing*  of  the  average  gravitational  field  facilitates  a 
comparison  between  the  measured  and  predicted  values  of  g(z). 
Eckhardt  et  al  in  fact  demonstrated  not  only  the  feasibility  of 
tower  measurements,  but  also  the  ability  of  this  technique  to 
set  useful  limits  on  possible  deviations  from  newtonian  gravity. 

In  their  original  analysis29,  Eckhardt  era/,  found  a  discrepancy 
of  (-500  ±35)  p.Gal  (lp.Gal=  10-8  m  s“2)  between  the  measured 
and  predicted  values  of  g(z)  at  z  =  562.24  m.  The  negative  sign 
corresponds  to  the  effect  expected  from  a  new  attractive  non- 
newtonian  force,  which  Eckhardt  er  al.  termed  the  ‘sixth*  force. 
But  it  was  subsequently  pointed  out  by  Bartlett  and  Tew30  that 
the  surface  gravity  database  used  in  ref.  29,  which  was  compiled 
by  the  Defense  Mapping  Agency,  could  have  been  biased  by 
oversampling  of  the  higher  elevation  terrain  compared  with  that 
at  lower  elevation.  Low-lying  terrain  is  often  relatively  inaccess- 


FK5.  3  The  Eotvos  ‘single-arm’  apparatus,  shown  to  scale.  This  drawing  is 
based  in  part  on  details  supplied  in  ref.  19.  and  on  an  original  torsion  balance 
which  survives  in  the  museum  of  the  Geophysical  Observatory  in  Tihany. 
Hungary.  See  ref.  20  for  additional  background. 


ible,  as  would  be  the  case  if  a  stream  were  running  through  the 
region.  As  the  gravimeters  used  in  such  experiments  are  expen¬ 
sive,  it  is  not  surprising  that  fewer  measurements  would  be  made, 
for  instance,  at  the  edge  of  a  body  of  water,  than  at  a  nearby 
road  running  above  the  water.  A  difference  of  I  m  between  the 
average  elevation  of  the  actual  topography  and  that  of  the 
sampled  points  could  lead  to  a  discrepancy  of —309  p.Gal,  which 
is  comparable  to  the  magnitude  of  the  effect  observed  in  ref. 
29.  Eckhardt  et  al.  have  recently  re-examined  the  effects  of  the 
terrain  in  their  experiment,  and  now  find  a  result  consistent  with 
newtonian  gravity31.  Similar  null  results  have  also  been  obtained 
recently  in  two  other  tower  experiments32-33.  Speake  et  al?3 
devoted  considerable  effort  to  understanding  how  tower  motion 
affects  the  performance  of  the  LaCoste- Romberg  gravimeters 
used  in  these  experiments.  Their  careful  analysis  supported  the 
original  claim  by  Eckhardt  et  al.  that  g(z)  could  be  measured 
on  towers  with  sufficient  precision  to  allow  a  meaningful  com¬ 
parison  with  the  predictions  of  newtonian  gravity.  The  limits 
implied  by  these  experiments  are  partially  responsible  for  filling 
in  the  ‘geophysical  window'  in  Fig.  1. 

Another  class  of  geophysical  experiments  which  tests 
newtonian  gravity  over  the  same  distance  scales  comprises  the 
pumped  lake  experiments  of  Muller  et  al.34  and  Moore  et  at.33. 
in  these,  the  change  in  the  local  value  of  g  is  measured  as  a 
function  of  the  height  to  which  the  reservoir  is  filled,  and  the 
experiments  can  thus  be  viewed  as  ‘weighing’  a  slab  of  water 
whose  lateral  extent  is  hundreds  of  metres.  The  lake  experiments 
are  subject  to  a  number  of  systematic  uncertainties,  and  thus 
far  have  not  achieved  the  sensitivity  of  the  tower  experiments. 

Similarly,  Zumberge  et  aL “  very  recently  used  the  ocean  as 
an  attracting  mass  to  determine  the  newtonian  constant  G  for 
matter  interacting  over  a  distance  scale  of  —5,000  m.  This  type 
of  Airy  experiment,  although  technically  complex,  offers  a  num¬ 
ber  of  advantages  over  similar  measurements  in  mines,  lakes 
and  boreholes.  First,  ocean  experiments  allow  G  to  be  studied 
over  distance  scales  that  are  larger  than  those  readily  accessible 
in  land-based  experiments.  Second,  the  density  of  the  5,000-m 
shell  of  matter  (see  Fig.  2)  is  not  only  much  more  uniform  than 
in  a  mine,  but  varies  in  a  known  and  reasonably  smooth  way 
with  depth.  ( Mines,  after  all,  are  situated  where  they  are  precisely 
because  of  the  inhomogeneities  they  contain  in  the  form  of 
relatively  dense  ores.)  By  combining  various  sea-surface 
measurements,  Zumberge  et  al  determined  G  from  a  generalized 
version  of  equation  (4).  Their  result,  G  =  (6.677  ±0.013)  x  l(Tn 
N  m2  kg'2  agrees  with  the  laboratory  value. 

This  experiment  helps  fill  in  a  large  part  of  the  geophysical 
window,  and  sets  useful  constraints  down  to  values  of  A  =  0.1  m. 
As  noted  previously,  experiments  typically  give  their  most  strin¬ 
gent  limits  for  values  of  A  comparable  to  the  characteristic 
separation  of  the  test  masses  in  that  experiment.  It  might  then 
seem  surprising  that  an  experiment  carried  out  on  a  distance 
scale  of  —5,000  m  can  give  a  useful  limit  down  to  —0.1  m.  But 
the  constraint  in  Fig.  1  is  obtained  by  comparing  the  Zumberge 
value  for  G  to  the  Cavendish  value  of  G0,  which  is  itself 
determined  over  a  scale  of  a  few  centimetres.  This  discussion 
applies  (with  appropriate  changes)  to  the  mine/borehole  and 
lake  experiments. 

The  other  regimes  in  which  the  limits  on  £  have  improved 
are  at  small  A  (laboratory  scales)  and  at  larger  A  (satellite  and 
astronomical  scales).  At  present  the  best  laboratory  limits  come 
from  the  null  experiments  of  Hoskins  et  al 37  and  Chen  et  alis 
These  experiments  arose  in  part  from  work  by  Long39,  who  in 
1976  claimed  to  see  deviations  from  the  inverse-square  law.  To 
illustrate  the  novel  technology  used  in  some  of  these  experi¬ 
ments,  we  will  briefly  describe  the  null  experiment  of  Hoskins 
et  al,  which  tests  Newton’s  law  in  the  2-5-cm  range.  These 
authors  note  that  if  the  1/r2  law  is  exact  then  a  test  mass  located 
inside  an  infinitely  long  hollow  cylinder  would  experience  no 
gravitational  force,  just  as  if  it  were  inside  a  spherical  shell. 
Although  there  are  corrections  for  end  effects  in  a  finite  cylinder. 
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these  are  small  enough  to  allow  detection  of  the  force  that  would 
arise  from  a  deviation  from  the  1/r2  law.  By  monitoring  the 
torque  on  a  test  mass  inside  the  cylinder  as  the  distance  between 
the  mass  and  the  cylinder  is  varied,  limits  can  be  set  on  possible 
non-newtonian  forces.  Another  novel  technology  is  the  laplacian 
detector  of  Paik40,41,  which  directly  tests  the  implication  of  the 
1/r*  law  that  the  gravitational  potential  <t>(r)  is  a  solution  of 
the  equation  V2d>(r)  =  0.  Although  the  existing  limits  from  this 
experiment  are  less  stringent  than  those  of  refs  37  and  38,  this 
technology  offers  the  prospect  of  considerable  improvements  in 
sensitivity.  Moreover,  the  laplacian  detector  can  be  used  to  study 
newtonian  gravity  in  space,  and  efforts  along  this  line  are  cur¬ 
rently  under  way. 

For  values  of  A  in  the  range  106-108  m  there  has  been  a 
considerable  improvement  in  £  which  resulted  from  a  com¬ 
parison  by  Rapp42  of  the  values  GM ®  inferred  from  Earth- 
surface  gravity  measurements,  and  from  a  study  of  the  orbits 
of  the  LAG  EOS  satellite  and  of  the  Moon  about  the  Earth.  The 
improved  limits  near  A  =  10"  m  come  from  the  analysis  in  ref. 
9  of  planetary  data  on  G(r),  as  well  as  data  on  anomalous 
planetary  precessions.  To  understand  how  these  data  can  be 
used  to  set  limits  on  a  and  A,  we  recall  from  equation  (3 )  that 
in  the  presence  of  a  new  force  the  newtonian  constant  Gx  is 
replaced  by  the  function  G(r).  It  follows  that  under  these 
circumstances  Kepler's  third  law  assumes  the  form 

Op=  G(ap)M&(Tp/2ir)2  (8) 

where  Tp  is  the  period,  and  ap  is  the  physically  measured 
semi-major  axis,  of  planet  p.  Given  a  set  of  values  of  ap  and  Tp 
for  the  various  planets,  G(ap)Me  can  be  determined  and  its 
constancy  tested.  The  curves  labelled  'Planetary'  in  Fig.  1  derive 
from  an  even  more  sensitive  test,  namely  the  anomalous  pre¬ 
cession  of  the  perihelion  of  the  orbit  induced  by  a  non-newtonian 
coupling.  For  the  Yukawa  interaction  in  equation  (2)  the 
anomalous  precession  Sip,  is  given  by 

8<P,  s  irat(a/ A)2  e~u/4  (9) 

where  a  is  the  mean  value  of  the  semi-major  axis  of  the  orbit. 
Because  a  precession  of  the  perihelion  is  also  predicted  by 
general  relativity  (GR),  one  can  extract  a  limit  on  Sip,  (and 
thereby  on  a  and  A)  only  by  assuming  that  GR  is  correct,  for 
which  ample  independent  evidence  exists43. 

In  summary,  much  progress  has  been  made  since  1981,  but 
additional  work  needs  to  be  done  over  the  distance  scales  in 
the  geophysical  window.  At  present  Eckhardt  and  coworkers 


are  in  the  process  of  doing  a  new  experiment  using  the  2,000-m 
WABG  television  tower  near  Indianola,  Mississippi.  This  site 
was  selected  in  part  because  the  flatness  of  the  surrounding 
terrain  should  help  to  minimize  the  problem  of  terrain  bias 
encountered  in  their  earlier  work.  Working  in  the  other  direction, 
the  terrain  problem  noted  by  Bartlett  and  Tew30  apparently 
accounts  for  most  of  the  anomalous  results  originally  reported 
by  Stacey  ei  al.  But  results  reported  by  Ander  el  al.  from 
Greenland  require  further  study  to  determine  whether  the 
apparent  anomalous  accelerations  g(  z )  in  that  experiment  could 
arise  from  subsurface  concentrations  of  relatively  dense 
material. 

Composition-dependent  experiments 

We  turn  next  to  describe  the  recent  progress  in  composition- 
dependent  experiments.  Before  1986,  the  only  reliable  data  came 
from  the  original  experiment  of  Eotvos,  Pekar  and  Fekete 
(EPF)*1''20,  which  primarily  compared  the  accelerations  of 
various  pairs  of  materials  falling  to  the  Earth,  and  from  three 
later  experiments26  28  which  used  the  Sun  as  a  source.  The 
Eotvos  apparatus  (Fig.  3)  was  used  to  measure  the  differential 
torque  on  a  torsion  balance  for  1 1  pairs  of  chemically  different 
materials'9.  As  the  experiment  was  originally  intended  to  test 
the  principle  of  universality  of  free  fall,  it  was  assumed  that  any 
acceleration  difference  between  the  test  samples  would  arise 
from  the  force  exerted  on  these  materials  from  the  Earth  as  a 
whole.  But  it  was  noted  in  ref.  1  that  their  results  could  also  be 
used  to  set  limits  on  interactions  arising  from  a  local  source 
(such  as  a  mountain,  cliff,  basement),  if  the  chemical  composi¬ 
tion  and  matter  distribution  of  the  source  were  sufficiently  well 
known.  References  7,  19  and  20  describe  the  original  EPF 
experiment  in  more  detail. 

The  original  Eotvos  experiment  can  be  used  to  set  limits  on 
intermediate-range  composition-dependent  forces,  but  because 
various  details  of  this  experiment  remain  unknown  to  us,  these 
limits  are  necessarily  somewhat  model-dependent.  Thus,  as  a 
practical  matter,  the  only  unambiguous  limits  on  composition- 
dependent  forces  before  1986  applied  to  those  for  which  A  » 
1  au,  which  derive  from  the  experiments  of  refs  26-28.  We  see 
from  Fig.  4  that  the  whole  distance  scale  from  a  few  centimetres 
to  I  au  has  now  been  filled  in  by  high-precision  experiments 
using  both  laboratory  and  geophysical  sources.  (As  in  Fig.  1, 
we  show  only  the  envelope  of  the  limits  derived  from  the  various 
experiments,  and  not  the  individual  results.)  Most  of  the  new 
experiments  use  torsion  balances  similar  to  those  developed  by 
Eotvos  to  compare  the  accelerations  of  their  test  masses.  We 
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FIG.  4  Constraints  on  the  coupling  constant  f  as  a  function  |()  s 

of  the  range  A  from  composition-dependent  experiments, 

assuming  a  coupling  to  B.  With  the  exception  of  the  region 

denoted  by  'Solar  Eotvos'.  all  the  constraints  shown  have  w>  I0'6 

been  obtained  since  1986,  and  most  of  the  results  shown  + 

are  from  the  Eot-Wash  experiments45.  The  hatched  region  7 

denotes  values  of  A  where  the  limits  on  f  depend  on  detailed 

models  of  the  Earth  and  are  at  present  uncertain.  The  dotted 

contour  is  our  'best  guess'  for  the  constraint  in  this  region.  I0"8 
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show  in  Fig.  5  the  apparatus  used  Adelberger  el  al.4>  (known 
as  the  Edt-Wash  collaboration),  who  have  carried  out  the  most 
extensive  searches  for  composition  dependent  effects  over  a 
wide  range  of  distance  scales.  As  can  be  seen  in  Fig.  S,  all  of 
the  test  masses  in  the  Edt-Wash  balance  are  located  in  the  same 
vertical  position,  in  contrast  to  the  case  for  the  Eotvos  balance. 
The  Eotvos  apparatus  was  originally  designed  to  measure  gravity 
gradients,  which  required  a  vertical  separation  of  the  masses. 
This  separation  is,  however,  undesirable  in  searches  for  compo¬ 
sition-dependent  effects,  as  it  produces  a  substantial  systematic 
bias  which  must  be  compensated  for.  The  Eot-Wash  balance 
introduced  other  innovations,  such  as  compensating  masses  to 
cancel  (vertical)  gravity  gradients,  a  high-precision  turntable  to 
rotate  the  apparatus  smoothly  and  an  optical  read-out  system. 
Adelberger  and  collaborators  have  compared  the  accelerations 
of  different  pairs  of  test  masses  towards  various  sources,  includ¬ 
ing  a  large  mass  of  lead  bricks,  the  hillside  adjacent  to  their 
apparatus,  and  the  Earth.  In  all  cases  their  results  have  been 
compatible  with  the  weak  equivalence  principle,  which  is  to  say 
that  they  see  no  evidence  for  a  composition-dependent  fifth 
force.  As  an  example,  their  most  recent  hillside  limit  on  the 
acceleration  difference  Aox  of  aluminium  and  beryllium  is4S 

AaJBe-AI]  =  [(-2.0±2.2)E  +  (0.9±2.1)JV] 

x  10'"  cm  s'2  (10) 

where  E  and  N  are  unit  vectors  which  point  east  and  north 
respectively.  For  further  details  of  the  Edt-Wash  experiments, 
see  ref.  45. 

A  complementary  design  for  a  torsion  balance  has  been  intro¬ 
duced  by  Boynton  el  aL46.  These  authors  measure  the  period 
7(0)  of  a  ring  composed  of  two  dissimilar  materials  oscillating 
in  a  horizontal  plane,  as  a  function  of  the  initial  orientation 
angle  0  of  the  ring  with  respect  to  their  source,  a  cliff  located 


FIG.  5  Diagram  of  the  Eot-Wash  torsion  balance  of  Adelberger  ef  at.  adapted 
from  ref.  45.  In  contrast  to  the  original  Eotvos  apparatus  (Fig.  3),  the 
Edt-Wash  balance  has  the  different  test  masses  very  symmetrically 
arranged.  The  experimental  site  is  the  Nuclear  Physics  Laboratory,  located 
on  a  hillside  at  the  University  of  Washington  in  Seattle. 
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near  Mount  Index  in  Washington.  It  can  be  shown7  that  a 
dependence  of  T(0)  on  the  orientation  of  their  ring  (which  is 
a  composition  dipole)  is  a  signal  for  a  composition-dependent 
fifth  force,  just  as  the  shift  in  the  equilibrium  position  of  the 
balance  is  for  the  EPF  and  Edt-Wash  experiments.  The  first 
results  published  by  Boynton  el  aL  found  a  marginally  significant 
signal  that  could  be  interpreted  as  evidence  for  a  fifth  force. 
Their  signal  was  a  nonzero  value  of  the  difference  A7(0)/  7(0) 
where  A 7(0)  =  7(0)-  7(0  +  ir),  and  where  0  =  0  corresponded 
to  the  aluminium  half  of  their  beryllium-aluminium  ring  being 
adjacent  to  the  cliff.  Boynton  et  aL  found  for  A7( 0) 

A  7(0)/  7(0)  =  (-4.6±l.l)xlO‘cos0 

(11) 

+  (+0.1  ±  1.2)  sin  0 

But  a  later  experiment47  using  a  copper-polyethylene  ring  found 
no  evidence  for  a  fifth  force.  It  is  not  yet  clear  how  we  should 
interpret  the  result  in  equation  (11),  and  Boynton  et  aL  are 
continuing  their  efforts  with  an  improved  apparatus.  The  result 
in  equation  ( 1 1)  is  not  compatible  with  the  more  stringent  limits 
of  experiments  that  obtained  null  results  (see  Table  1  of  ref.  3). 

in  addition  to  variants  of  torsion  balances,  a  number  of  novel 
techniques  have  also  been  introduced,  such  as  the  floating-ball 
experiments  of  Thieberger48  and  of  Bizzeti  et  aL*9,  and  the 
free-fall  experiments  of  Niebauer  et  aL 24  and  Kuroda  and  Mio25. 
Thieberger’s  experiment  is  shown  schematically  in  Fig.  6.  The 
central  feature  is  a  hollow  copper  ball  floating  in  water.  A 
composition-dependent  fifth-force  F3  (arising  in  his  experiment 
from  a  cliff  at  the  Palisades,  New  Jersey)  would  act  differently 
on  the  test  masses,  which  are  the  copper  shell  and  the  water  it 
displaces.  This  unbalanced  force  causes  the  ball  to  move,  and 
from  its  velocity  we  can  infer  the  magnitude  of  F; .  The  advantage 
of  such  an  apparatus  is  that  the  high  degree  of  symmetry  of  the 
test  masses  makes  this  experiment  relatively  insensitive  to  gravity 
gradients,  which  could  simulate  the  presence  of  a  fifth-force.  In 
practice,  the  symmetry  of  Thieberger’s  apparatus  is  broken  by 
a  small  pin  which  protrudes  from  the  sphere  and  which  is  needed 
to  make  the  sphere  buoyant.  Bizzeti  et  aL  have  carried  out  a 
similar  experiment  (in  a  monastery  near  Vallambrosa,  Italy) 
with  an  apparatus  of  even  higher  symmetry.  They  use  a  solid 
nylon  sphere  floating  in  a  solution  containing  various  salts,  and 
the  ball  is  kept  buoyant  by  a  gradient  in  the  density  of  the 
solution. 

Historically,  Thieberger’s  experiment  was  the  first  test  of  the 
fifth  force  hypothesis  to  search  for  a  composition-dependent 
effect.  In  fact  Thieberger  detected  such  an  effect,  in  the  form  of 
a  steady  drift  of  his  sphere  across  the  tank.  The  average  velocity 
of  his  sphere,  v  =  (4.7  ±0.2)  mm  h'1,  could  be  accounted  for  by 
the  model  in  equations  (5)-(7)  with  aA  =  ( 1 .2  ±  0.4)m.  By  con¬ 
trast,  Bizzeti  et  aL  have  set  an  upper  limit  of  v<  0.010 
mm  h_,( lcr)  on  a  drift  of  their  sphere,  and  the  corresponding 
limit  for  a  coupling  to  B  is  fA  <  0.30  m  (lcr)  for  A  ^  1  km. 

Another  interesting  route  which  is  being  explored  is  the  use 
of  gravity  wave  antennas  as  fifth  force  detectors 50,51 .  In  simple 
conceptual  terms,  such  an  experiment  consists  of  a  cylindrical 
detector  and  a  nearby  rotating  composition  dipole  whose  axis 
passes  through  the  centre-of-mass  of  the  rotating  masses.  In  the 
absence  of  a  fifth  force  the  rotor  would  drive  the  detector  at 
twice  the  frequency  w  of  rotation  (assuming  that  the  axis  of 
rotation  is  precisely  aligned  with  the  centre-of-mass  of  the  rotor), 
in  the  presence  of  a  composition-dependent  fifth  force,  there 
would  be  an  additional  force  at  w,  and  if  w  corresponds  to  the 
resonant  frequency  of  the  bar,  then  limits  can  be  set  on  £.  At 
present  these  limits  are  not  yet  competitive  with  those  from 
other  experiments,  but  as  the  technology  of  gravity-wave  detec¬ 
tors  improves,  so  will  the  sensitivity  of  such  experiments. 

In  addition  to  composition-dependent  experiments  that  use 
geophysical  sources,  other  experiments  constrain  £  for  both 
smaller  and  larger  values  of  A.  Using  Fig.  4,  we  can  summarize 
the  present  situation  as  follows.  For  small  values  of  A,  experi- 
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FIG.  6  Schematic  drawing  of  the  forces  on  a  hollow  copper  shell  floating  in 
water.  In  the  absence  of  a  composition-dependent  fifth  force  F,,  the  weight 
Mg  is  supported  by  the  buoyant  force  F„.  The  presence  of  F,  (due  to  an 
external  source)  leads  to  an  unbalanced  force  which  moves  the  shell.  See 
refs  48  and  49  for  more  details. 

merits  have  typically  involved  laboratory-size  sources  construc¬ 
ted  using  blocks  of  lead45-52'55.  The  obvious  advantage  of  such 
experiments  is  that  the  source  has  known  dimensions  and  com¬ 
position,  but  their  disadvantage  is  that  the  sensitivity  of  such 
experiments  drops  quickly  as  the  range  of  the  putative  fifth  force 
increases.  As  noted  above,  for  larger  A  the  best  limits  on  £  come 
from  experiments  using  geophysical  sources,  such  as  a  canal 
lock56,  a  mountain  or  a  cliff45’46'4”’49,57.  These  have  the  advantage 
of  using  much  larger  sources  than  laboratory  experiments,  albeit 
ones  whose  compositions  are  less  well  known. 

As  A  increases,  there  is  a  region  in  which  the  best  limits  on 
£  come  from  Galileo-type  experiments24'25,  which  compare  the 
accelerations  of  two  test  masses  falling  freely  towards  the  Earth. 
These  experiments  are  inherently  less  sensitive  than  those  using 
torsion  balances,  but  are  in  practice  easier  to  analyse  in  the 
region  104  m«  A  «  106  m.  This  is  a  consequence  of  the  fact  that 
Galileo  experiments  are  sensitive  to  the  component  of  force 
parallel  to  g,  (‘Fj’)  whereas  torsion  balance  experiments  are 
necessarily  sensitive  only  to  the  component  of  F  perpendicular 
to  g.  Over  this  intermediate  range,  F„  is  significantly  less  sensitive 
to  fine  details  of  the  underlying  structure  of  the  Earth  than  is 
F± ,  and  hence  is  much  more  precisely  known.  As  the  Galileo 
experiments  can  be  interpreted  with  much  higher  certainty  than 
the  torsion  balance  experiments,  they  consequently  give  rise  to 
the  best  limits  in  this  region  of  A.  Further  improvements  in  the 
analysis  of  Fx  over  these  distance  scales  may  eventually  allow 
limits  from  torsion  balance  measurements  to  surpass  those  from 
Galileo  measurements,  which  are  indicated  by  the  hatched 
region  in  Fig.  4. 

As  A  increases  still  further,  terrestrial  Eotvos  experiments 
again  surpass  free-fall  experiments  in  the  limits  they  give  for  £, 
because  F±  now  depends  only  on  global  properties  of  the  struc¬ 
ture  of  the  Earth.  Finally,  for  A  *»  1  au  the  Eotvds-type  experi¬ 
ments  using  the  Sun  as  a  source  are  the  most  sensitive,  as  we 
noted  earlier.  We  see  that  at  present,  the  distance  scale  from 
roughly  1  m  to  1  AU  is  now  covered  by  overlapping  experiments 
using  a  variety  of  technologies. 

When  studying  composition-dependent  effects  it  is  also 
necessary  to  take  account  of  the  possibility  that  a  given  pair  of 
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materials  j  and  /  may  coincidentally  have  the  same  values  of 
B/fi  (ratio  of  baryon  number  to  m/mH),  and  hence  have  the 
same  accelerations  even  in  the  presence  of  a  fifth  force.  Related 
to  this  is  the  possibility  that  the  charge  that  determines  the 
strength  of  the  fifth  force  interaction  is  not  necessarily  B,  but 
may  be  some  linear  combination  of  the  possible  quantum  num¬ 
bers  which  describe  a  sample  of  bulk  matter.  If  N  and  Z  denote 
the  numbers  of  neutrons  and  protons  respectively,  then  the  most 
general  additive  charge  Q  is 

Q  =  cos  0sB  +  sin  8SIZ  (12) 

where  (as  before)  B  =  N  +  Z,  Iz  =  N  -Z  and  05  is  an  appropri¬ 
ate  ‘mixing  angle’.  The  previous  formalism  can  then  be  taken 
over  with  B/fi  replaced  by  Q/n.  It  follows  from  equation  (12) 
that  at  least  two  pairs  of  materials  must  be  compared  to  rule 
out  the  possibility  of  an  accidental  cancellation  in  the  charges 
of  the  test  masses,  and  this  has  been  checked  explicitly  by  the 
Eot-Wash  collaboration45.  It  is  also  possible  that  for  some 
choices  of  Q  (for  example  Qa  lz)  the  source  strength  (<?/ft)i 
could  also  vanish.  For  Q=  lz  this  is  an  interesting  possibility, 
because  most  geophysical  sources  have  near-zero  values  of  lz. 
(An  exception  is  water,  which  was  the  source  in  the  experiment 
of  Bennett56.)  The  possibility  that  some  geophysical  experiments 
had  obtained  null  results  because  Q  =  lz  =  0  was  excluded  by 
various  laboratory  experiments45,52'55  which  capitalized  on  the 
fact  that  a  lead  source  has  a  relatively  large  value  of  lz. 

Elementary  particle  experiments 

The  energy-dependent  effects  reported  by  Aronson  et  al ,16  which 
provided  part  of  the  motivation  for  the  fifth-force  hypothesis, 
have  also  been  tested  experimentally.  Grossman  et  al.  searched 
for  an  energy  dependence  of  the  Ks  lifetime  rs  *=  h/l's  by  com¬ 
paring  its  value  measured  in  the  momentum  range  (100-350) 
GeV/c  to  the  laboratory  value  determined  in  the  range  (2-5) 
GeV/c.  Their  result,  r*  =  (0.8920 ±  0.0044) x  10~I0s,  is  in 
excellent  agreement  with  the  laboratory  value  rs  = 
(0.8922  ±  0.0020)  x  10" 10  s.  Similarly,  Coupal  et  aL 59  found  |  = 

(2.28  ±  0.06)  x  10' 3  in  the  momentum  range  (30-150)  GeV/c, 
which  agrees  well  with  the  low-energy  laboratory  value  |n±l  = 
(2.268 ± 0.023)  x  10~3.  But  as  noted  in  ref.  60  the  Coupal  result 
for  |tj*|  is  also  consistent,  within  their  quoted  errors,  with  that 
of  Aronson  et  at16.  Recently,  Carosi  et  aL 61  have  reported 
high-energy  determinations  of  4>*  and  </>oo.  They  obtain 
<t>±  =  46.9° ±1.4° ±0.7°,  <£00  =  47.1°±2.1°±  1.0°,  70GeV*sEK« 
170  GeV,  where  the  first  error  is  statistical  and  the  second  is 
systematic.  These  results  are  in  good  agreement  with  the  corres¬ 
ponding  low-energy  values62. 

Another  implication  of  the  original  model  of  a  coupling  to 
hypercharge  has  also  been  called  into  question  by  recent  experi¬ 
ments.  This  is  the  suggestion63'66  that  the  quanta  of  this  field, 
called  hyperphotons  (yY),  could  be  detected  in  decays  such  as 
K+  -» ir+  +  yY .  The  signal  for  such  a  decay  would  be  the  detection 
of  a  it*  with  an  appropriate  momentum,  |p„|  =  227  MeV/c,  and 
no  other  detected  particles.  This  is  a  fairly  clean  signal  which 
has  been  looked  for  in  several  recent  experiments.  The  current 
best  limits  on  this  branching  ratio  are 

T(K4  -» it*  +  undetected  neutrals) 
r(K+-»all  final  states) 

(13) 

4.6  x  10  8, 90%  confidence  limits6 
6.4  x  10'\  90%  confidence  limits6* 

Several  theoretical  calculations  have  used  the  experimental  limit 
in  equation  (13)  above  to  infer  limits  on  the  strength,  f2,  of  a 
possible  coupling  to  hypercharge.  Although  the  results  of  these 
calculations  differ  considerably,  for  a  coupling  to  hypercharge 
(  Y  =  B  +  S)  they  all  agree  that  the  inferred  values  of  /2  are 
probably  too  small  to  be  compatible  with  either  the  energy- 
dependence  of  the  kaon  parameters  claimed  by  Aronson  et  al. 
or  with  the  value  of  /2  required  to  explain  the  Eotvos  data19. 
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We  emphasize  that  this  conclusion  applies  specifically  to  the 
case  where  yv  is  a  vector  particle  and  couples  to  hypercharge; 
for  a  scalar  field  there  are  no  meaningful  limits  from  kaon  decay. 
Similarly,  if  yY  couples  to  some  more  general  linear  combination 
of  B  and  5,  the  constraints  arising  from  K*  decay  are  not  in 
direct  conflict  with  other  data  (see  ref.  60  for  more  details). 

Another  class  of  elementary-particle  experiments  that  can  be 
used  to  explore  models  of  the  putative  fifth  force  consists  of 
those  comparing  the  free-fall  accelerations  of  particles  and  their 
antiparticles.  In  the  1960s  Fairbank  and  coworkers64  attempted 
to  measure  the  gravitational  accelerations  of  e*  and  e~,  but 
these  efforts  were  thwarted  by  the  presence  of  residual  elec¬ 
tromagnetic  fields  produced  by  surface  effects  in  their  apparatus. 
Goldman,  Hughes  and  Nieto  (GHN)70  subsequently  noted  that 
such  effects  would  be  relatively  less  important  in  an  experiment 
comparing  the  accelerations  of  p  and  p,  because  of  the  larger 
masses  of  these  particles.  The  impetus  for  carrying  out  such  an 
experiment  increased  following  the  suggestion  of  a  possible  fifth 
force,  as  the  proposed  (vector)  hypercharge  field  would  be 
expected  to  produce  a  difference  in  the  apparent  gravitational 
accelerations  of  p  and  p.  In  most  models  this  difference  would 
be  too  small  to  be  detected,  but  GHN  also  noted  that  some 
models  containing  both  scalar  and  vector  fields  predicted  an 
effect  large  enough  to  be  measurable. 

This  claim  has  been  called  into  question  by  several  recent 
analyses71'74.  As  was  first  noted  by  Schiff75,  if  matter  and  anti¬ 
matter  behaved  differently  in  a  gravitational  field,  this  could 
show  up  as  an  anomaly  in  the  Eotvos  experiment.  Using  a  similar 
argument,  Adelberger  el  aL7'  demonstrate  that  a  fifth-force 
model  that  could  produce  an  acceleration  difference  between  p 
and  p  would  also  produce  an  acceleration  difference  between 
the  test  masses  in  a  typical  Eotvos  experiment.  Given  the  extra¬ 
ordinary  sensitivity  of  the  current  Eotvos  experiments,  Adelber¬ 
ger  el  at.  argue  that  existing  limits  from  their  Eot-Wash  experi¬ 
ments  preclude  seeing  any  acceleration  difference  between  p 
and  p  at  the  sensitivity  level  expected  in  these  experiments.  See 
refs  71-74  and  76  for  further  discussion. 

Pushing  the  limits 

The  experimental  searches  for  both  composition-independent 
and  composition-dependent  deviations  from  newtonian  gravity 
have  made  enormous  advances  since  1986,  as  shown  in  Figs  1 
and  4.  No  compelling  evidence  has  yet  emerged  that  would 
indicate  the  presence  of  a  fifth  force,  although  the  anomalies 
reported  in  the  original  Eotvos  experiment  remain  to  be  under¬ 
stood,  as  do  those  in  the  experiments  of  Thieberger48  and 
Boynton  el  al.*6.  On  the  experimental  side,  efforts  continue  to 
set  even  more  stringent  limits  on  possible  deviations  from 
newtonian  gravity,  motivated  in  part  by  the  recognition  that 
such  experiments  may  be  our  most  powerful  tool  in  exploring 
physics  at  the  Planck  scale.  The  connection,  to  which  we  have 
alluded  earlier,  between  such  experiments  and  physics  at  the 
Planck  scale  has  been  explored  in  several  models70-77'81.  An 
excellent  review  of  the  wide  class  of  theories  that  predict  the 
existence  of  weak  macroscopic  forces  has  been  given  recently 
by  Fujii81.  Although  it  would  be  difficult  to  summarize  the 
content  of  these  models  briefly,  the  main  lesson  is  simply  that 
there  are  many  theories  that  predict  such  forces.  It  follows  that 
if  additional  weak  forces  are  not  seen,  as  is  the  indication  from 
present  experiments,  then  the  assumptions  behind  a  wide  class 
of  theories  will  have  been  called  into  question. 

The  experiments  we  have  discussed  can  be  understood  as 
extending  the  maximum  energy  scale  that  can  be  explored  from 
4x  104  GeV  (at  the  proposed  Superconducting  Super-Collider), 
to  — 1014  GeV,  an  energy  beyond  the  capability  of  any  acceler¬ 
ator.  Viewed  in  this  way,  torsion  balances,  floating  balls  and 
tall  towers  may  be  the  ultimate  high-energy  physics  tools.  □ 

Ephraim  Fischbach  and  Carrick  Talmadge  are  at  the  Physics  Department 
Purdue  University.  West  Lafayette,  Indiana  47907.  USA. 
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ABSTRACT 

The  Phillips  Laboratory  and  Purdue  University  are  conducting  a  tower 
gravity  experiment  near  the  town  of  Inverness,  MS.  Gravity  is  measured  at  six 
elevations  on  the  610  m  WABG-TV  tower  as  well  as  on  the  surface  in  an  8  km 
radius  about  the  tower.  These  data  are  combined  with  archived  data  extending  to 
300  km.  Using  previously  devised  techniques,  the  surface  data  are  analytically 
continued  and  compared  with  the  observations.  The  current  difference  at  the 
highest  tower  elevation  surveyed  so  far,  493  m,  is  34  pGal. 
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INTRODUCTION 


In  December  1991  the  Phillips  Laboratory  (PL)  Geophysics  Directorate 
(formerly  AFGL)  published  the  final  results  of  a  search  for  non-Newtonian 
gravity.  The  experiment  conducted  on  the  WTVD  tower  in  Clayton,  NC,  led  to  a 
null  result.  Many  difficulties  were  encountered  in  the  North  Carolina  experiment, 
not  the  least  of  which  were  systematic  effects  due  to  improperly  modeled  terrain. 
This  led  to  a  bias  which  resulted  from  the  fact  that  the  gravity  survey  elevations 
were  not  representative  of  the  actual  terrain.  These  problems  were  compounded 
by  the  sparsity  of  gravity  data  between  5  and  10  km  from  the  tower,  and  by  the 
inaccessibility  of  those  areas  due  to  dense  regions  of  trees.  Because  there  were 
some  lingering  uncertainties  involving  the  final  WTVD  results,  we  embarked  on  a 
follow-on  tower  experiment  in  an  area  where  all  known  uncertainties  could  be 
minimized. 

After  a  lengthy  search,  PL  selected  the  610  m  WABG-TV  tower  located 
northeast  of  Inverness,  MS.  Inverness  is  about  50  km  east  of  the  Mississippi  River 
and  about  350  km  north  of  the  Gulf  of  Mexico.  The  area  surrounding  the  tower  is 
extremely  flat  out  to  a  distance  of  40  km.  Also,  the  area  is  free  of  any  type  of 
forests,  permitting  gravity  measurements  to  be  made  at  most  desired  locations. 
Furthermore,  the  existing  gravity  data  in  the  area  are  much  more  extensive  than 
that  in  Clayton,  NC.  We  believe  that,  given  these  advantages,  we  shall  be  able  to 
resolve  many,  if  not  all,  of  the  previous  difficulties. 

SURFACE  GRAVITY  SURVEY  AND  ANALYTIC  CONTINUATION 

Based  on  existing  data  in  the  tower  area,  supplied  by  the  Defense  Mapping 
Agency  (DMA),  an  inner  zone  survey  out  to  8  km  was  deemed  sufficient.  The 
survey  plan  called  for  a  set  of  concentric  rings  with  up  to  ten  points  in  each  ring. 
The  ring  spacing  was  such  that  each  of  the  inner  ten  rings  contributed  nominally 
equal  weight  in  the  analytic  continuation  at  the  top  of  the  tower.  This  led  to  ring 
radii  of  150,  300,  450,  600,  800,  1050,  1400,  1900,  2600,  and  3600  m  from  the 
tower.  The  remaining  rings  were  placed  at  distances  of  4900,  6400,  and  8100  m. 
PL  and  Purdue  began  the  near-tower  gravity  survey  in  the  fall  of  1991.  The 
points  were  positioned  using  a  combination  of  the  Global  Positioning  System  (GPS) 
and  trigonometric  leveling  using  an  Electronic  Distance  Meter  (EDM). 
Positioning  all  the  points  using  GPS  was  not  possible  due  to  interference  from  the 
transmitter  near  the  tower1).  The  EDM  points  were  positioned  to  an  accuracy 
(relative)  of  1  m  in  the  horizontal  and  2  cm  in  the  vertical;  the  GPS  points  were 
accurate  (relatively)  to  3  cm  in  the  horizontal  and  4  cm  in  the  vertical.  A  total  of 
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35S  observations  contained  in  25  survey  loops  resulted  in  123  gravity  points 
surveyed  in  the  region  over  a  period  of  six  months.  A  least-squares  adjustment 
was  performed  on  the  data  along  with  corrections  for  earth  tide,  gravimeter  drift, 
and  scale  factor.  The  resultant  rms  error  is  13  pGal  with  no  individual  errors 
greater  than  30  pGal. 


LONGITUDE  (KM) 

Fig.  1  Survey  determined  elevations  minus  digitized  elevations 
contoured  at  0.5  m.  Dots  represent  locations  of  survey  points. 

The  next  phase  was  the  removal  of  the  high  wave  number  component  of  the 
gravity  field  from  the  surface  measurements.  Using  U.S.  Geological  Survey 
topographic  maps,  with  a  reported  accuracy  of  76  cm,  we  digitized  the  elevations 
inside  a  10  x  10  km  region.  We  interpolated  the  digital  terrain  to  die  points  of  our 
survey,  determined  elevations,  and  compared  the  results.  The  USGS  maps  are 
good,  but  not  perfect;  the  rms  difference  between  our  elevations  and  those  of  the 
maps  is  85  cm  (Figure  1).  Especially  striking  is  the  effect  of  the  catfish  farms 
(large  pools  of  raised  earth  where  local  farmers  breed  catfish)  just  south  of  the 
tower  that  do  not  appear  on  the  USGS  maps  and  thus  show  up  as  large  elevation 
differences.  Using  the  comparison  results,  we  corrected  the  digitized  elevations 
and  computed  terrain-corrected  Bouguer  anomalies  for  all  points  inside  of  10  km. 
Due  to  the  benign  nature  of  the  surrounding  terrain  the  rms  terrain  correction  to 
the  Bouguer  anomalies  was  a  mere  0.4  uGal.  The  resultant  gravity  field  is  very 
smooth  with  an  anomalous  horizontal  gradient  of  about  6.7  E  in  contrast  to  about 
13  E  in  the  North  Carolina  area  (Figure  2). 
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LONGITUDE  (KM) 

Fig.  2  Terrain -corrected  Bouguer  anomalies.  Dots  show  the 
location  of  survey  points;  contour  interval  is  1  mGal. 

We  merged  our  corrected  data  with  archived  data  obtained  from  DMA 
which  were  accurate  to  1-2  mGal  and  extended  to  300  km  from  the  tower.  We 
selected  7781  points  to  be  used  in  the  analytic  continuation  from  a  total  of  50292. 
We  obtained  digital  data  from  DMA  for  use  in  terrain  correcting  all  gravity 
points,  but  the  voluminous  amount  of  data  made  for  large  computational  and 
storage  requirements.  So,  given  that  the  terrain  corrections  inside  10  km  were 
very  small,  we  computed  simple  Bouguer  anomalies  for  the  DMA  points  that  fell 
outside  of  10  km  from  the  tower.  The  DMA  points  within  10  km  of  the  tower 
were  also  terrain  corrected  using  the  digitized  elevation  data. 

For  points  outside  10  km,  we  used  a  procedure  similar  to  the  one  in  North 
Carolina.  There,  we  found  that  the  DMA  data  were  biased  towards  the  higher 
elevations.  Given  that  the  bias  extended  to  20  km  from  the  tower,  we  assumed  a 
constant  bias  of  7  m  from  20  km  out  to  200  km.  This  led  to  a  constantly  sloping 
residual  of  42  uGal  at  the  top  of  the  tower  and  zero  at  the  base2).  Presumably  the 
DMA  data  around  the  WABG  tower  also  contained  some  outer  zone  terrain  bias 
despite  the  flat  terrain. 

Using  the  USGS  and  DMA  digitized  data  together,  we  computed  mean 
elevations  out  to  40  km  from  the  tower.  We  also  computed  mean  elevations  of  the 
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gravity  data  (£40  km)  and  compared  the  results  (Figure  3).  The  figure  clearly 
shows  a  terrain  bias  beyond  10  km;  the  bias  continues  beyond  the  boundaries  of  the 
figure,  growing  to  as  large  as  5  m.  So,  as  in  the  WTVD  experiment,  the  data  are 
biased  towards  higher  elevations  even  though  in  the  inner  survey  area  (<  10  km) 
we  were  very  careful  to  insure  unbiased  data.  We  corrected  for  the  bias  out  to  40 
km  based  on  the  results  of  the  comparison  between  digital  elevations  and  gravity 
elevations  shown  in  Figure  3.  We  then  assumed  that  the  5  m  bias  at  40  km  is 
constant  out  to  300  km  (the  full  extent  of  the  data). 


RADIUS  (KILOMETERS) 

Fig.  3  Azimuthally  averaged  elevations.  Triangles  are  the 
elevations  of  the  gravity  points.  Boxes  are  the  digitized 
elevations  of  USGS  (<10  km)  and  DMA  (10  •  35  km). 

We  analytically  continued  the  surface  Bouguer  anomalies  using  a 
combination  of  a  summation  of  Fourier-Bessel  series  as  a  reference  field,  and  a 
numerical  integration  technique  for  the  residuals  from  the  reference  field3).  We 
performed  a  four-step  nested  symmetric  Bessel  function  fit.  Residuals  were 
computed  after  each  fit  which  served  as  input  to  the  succeeding  fit  at  progressively 
smaller  distances;  this  allowed  for  the  resolution  of  higher  wave  numbers.  Step  1 
is  a  35  parameter  fit  extending  to  300  km;  step  2  is  a  6  parameter  fit  extending  to 
7  km;  step  3  is  a  3  parameter  fit  extending  to  1.3  km;  and  the  final  step  is  a  2 
parameter  fit  extending  to  0.3  km.  Weights  were  computed  for  residuals  from  this 
reference  field  which  were  then  analytically  continued.  The  terrain  had  been 
removed  prior  to  analytic  continuation,  so  its  effect  was  then  added  to  the 
predicted  values  at  the  various  tower  elevations. 
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TOWER  DATA  AND  PRELIMINARY  RESULTS 


The  tower  gravity  experiment  is  currently  incomplete.  We  have  measured 
gravity  at  five  elevations  on  the  WABG  tower  using  the  LaCoste-Romberg 
gravimeter,  G-152.  The  data  were  collected  in  four  loops  with  a  total  of  only  nine 
observations.  The  tower  elevations  were  determined  to  an  accuracy  of  2  cm  using 
an  EDM.  All  measurements  have  been  made  in  less  than  ideal  conditions  with 
wind  speeds  exceeding  15  km/hr.  We  estimate  the  accuracy  of  the  tower  data  to  be 
in  the  range  20-25  pGal  even  though  the  rms  errors  are  on  the  order  of  10  pGal. 
An  attempt  to  measure  at  a  sixth  elevation,  571  m  above  ground  level,  failed.  At 
571  m,  both  the  galvanometer  and  the  reading  line  on  the  gravimeter  were 
disabled,  for  reasons  that  are  as  yet  unclear.  One  possibility  is  the  presence  of  a 
very  strong  magnetic  field,  although  we  cannot  rule  out  other  effects  such  as  radio 
frequency  interference  (RFI).  The  WABG  preliminary  results  (shown  in  Table  1) 
are  plotted  alongside  those  of  the  WTVD  tower  at  commensurate  elevations  in 
Figure  4.  The  agreement  is  good  and,  with  the  exception  of  the  94  m  level,  the 
two  results  agree  to  within  16  pGal. 


Table  1.  Preliminary  WABG  Analytic  Continuation  Results 


Elevation 

Observed 

Predicted 

Observed-Predicted 

(m  above  ground) 

(mGal) 

(mGal) 

(mGal) 

0.000 

9.445±.009 

9.4341.054 

0.0111.055 

93.845 

9.379±.022 

9.3631.025 

0.0161.033 

194.363 

9.300±.022 

9.3051.019 

-0.0051.029 

292.564 

9.233±.023 

9.2431.020 

-0.0101.030 

388.511 

9.148±.023 

9.1791.023 

-0.0311.033 

493.589 

9.0781.024 

9.1121.026 

-0.0341.035 

SUMMARY  AND 

FUTURE  PLANS 

These  results  leave  us  with  several  tasks  to  perform:  1)  obtain  more  tower 
data  during  better  weather  conditions  (wind  speeds  <15  km/hr);  2)  terrain  correct 
the  gravity  data  out  to  40  km  or  beyond;  3)  continue  error  analysis  and  obtain 
improved  error  estimates  for  analytically  continued  values;  and  4)  resolve  the 
problem  at  the  571  m  elevation  so  that  gravity  data  can  be  collected.  Previous  tests 
have  shown  that  RFI  disables  the  galvanometer  but  has  no  effect  on  the  reading 
line.  In  addition,  mu-metal  shielding  around  the  gravimeter  should  protect  it  from 
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stray  magnetic  fields.  So,  either  the  previous  tests  are  somehow  incomplete,  or 
there  is  some  other  yet  unknown  cause.  Once  all  the  above  tasks  have  been 
completed  to  our  satisfaction  we  should  be  able  to  present  our  final  results  for  the 
WABG  tower  and  for  the  topic  of  tower  gravity  in  general. 


HEIGHT  ABOVE  GROUND  (METERS) 

Fig.  4  Observed  minus  model  results  for  the  two  tower 
experiments  and  their  associated  errors.  The  boxes  are  final 
WTVD  results  and  the  diamonds  are  preliminary  WABG  results. 
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ABSTRACT 

In  November  1991,  the  Phillips  Laboratory  (PL),  Geophysics  Directorate  in  cooperation 
with  Purdue  University,  conducted  a  Global  Positioning  System  (GPS),  survey  near  Inverness, 
MS.  The  survey  was  carried  out  around  the  610  m  WABG-TV  transmitting  tower.  The  data  were 
used  in  the  ongoing  program  by  PL  to  study  possible  departures  from  Newtonian  gravity.  Data 
were  collected  within  an  8  km  radius  of  the  tower.  Inside  2  km,  half  of  the  six  GPS  receivers  were 
unable  to  lock  onto  incoming  GPS  signals.  We  believe  that  some  form  of  tower  transmissions, 
possibly  microwave,  interfered  with  the  signal  transmitted  by  the  GPS  satellites. 
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INTRODUCTION 


In  November  1991,  the  Geophysics  Directorate  of  the  Phillips  Laboratory  in  collaboration 
with  Purdue  University  performed  a  GPS  survey  near  the  town  of  Inverness,  MS.  The  survey 
was  designed  to  establish  accurate  geodetic  positions  where  gravity  measurements  could 
subsequently  be  made.  At  the  center  of  the  survey  area  is  the  610  m  WABG-TV  television 
transmitting  tower  which  is  the  main  focus  of  the  ongoing  non-Newtonian  gravity  program  at  PL 
(Jekeli  et  ah,  1990,  Romaides  et  al .,  1991).  The  reason  for  using  GPS  is  that  it  provides  one  of 
the  most  accurate  methods  of  position  determination  for  the  surface  gravity  measurements.  Gravity 
measurements  are  also  being  carried  out  on  this  tower  as  well  as  on  the  ground.  These  ground 
data,  along  with  previously  collected  data  (obtained  from  the  Defense  Mapping  Agency  Gravity 
Library),  will  be  used  to  analytically  predict  what  gravity  should  be  at  various  elevations  on  the 
tower  in  the  absence  of  non-Newtonian  forces.  The  predicted  values  will  then  be  compared  to  the 
tower  measurements.  Any  error  in  the  latitude,  longitude,  or  elevation  of  the  surface  gravity  also 
adds  to  the  error  of  the  predicted  values.  For  the  latitude  and  longitude  the  error  is  on  the  order  of 
2  pGal/m,  and  for  the  elevation  3  pGal/cm  (1  pGal  =  10*6  cm  s'2).  Given  an  allowed  prediction 
error  of  25  pGal,  the  latitudes  and  longitudes  must  be  known  to  an  accuracy  of  20  -  50  cm,  and  the 
elevation  2  -  5  cm  depending  on  the  relative  distance  to  the  tower.  GPS  provides  a  convenient  and 
rapid  way  fen:  the  survey  to  be  completed  given  the  unpredictability  of  the  weather  in  that  area.  The 
survey  consisted  of  134  points  extending  to  a  radius  of  8  km  from  the  tower,  but  for  the  purposes 
of  this  paper  we  shall  focus  on  the  inner  2  km  of  the  survey.  It  is  in  this  area  that  we  believe  the 
television  signal  interfered  with  the  reception  of  the  GPS  signal,  severely  hindering  our  survey 
efforts  in  that  region,  and  leading  to  delays  in  the  project  completion. 

THE  INITIAL  GPS  SURVEY 

The  survey  was  divided  into  concentric  rings  where  the  inner  2  km  radii  were:  150,  300, 
450,  600,  800,  1050,  1400,  and  1900  m  from  the  tower  (rings  A  through  H  respectively).  Each 
ring  contained  10  points  with  every  other  ring  offset  by  18  degrees  (Figure  1).  We  planned  the 
survey  employing  a  six  GPS  receiver  configuration.  The  receivers  were  Ashtech  Inc.  dual- 
frequency  GPS  receivers,  two  M-XII  models  and  four  LD-XH  models.  The  two  M-XII  models 
were  our  own  and  the  four  LD-XII  were  rented  from  Ashtech.  We  labeled  GPS  receivers  1, 2, 3, 
1  A,  2A,  and  3A  strictly  for  identification  purposes.  We  had  two  survey  teams  each  responsible  for 
three  receivers;  one  team  had  receivers  1, 2,  and  3  and  the  other  1A,  2A,  and  3A.  At  the  time  of 
the  survey  there  were  favorable  constellations  during  the  day  that  permitted  five  observing  sessions 
of  1  hour  and  45  minutes  with  30  minutes  allowed  between  setups.  During  the  sessions  there  were 
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four  satellites  available  most  of  the  time  with  five  (and  possibly  six)  satellites  available  for  about 
half  of  the  time  in  each  session.  The  survey  was  designed  to  form  closed  loops  of  six-sided 
polygons  with  each  receiver  placed  at  a  vertex.  The  plan  was  to  do  two  concentric  rings  a  day 
working  outward  from  the  tower  and  repeating  the  outer  of  the  two  rings  on  each  succeeding  day. 
In  each  session  three  receivers  would  be  set  up  on  the  inner  ring  and  three  on  the  outer  ring.  At  the 
end  of  each  session,  four  of  the  six  receivers  would  be  moved  to  new  locations  with  the  other  two 
remaining  fixed  to  form  the  next  polygon.  In  this  way  five  sessions  would  allow  for  the 
completion  of  two  rings;  also  the  first  two  points  observed  in  session  1  would  be  the  same  as  the 
last  two  points  observed  in  session  5.  (See  Figure  2  for  receiver  configuration  and  session 
observations.)  Because  of  some  practical  considerations  involving  accessibility,  the  placement  of 
the  points  was  not  exactly  as  planned  (see  Figure  3). 

On  13  November  1991  (Julian  day  317)  we  began  our  GPS  survey  setting  up  on  the  A  and 
B  rings  (ISO  and  300  m)  with  our  six  receivers.  After  about  10  minutes  we  noticed  that  only  two 
of  the  six  receivers  had  locked  onto  any  satellites,  the  other  four  were  in  a  search  mode  with  none 
of  them  able  to  lock.  (From  this  point  on,  we  will  refer  to  the  receiver  and  antenna 
interchangeably  even  though  it  is  the  antenna  that  locks  onto  the  GPS  satellite  signal  and  not  the 
receiver  itself.)  The  only  two  that  locked  were  at  the  points  designated  B8  and  B9  (See  Figure  2), 
and  even  those  two  did  not  immediately  lock  onto  the  satellite  signals  as  was  usually  the  case. 
Even  after  an  hour  the  situation  did  not  change;  the  receivers  on  points  B8  and  B9  continued  to 
collect  data  while  the  other  four  obtained  no  locks.  Since  there  were  no  obstructions  to  hinder 
incoming  satellite  signals,  the  implication  was  that  something  associated  with  the  television 
transmitter  was  creating  interference.  In  an  attempt  to  determine  the  extent  of  the  problem,  receiver 
3  (selected  at  random)  was  placed  at  several  selected  points  to  ascertain  whether  a  signal  could  be 
detected  anywhere  in  the  vicinity  of  the  tower.  When  this  receiver  was  located  on  the  C  and  D 
rings  located  450  m  and  600  m  respectively  from  the  tower,  no  satellite  locks  were  obtained,  just 
as  in  the  case  of  the  A  and  B  rings.  However  in  the  E  ring  (800  m  from  the  tower),  the  receiver 
was  able  to  lock  onto  four  satellites  in  one  position  (located  SSW  from  the  tower)  but  not  in 
another  position  located  NW  of  the  tower.  Finally  at  the  F  ring,  1050  m  from  the  tower,  the 
receiver  was  able  to  lock  onto  all  available  satellites  in  all  the  azimuthal  orientations  that  were 
tested.  Due  to  the  stringent  time  constraints  we  were  working  under,  a  thorough  analysis  of  the 
interference  problem  was  not  possible  at  that  time.  So  based  on  the  preliminary  analysis  we 
deemed  it  best  to  begin  the  survey  at  the  F  and  G  rings,  and  to  then  try  to  obtain  positions  and 
elevations  for  the  inner  rings  using  some  other  method  such  as  trigonometric  leveling. 

On  14  November  1991  (day  318)  the  receivers  were  placed  on  the  F  and  G  rings  in  the 
configuration  called  for  in  the  survey  plan  (Figure  4).  In  the  first  session  the  six  receivers  were  set 
on  points  GO,  Gl,  G2,  F9,  F0,  and  FI.  Problems  again  surfaced  with  only  three  of  the  six 
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receivers  being  able  to  lock  onto  the  available  satellites.  After  five  complete  observing  sessions  on 
day  318,  we  had  successfully  collected  data  for  seven  out  of  ten  points  in  the  F  ring  and  five  out  of 
ten  points  in  the  G  ring.  On  the  next  day  (day  319)  our  receivers  were  set  up  on  the  G  and  H  rings 
beginning  with  points  HO,  HI,  H2,  GO,  Gl,  and  G2.  This  time  four  of  the  six  receivers  were  able 
to  lock  onto  available  satellites.  As  day  319  progressed  the  situation  improved,  and  by  the  end  of 
the  day  we  had  obtained  data  for  six  out  of  ten  points  in  the  G  ring  (one  more  than  the  previous 
day)  and  eight  out  of  ten  points  in  die  H  ring.  Finally  on  day  320  we  collected  data  for  nine  out  of 
ten  points  in  the  H  ring  and  all  ten  points  in  the  I  ring  (2400  m).  Interestingly,  for  the  H  ring  we 
collected  no  data  for  H2  on  day  319  but  were  successful  on  the  following  day.  Similarly  data  were 
collected  for  H4  on  day  319  but  not  on  day  320.  The  situation  was  different  for  H2  and  H4.  In 
both  of  these  cases  the  receivers  were  able  to  lock  onto  only  one  or  two  satellites  even  though  in 
each  of  the  two  observing  sessions  there  were  at  least  four  and  possibly  five  satellites  available. 
The  H  ring,  situated  1.9  km  from  the  tower,  was  the  furthest  ring  from  where  any  of  our  receivers 
had  trouble  locking  onto  satellite  signals.  The  remainder  of  the  survey  out  to  8  km  was  completed 
without  problems.  Also  where  possible  we  surveyed  points  in  the  inner  rings  using  an  Electronic 
Distance  Meter  (EDM).  We  were  able  to  position  33  points  in  this  manner.  The  EDM 
measurements  were  facilitated  by  the  generally  flat  and  featureless  terrain.  With  the  exception  of  a 
few  trees  there  were  no  significant  obstructions  that  could  account  for  the  reception  problems  that 
we  encountered. 

PRELIMINARY  ANALYSIS 

After  the  initial  GPS  survey,  we  decided  to  do  a  more  rigorous  analysis  of  this  interference 
problem.  We  computed  WGS84  positions  for  the  GPS  surveyed  points  as  well  as  for  the  points 
on  the  inner  rings  that  we  had  measured  using  trigonometric  leveling.  For  the  points  where 
obtaining  GPS  positions  were  unsuccessful,  we  used  a  U.S.  Geological  Survey  topographic  map 
to  obtain  approximate  positions.  Figure  3  shows  the  result  of  this  exercise.  The  unmarked  points 
are  those  with  successfully  obtained  GPS  positions.  Those  identified  in  italics  are  points  where  we 
were  unable  to  obtain  satellite  locks  at  any  time,  and  the  remainder  are  points  where  the  receivers 
were  able  to  lock  at  one  time  but  unable  to  lock  at  another  time.  As  can  be  seen  from  Figure  5,  die 
problem  areas  do  not  appear  to  be  random,  rather  there  appear  to  be  preferred  directions.  First, 
there  is  only  one  point  west  of  the  tower  (beyond  the  E  ring),  G7,  where  we  were  unable  to  obtain 
a  GPS  measurement,  most  of  the  problems  occurring  to  the  east.  But  even  in  the  east  there 
appeared  to  be  preferred  directions.  The  largest  areas  of  difficulty  seemed  to  be  concentrated  in 
three  azimuthal  regions:  NE,  ESE,  and  SSE. 
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We  examined  the  GPS  satellite  constellations  for  the  points  Gl,  G3,  G5,  and  G7.  These 
were  the  points  that  were  observed  twice  but  on  both  occasions  the  receivers  were  unable  to  obtain 
satellite  locks.  Even  though  the  observing  sessions  were  only  separated  by  a  day,  all  four  points 
were  observed  at  different  times  during  the  second  day  with  a  somewhat  different  satellite 
constellation.  However  there  did  not  appear  to  be  any  differences  in  the  satellite  constellations  that 
could  give  rise  to  any  reception  problem.  Furthermore,  stations  where  data  were  not  collected 
were  grouped  with  others  in  the  same  observing  session  where  GPS  signals  were  successfully 
obtained.  This  conflicting  behavior  of  the  different  stations  was  unexpected  since  stations  in  any 
given  observing  session  were  seeing  essentially  the  same  satellite  constellation  even  though  they 
were  separated  by  a  small  distance.  A  similar  exercise  was  attempted  with  the  points  where  data 
had  been  collected  on  one  session  but  not  on  another  session.  Of  particular  interest  were  the  points 
G4  and  GO.  At  G4,  data  were  successfully  collected  on  day  319  but  not  on  day  318.  Examination 
of  the  satellite  sky  plots  for  G4  on  days  318  and  319  again  revealed  no  clues  to  the  problem.  Due 
to  the  similarity  of  the  satellite  configurations,  the  only  other  possibility  was  a  variation  of  the 
television  signal  at  different  times  of  the  day.  Another  example  is  GO,  where  data  were 
successfully  collected  at  the  end  of  day  318,  1720  -  1928  hrs,  but  not  at  the  start  of  the  session, 
0810  -  1005  hrs.  Again,  the  two  satellite  sky  plots  for  the  sessions  1  and  5  respectively,  revealed 
nothing  significant.  The  only  conclusive  feature  of  this  puzzle  was  that  west  of  the  tower  there 
was  only  one  station  that  caused  a  problem;  the  preponderance  of  the  difficulties  were  to  the  east. 
It  is  interesting  to  note  that  in  addition  to  the  610  m  television  transmitter,  there  is  another  smaller 
television  tower  (400  m)  located  about  500  m  northwest.  The  region  northwest  of  the  tower  is  one 
of  only  two  areas,  the  other  being  the  southwest,  where  we  had  no  problem  receiving  satellite 
signals.  Also  during  the  first  attempts  at  observing  on  the  A  and  B  rings,  the  only  two  points 
where  data  were  collected  were  the  points  B8  and  B9,  both  located  in  the  northwest  quadrant.  One 
final  note,  every  day  we  observed  the  rings  in  a  clockwise  direction  starting  in  the  north  where  the 
first  observations  were  made  between  0730  and  0800  hrs,  and  ending  in  the  north  with  final 
observations  made  between  1830  and  1930  hrs.  This  would  mean  that  all  afternoon  observations 
were  made  in  the  west  and  all  morning  observations  were  made  in  the  cast.  Once  again, 
examination  of  the  satellite  sky  plots  for  the  morning  and  afternoon  of  day  318  did  not  reveal  any 
good  reason  for  the  preferred  directions,  so  again  we  suspected  variations  in  the  transmitting 
signal. 

THE  FOLLOW-UP  SURVEY 

As  a  result  of  the  problems  we  encountered  with  GPS,  we  initially  fell  26  points  short  of 
our  objective  of  134  points.  Hence  in  March  1992  we  returned  to  Inverness,  MS  both  to  finish  our 
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survey  and  to  more  clearly  define  some  of  the  problems  we  had  been  having.  We  also  began  our 
gravity  survey  at  the  known  GPS  points.  This  time  we  had  only  our  own  two  M-XII  receivers 
labeled  2A  and  3.  Without  the  luxury  of  six  receivers  the  only  plan  was  to  set  up  one  receiver  on  a 
previously  surveyed  point  and  to  set  up  the  other  receiver  on  an  unknown  point.  Then  as  a  check 
the  receiver  on  the  known  point  would  be  moved  to  another  known  point  In  this  way  coordinates 
for  the  unknown  point  could  be  computed  from  two  known  points  thereby  providing  a  check  on 
the  accuracy.  The  key  to  any  experimental  anomaly  is  repeatability,  and  so  the  question  was 
whether  we  could  repeat  the  observational  problems  we  had  encountered  the  previous  November. 
We  began  our  survey  on  12  March  1992  (day  072)  and  at  first  data  acquisition  went  smoothly  for 
the  most  pan.  In  the  month  of  March  the  satellite  constellation  was  favorable  in  the  morning  and 
from  the  middle  of  the  afternoon  in  to  the  evening.  We  used  the  small  window  where  there  was  a 
dearth  of  available  satellites  (1200  - 1500  hrs)  to  continue  our  testing.  Using  receiver  3  we  tested 
several  areas,  including  those  that  had  previously  proven  problematical,  but  at  no  location  were  we 
unable  to  lock  onto  available  satellites.  At  one  point,  we  even  placed  receiver  3  at  a  distance  of  3  m 
from  the  tower  and  it  was  able  to  lock  onto  all  available  satellites.  It  would  appear  that  more 
questions  were  raised  than  were  answered,  but  as  the  survey  progressed  to  day  073  similar 
problems  began  to  surface.  This  time  with  only  two  receivers  and  no  rigid  time  constraints  we 
finally  began  to  unravel  the  mystery. 

It  turned  out  that  only  one  of  our  two  receivers,  2A,  was  experiencing  the  interference 
problems;  the  other  receiver,  3,  had  no  problem  locking  onto  satellites  in  any  direction  or  at  any 
distance  from  the  tower.  It  would  appear  the  reason  receiver  3  had  problems  on  day  317  during 
the  initial  tests  but  experienced  no  such  difficulties  in  March  1992,  is  simply  that  we  did  not  allow 
the  receiver  enough  time  to  acquire  locks  during  the  early  tests.  We  interchanged  the  antennas  and 
the  receivers  and  the  problem  tracked  the  antenna.  This  was  something  we  had  not  expected.  The 
thought  of  a  malfunctioning  antenna  (or  receiver)  had  occurred  to  us  but  it  was  surprising  that  an 
antenna  would  function  correctly  at  distances  greater  than  2  km  from  a  television  transmitter,  yet 
experience  difficulties  at  closer  distances.  The  antennas  for  all  six  receivers  were  essentially 
identical;  if  one  or  more  could  receive  GPS  signals  in  the  presence  of  a  television  transmitter  it 
seemed  unlikely  that  the  others  would  behave  differently.  Yet  this  was  exactly  what  was 
happening,  although  the  solution  was  still  not  that  simple.  Even  though  receiver  2A  was 
experiencing  difficulties  locking  onto  satellites,  the  problem  was  still  directional  to  a  certain  extent. 
In  the  E  ring  for  the  points  we  tested,  receiver  2A  was  able  to  lock  onto  satellites  for  points  E5  and 
E6  but  not  for  points  El,  E7,  E8,  and  E9.  In  the  F  ring  we  tested  all  but  one  point  Satellite  locks 
were  obtained  for  data  points  FI,  F3,  F5,  F6,  F7,  and  F9  but  not  for  points  F0,  F2,  and  F4  (see 
Figure  5).  So  once  again  there  appeared  to  be  some  directional  variations  in  the  television  signal 
that  were  somehow  interfering  with  the  GPS  signals. 


132 


Armed  with  this  information  we  were  able  to  complete  most  of  the  required  survey.  This 
was  done  simply  by  planning  the  survey  so  that  receiver  3  was  always  closest  to  the  tower  or  at  the 
problem  areas  (e.g.  F2  or  G5),  and  2A  was  further  away  from  the  tower  and  thus,  from  any 
interference.  Unfortunately  three  of  the  GPS  stations  previously  surveyed  had  beet  disturbed  or 
destroyed  before  we  had  had  the  opportunity  to  obtain  gravity  measurements,  so  some  of  these 
needed  to  be  resurveyed  as  well.  Also  with  only  two  receivers,  achieving  the  same  level  of 
redundancy  and  checking  became  more  time  consuming  than  had  been  anticipated.  In  addition  to 
the  GPS  surveying,  a  great  deal  of  time  was  spent  completing  the  majority  of  the  necessary  gravity 
measurements.  The  heavy  work  schedule  meant  that  by  the  end  of  this  trip  we  were  still  about  12 
points  short  of  our  final  objective  thus  necessitating  a  third  trip.  Prior  to  returning  to  Inverness  it 
occurred  to  us  that  perhaps  RF  interference  was  affecting  one  of  our  receivers.  It  was  possible  the 
shielding  around  the  antenna  cable  for  receiver  3  was  superior  to  that  of  receiver  2A.  We 
constructed  some  heavy-duty  (RG-214)  cables  to  be  used  in  place  of  the  existing  antenna  cables. 
We  tested  all  cables  before  leaving  to  ensure  they  were  functioning  properly. 

We  returned  for  a  third  time  to  Mississippi  on  5  April  to  complete  the  GPS  survey.  This 
time  in  addition  to  our  two  M-XII  receivers  we  rented  two  LD-XII  receivers  from  Ashtech  Inc.  As 
it  turned  out,  one  of  the  GPS  receivers  (receiver  2)  was  the  same  one  we  had  rented  in  November 
1991.  The  other  receiver  was  totally  new  to  us.  This  provided  a  unique  opportunity  to  do  some 
additional  testing  on  the  antennas  in  and  around  the  tower.  When  we  began  the  survey  on  6  April 
(day  097)  it  became  obvious  that  the  only  receiver  having  any  problems  locking  onto  satellite 
signals  was  2A.  We  then  decided  to  test  for  RF  interference  from  the  tower  using  the  heavy-duty 
cables  that  we  had  constructed  to  replace  the  RG-58  cables  we  were  previously  using.  We  set  up 
all  four  receivers  at  E8  employing  the  RG-214  cable  on  2A.  The  results  were  the  same,  the  2A 
receiver  was  unable  to  lock  onto  any  satellites  while  the  remaining  three  receivers  had  no  such 
difficulty.  This  was  the  last  test  we  did  with  the  GPS  antennas.  By  this  time  plowing  was  well 
under  way  in  the  fields  surrounding  the  tower,  which  destroyed  several  additional  GPS  points  and 
brought  the  total  number  of  points  requiring  surveying  to  27.  We  proceeded  to  collect  the  final  23 
GPS  points  (five  less  than  we  had  planned)  and  completed  the  gravity  survey  without  further 
problems. 

THE  FINAL  ANALYSIS 

In  the  end,  129  out  of  a  proposed  134  GPS  points  were  measured  in  and  around  the 
WABG-TV  tower  in  Inverness,  MS.  The  five  remaining  points  were  never  surveyed  due  to  the 
inaccessibility  of  the  areas.  The  data  were  obtained  on  three  separate  trips:  November  1991,  and 
in  March  and  April  1992.  Given  that  fact,  the  resolution  to  the  interference  problem  was  not 
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intuitively  obvious  at  the  the  time  of  the  survey  despite  several  pieces  of  evidence  to  that  effect 
Using  all  the  data  at  our  disposal  some  missing  steps  can  be  retraced,  and  some  definitive 
conclusions  can  be  drawn. 

Looking  back  to  November  1991  (day  318)  we  can  make  a  few  conclusive  statements.  We 
now  know  that  receivers  1, 2,  and  3  were  placed  on  the  F  ring  and  receivers  1A,  2A,  and  3A  were 
placed  on  the  G  ring.  Examination  of  the  survey  plan  for  the  F  ring  indicates  that  only  one  of  the 
three  receivers  was  unable  to  lock  onto  the  GPS  signals.  The  survey  plan  for  the  F  ring  was  as 
follows: 

Session  Stations 

1  F9,  FO,  FI 

2  FI,  F2,  F3 

3  F3,  F4,  F5 

4  F5,  F6,  F7 

5  F7,  F8,  F9. 

Given  the  clockwise  direction  we  were  working,  it  is  easy  to  see  how  the  one  problem  receiver 
could  have  been  placed  on  FO,  F2,  and  F4  during  sessions  1,  2,  and  3.  Knowing  that  receiver  3 
was  operating  correctly,  and  having  tested  receiver  2  in  April,  we  concluded  that  receiver  1  was  the 
one  that  was  unable  to  lock  onto  the  GPS  signals  at  various  azimuthal  directions  in  the  tower 
vicinity.  There  were  no  problems  with  any  of  die  receivers  in  sessions  4  and  5.  We  also  know 
that  receivers  1, 2,  and  3  were  placed  on  B7,  B8,  and  B9  on  day  317  when  these  problems  were 
first  encountered.  The  above  conclusions  are  also  consistent  with  day  317  of  the  survey  when  two 
receivers  in  the  B  ring  (undoubtedly  2  and  3)  were  able  to  lock  onto  signals  and  the  third  was  not 
For  the  G  ring,  the  situation  is  less  clear.  The  survey  plan  for  the  G  ring  on  day  318  is  given 
below: 


Session  Stations 

1  GO,  Gl,  G2 

2  G2,  G3,  G4 

3  G4,  G5,  G6 

4  G6,  G7,  G8 

5  G8,  G9,  GO. 

From  past  experience  we  knew  receiver  2A  was  having  problems.  What  we  do  not  know  for  sure 
is  which  of  the  other  two,  1A  or  3A,  experienced  similar  difficulties.  On  day  318  we  failed  to 


134 


collect  data  for  points  GO  (session  1),  Gl,  G3,  G4,  GS,  and  G7.  So  in  the  first  three  sessions  two 
out  of  three  receivers  were  not  collecting  data,  while  in  session  4  the  problem  was  limited  to  one 
receiver.  There  are  two  puzzling  aspects  of  this  conclusion.  Firstly,  for  station  GO,  data  were 
successfully  collected  at  the  end  of  the  day  (session  5)  but  not  at  the  beginning  (session  1). 
Secondly,  on  the  following  day,  day  319,  only  one  receiver  in  each  of  the  first  three  sessions  failed 
to  collect  data  whereas  on  day  318  it  had  been  two  receivers  that  had  failed.  The  one  problematic 
receiver  on  day  319  was  obviously  receiver  2A.  Unfortunately  we  did  not  begin  keeping  track  of 
where  the  receivers  were  placed  until  day  320.  If  one  assumes  some  directional  nature  of 
interference,  it  is  possible  to  construct  a  scenario  where  the  receivers  are  placed  on  different 
stations  on  days  318  and  319  thus  causing  the  inconsistency  between  the  two  days  of  observation. 
Once  again  we  were  faced  with  the  conclusion  that  either  receiver  1 A  or  3A  was  having  difficulty 
locking  onto  GPS  signals  at  a  few  azimuthal  orientations.  One  final  note  on  these  receivers:  on 
day  317  all  three  were  placed  on  the  A  ring  (A7,  A8,  and  A9),  and  after  an  hour  none  of  them  had 
collected  any  data.  So  it  would  appear  all  three  receivers  had  some  degree  of  difficulty  in  receiving 
satellite  signals. 

On  day  319  receivers  1,  2,  and  3  were  placed  on  the  H  ring  beginning  with  HO,  HI,  and 
H2  and  ending  with  H8,  H9,  and  HO.  The  results  of  the  data  collection  were:  no  data  at  H2  and  at 
the  second  observation  of  HO.  Once  again  we  believed  the  problem  was  with  receiver  1.  Since  we 
were  solving  for  latitude,  longitude  and  elevation,  a  minimum  of  four  GPS  satellites  locks  were 
required  before  the  receiver  would  begin  collecting  data.  As  previously  stated,  for  the  H  ring,  and 
to  a  certain  extent  the  G  ring,  the  receiver  was  able  to  lock  onto  one  or  two  of  the  four  or  five 
satellites  available  which  was  not  enough  to  initiate  data  collection.  This  supported  the  theory  that 
at  greater  distances  from  the  tower  the  problems  began  to  subside. 

What  could  cause  three  out  of  six  GPS  receivers  to  be  unable  to  lock  onto  signals  from  the 
GPS  satellites  in  the  vicinity  of  a  television  transmitter  with  no  visible  obstructions?  The  WABG- 
TV  tower  is  610  m  tall  and  transmits  at  100  kW  with  a  frequency  of  88  MHz.  According  to  the 
engineers  at  the  TV  station,  the  transmitter  is  a  circularly  polarized  antenna  transmitting  horizontally 
and  vertically  downward  with  no  measurable  decrease  in  the  signal  strength  up  to  a  distance  of  80 
km.  This  means  that  there  is  no  variation  of  the  television  signal  in  any  azimuthal  orientation.  The 
only  component  of  the  tower  that  could  have  any  directional  effect  are  the  five  4  watt  microwave 
transmitters  located  between  50  and  100  m  above  ground  level.  Three  of  these  microwave  dishes 
transmit  to  Greenville,  and  the  other  two  to  Greenwood.  Greenville  is  WNW  of  the  tower  and 
Greenwood  is  NE  of  the  tower.  Interestingly,  these  are  close  to  two  of  the  four  directions  where 
we  had  difficulty  in  obtaining  satellite  locks.  Perhaps  this  was  all  coincidence  as  the  microwave 
transmitters  still  did  not  explain  the  problems  to  the  east  and  the  south.  Further,  incoming  GPS 
signals  arrive  from  all  directions  and  not  just  from  overhead  so  even  though  we  were  situated 
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directly  under  the  microwave  beam,  there  should  not  have  been  any  reception  problems.  As 
previously  stated,  there  is  another  transmitting  tower  northwest  of  WABG-TV.  This  is  the  400  m 
tall  WMAO-TV  5  MW  transmitter,  and  is  located  between  points  C9  and  CO.  This  tower  contains 
only  one  microwave  transmitter  and  it  transmits  in  the  direction  of  Jackson  (SSE  of  the  tower).  It 
is  again  interesting  to  note  that  SSE  was  one  of  the  four  problem  directions. 

CONCLUSIONS  AND  RECOMMENDATIONS 

What,  if  any,  conclusions  can  we  draw  from  our  experience  in  Inverness,  MS?  The  most 
obvious  conclusion  is  that  GPS  measurements  cannot  routinely  be  made  in  the  direct  vicinity  (<2 
km)  of  television  and  possibly  radio  transmitting  towers.  This,  however,  may  be  a  luxury  that  a 
surveyor  cannot  afford  if  a  GPS  survey  needs  to  be  conducted  at  an  airport  for  instance,  with 
numerous  transmitting  towers.  It  does  seem  possible  that  some  microwave  signals  emanating  from 
these  towers  create  signal  reception  difficulties  for  some  GPS  receivers.  In  three  out  of  the  four 
directions  where  we  experienced  signal  reception  difficulties  there  are  microwave  transmissions 
present.  The  problems  we  encountered  due  east  remain  largely  unexplained.  We  cannot  devise  a 
mechanism  to  explain  how  the  microwave  transmission  interferes  with  the  GPS  signal,  but  it  is 
possible  that  the  signal  to  noise  ratio  is  decreased  making  it  difficult  for  less  sensitive  GPS 
antennas  to  obtain  satellite  locks.  The  GPS  frequencies  are  at  1575.42  MHz  for  LI  and  1227.6 
MHz  for  L2.  It  is  possible  that  various  harmonics  of  the  microwave  transmissions  could  fall  near 
the  GPS  frequencies  thus  causing  some  degree  of  receiver  reception  problems  (Johannessen  et  cd ., 
1990).  The  problems  undoubtedly  will  vary  from  one  transmitter  to  another.  Also  it  has  been 
shown  that  GPS  signals  are  susceptible  to  jamming  whether  it  is  unintentional  or  deliberate 
(Johannessen,  1992;  Owen,  1992). 

Given  the  fact  that  the  GPS  receiver  manufacturers  contend  that  television  or  radio 
transmitters  should  pose  no  problems,  how  should  one  proceed?  If  a  situation  necessitates  a  GPS 
survey  near  a  transmitting  tower,  the  only  prudent  course  of  action  would  be  to  test  all  receivers 
thoroughly  prior  to  commencement  of  the  survey.  The  test  should  be  conducted  at  various  times  of 
the  day  and  in  all  azimuthal  orientations  extending  out  to  4  km.  Had  we  performed  these  tests  and 
been  aware  of  all  these  problems  in  November  1991,  it  is  quite  possible  we  would  have  been  able 
to  complete  the  entire  GPS  survey  (134  points)  without  any  additional  trips.  It  is  our  belief  that  for 
every  type  of  transmitting  tower  some  type  of  interference  is  present  Therefore  any  GPS  surveys 
that  are  performed  in  their  vicinity  will  experience  problems  similar  to  the  ones  we  encountered.  It 
is  up  to  the  surveyor  to  take  precautionary  steps  in  order  to  avoid  potential  problems,  or  using  a 
variation  to  the  old  salesperson  adage:  "let  the  surveyor  beware." 
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FIGURE  CAPTIONS 


Figure  1.  The  original  survey  plan  of  the  GPS  points  within  a  2  km  radius  of  the  tower. 

Figure  2.  The  original  plan  for  the  inner  three  rings  (A,  B,  and  C)  with  each  six-sidded 

polygon  representing  one  session  of  observation.  The  figure  depicts  two  days  of  observations 
consisting  of  five  1  hour  and  45  minute  sessions. 

Figure  3.  The  final  survey  plan  of  die  GPS  pants  within  2  km  of  the  tower. 

Figure  4.  The  survey  plan  for  the  F,  G,  and  H  rings  on  Julian  days  318  and  319  prior  to  the 

discovery  of  the  interference  problem. 

Figure  5.  The  survey  points  of  the  F,  G,  and  H  rings.  Points  without  designation  are  those 
where  GPS  data  was  successfully  collected.  The  points  labeled  in  italic  are  those  where  GPS  data 
was  not  originally  collected  despite  more  than  one  attempt  at  obtaining  data.  The  remaining  points 
are  ones  where  data  was  collected  on  one  day  but  not  on  another  day.  The  solid  lines  represent  the 
approximate  regions  where  reception  difficulties  were  prevalent 
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